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Foreword 


If you’re reading this book, you’re probably already enrolled in an introductory university course 
in Mechanical Engineering. The primary goals of this textbook are, to provide you, the student, 
with: 

1. An understanding of what Mechanical Engineering is and to a lesser extent what it is not 

2. Some useful tools that will stay with you throughout your engineering education and career 

3. A brief but significant introduction to the major topics of Mechanical Engineering and 
enough understanding of these topics so that you can relate them to each other 

4. A sense of common sense 

The challenge is to accomplish these objectives without diluting the effort so much that you can’t 
retain anything. 

In regards to item 2 above, many of my university courses I remember nothing about, even if I 
use the information I learned therein. In others I remember “factoids” that I still use. One goal of 
this textbook is to provide you with a set of useful factoids so that even of you don’t remember any 
specific words or figures from this text, and don’t even remember where you learned these factoids, 
you still retain them and apply them when appropriate. 

In regards to item 3 above, in particular the relationships between topics, this is one area where I 
feel engineering faculty (myself included) do not do a very good job. Time and again, I find that 
students learn something in class A, and this in formation is used with different terminology or in a 
different context in class B, but the students don’t realize they already know the material and can 
exploit that knowledge. As the old saying goes, “We get too soon old and too late smart...” 
Everyone says to themselves at some point in their education, “oh... that’s so easy... why didn’t the 
book [or instructor] just say it that way...” I hope this text will help you to get smarter sooner and 
older later. 

A final and less tangible purpose of this text (item 4 above) is to try to instill you with a sense of 
common sense. Over my 26 years of teaching, I have found that students have become more 
technically skilled and well rounded but have less ability to think and figure out things for 
themselves. I attribute this in large part to the fact that when I was a teenager, cars were relatively 
simple and my friends and I spent hours working on them. When our cars weren’t broken, we 
would sabotage (nowadays “hack” might be a more descriptive term) each others’ cars. The best 
hacks were those that were difficult to diagnose, but trivial to fix once you know what was wrong. 
We learned a lot of common sense working on cars. Today, with electronic controls, cars are very 
difficult to work on or hack. Even with regards to electronics, today the usual solution to a broken 
device is to throw it away and buy a newer device, since the old one is probably nearly obsolete by 
the time it breaks. Of course, common sense per se is probably not teachable, but a sense of common 
sense , that is, to know when it is needed and how to apply it, might be teachable. If I may be allowed 
an immodest moment in this textbook, I would like to give an anecdote about my son Peter. When 
he was not quite 3 years old, like most kids his age had a pair of shoes with lights (actually light- 
emitting diodes or LEDs) that flash as you walk. These shoes work for a few months until the heel 
switch fails (usually in the closed position) so that the LEDs stay on continuously for a day or two 
until the battery goes dead. One morning he noticed that the LEDs in one of his shoes were on 
continuously. He had a puzzled look on his face, but said nothing. Instead, he went to look for his 
other shoe, and after rooting around a bit, found it. He then picked it up, hit it against something 
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and the LEDs flashed as they were supposed to. He then said, holding up the good shoe, “this shoe 
- fixed... [then pointing at the other shoe] that shoe - broken!” I immediately thought, “I wish all 
my students had that much common sense..In my personal experience, about half of engineering 
is common sense as opposed to specific, technical knowledge that needs to be learned from 
coursework. Thus, to the extent that common sense can be taught, a final goal of this text is to try 
to instill this sense of when common sense is needed and even more importantly how to integrate it 
with technical knowledge. The most employable and promotable engineering graduates are the most 
flexible ones, i.e. those that take the attitude, “I think I can handle that” rather than “I can’t handle 
that since no one taught me that specific knowledge.” Students will find at some point in their 
career, and probably in their very first job, that plans and needs change rapidly due to testing failures, 
new demands from the customer, other engineers leaving the company, etc. 

In most engineering programs, retention of incoming first-year students is an important issue; at 
many universities, less than half of first-year engineering students finish an engineering degree. Of 
course, not every incoming student who chooses engineering as his/her major should stay in 
engineering, nor should every student who lacks confidence in the subject drop out, but in all cases 
it is important that incoming students receive a good enough introduction to the subject that they 
make an informed, intelligent choice about whether he/she should continue in engineering. 

Along the thread of retention, I would like to give an anecdote. At Princeton University, in one 
of my first years of teaching, a student in my thermodynamics class came to my office, almost in 
tears, after the first midterm. She did fairly poorly on the exam, and she asked me if I thought she 
belonged in Engineering. (At Princeton thermodynamics was one of the first engineering courses 
that students took). What was particularly distressing to her was that her fellow students had a 
much easier time learning the material than she did. She came from a family of artists, musicians 
and dancers and got little support or encouragement from home for her engineering studies. While 
she had some of the artistic side in her blood, she said that her real love was engineering, but she 
wondered was it a lost cause for her? I told her that I didn’t really know whether she should be an 
engineer, but I would do my best to make sure that she had a good enough experience in 
engineering that she could make an informed choice from a comfortable position, rather than a 
decision made under the cloud of fear of failure. With only a little encouragement from me, she did 
better and better on each subsequent exam and wound up receiving a very respectable grade in the 
class. She went on to graduate from Princeton with honors and earn a Ph.D. in engineering from a 
major Midwestern university. I still consider her one of my most important successes in teaching. 
Thus, a goal of this text is (along with the instructor, fellow students, and infrastructure) is to 
provide a positive first experience in engineering. 

There are also many topics that should be (and in some instructors’ views, must be) covered in an 
introductory engineering textbook but are not covered here because the overriding desire to keep 
the book’s material manageable within the limits of a one-semester course: 

1. History of engineering 

2. Philosophy of engineering 

3. Engineering ethics 

Finally, I offer a few suggestions for faculty using this book: 

1. Syllabus. Appendix A gives an example syllabus for the course. As Dwight Eisenhower 
said, “plans are nothing... planning is everything.” 

2. Projects. I assign small, hands-on design projects for the students, examples of which are 
given in Appendix B. 
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3. 


Demonstrations . Include simple demonstrations of engineering systems — thermoelectrics, 
piston-type internal combustion engines, gas turbine engines, transmissions, ... 

4. Computer graphics. At USC, the introductory Mechanical Engineering course is taught in 
conjunction with a computer graphics laboratory. 
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Chapter 1. What is Mechanical Engineering? 

'The journal of a thousand miles begins with one step. " 

- Lao Zhu 

Definition of Mechanical Engineering 

My favorite definition of Mechanical Engineering is 

If it needs engineering but it doesn't involve electrons, chemical reactions , arrangement of molecules, life forms, isn't 
a structure (building! bridge / dam) and doesn'tfly, a mechanical engineer will take care of it... but 

if it does involve electrons, chemical reactions, arrangement of molecules, life forms, is a structure or does fly, 
mechanical engineers may handle it anyway 

Although every engineering faculty member in every engineering department will claim that 
his/her field is the broadest engineering discipline, in the case of Mechanical Engineering that’s 
actually true because the core material permeates all engineering systems (fluid mechanics, solid 
mechanics, heat transfer, control systems, etc.) 

Mechanical engineering is one of the oldest engineering fields (perhaps Civil Engineering is even 
older) but in the past 20 years has undergone a rather remarkable transformation as a result of a 
number of new technological developments including 

• Computer Aided Design (CAD). The average non-technical person probably thinks that 
mechanical engineers sit in front of a drafting table drawing blueprints for devices having nuts, 
bolts, shafts, gears, bearings, levers, etc. While that image was somewhat true 100 years ago, 
today the drafting board has long since been replaced by CAD software, which enables a part to 
be constructed and tested virtually before any physical object is manufactured. 

• Simulation. CAD allows not only sizing and checking for fit and interferences, but the 
resulting virtual parts are tested structurally, thermally, electrically, aerodynamically, etc. and 
modified as necessary before committing to manufacturing. 

• Sensor and actuators. Nowadays even common consumer products such as automobiles have 
dozens of sensors to measure temperatures, pressures, flow rates, linear and rotational speeds, 
etc. These sensors are used not only to monitor the health and performance of the device, but 
also as inputs to a microcontroller. The microcontroller in turn commands actuators that adjust 
flow rates (e.g. of fuel into an engine), timings (e.g. of spark ignition), positions (e.g. of valves), 
etc. 

• 3D printing. Traditional “subtractive manufacturing” consisted of starting with a block or 
casting of material and removing material by drilling, milling, grinding, etc. The shapes that can 
be created in this way is limited compared to modern “additive manufacturing” or “3D printing” 
in which a structure is built in layers 

• Collaboration with other fields. Historically, a nuts-and-bolts device such as an automobile 
was designed almost exclusively by mechanical engineers. Modern vehicles have vast electrical 
and electronic systems, safety systems (e.g. air bags, seat restraints), specialized batteries (in the 
case of hybrids or electric vehicles), etc., which require design contributions from electrical, 
biomechanical and chemical engineers, respectively. It is essential that a modern mechanical 



engineer be able to understand and accommodate the requirements imposed on the system by 
non-mechanical considerations. 

These radical changes in what mechanical engineers do compared to a relatively short time ago 
makes the field both challenging and exciting. 

Mechanical Engineering curriculum 

In almost any accredited Mechanical Engineering program, the following courses are required: 

• Basic sciences - math, chemistry, physics 

• Breadth or distribution (called “General Education” at USC) 

• Computer graphics and computer aided design 

• Experimental engineering & instrumentation 

• Mechanical design - nuts, bolts, gears, welds 

• Computational methods - converting continuous mathematical equations into discrete 
equations for example 

• Core “engineering science” 

o Dynamics — essentially F = ma applied to many types of systems 
o Strength and properties of materials 
o Fluid mechanics 
o Thermodynamics 
o Heat transfer 
o Control systems 

• Senior “capstone” design project 

Additionally you may participate in non-credit “enrichment” activities such as undergraduate 
research, undergraduate student paper competitions in ASME (American Society of Mechanical 
Engineers, the primary professional society for mechanical engineers, the SAE Formula racecar 
project, etc. 



Figure 1. SAE Formula racecar project at USC (photo: http://www-scf.usc.edu/~scracer/) 
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Examples of industries employing MEs 


Many industries employ mechanical engineers; a few industries and the type of systems MEs 
design are listed below, 
o Automotive 

• Combustion 

• Engines, transmissions 

• Suspensions 

o Aerospace (w/ aerospace engineers) 

• Control systems 

• Heat transfer in turbines 


o 


o 


o 


o 


o 


o 


• Fluid mechanics (internal & external) 
Biomedical (w/ physicians) 

• Biomechanics — prosthesis 

• Flow and transport in vivo 
Computers (w/ computer engineers) 

• Heat transfer 

• Packaging of components & systems 
Construction (w/ civil engineers) 

• Heating, ventilation, air conditioning (HVAC) 

• Stress analysis 

Electrical power generation (w/ electrical engineers) 

• Steam power cycles - heat and work 

• Mechanical design of turbines, generators, ... 
Petrochemicals (w/ chemical, petroleum engineers) 

• Oil drilling - stress, fluid flow, structures 

• Design of refineries - piping, pressure vessels 
Robotics (w/ electrical engineers) 

• Mechanical design of actuators, sensors 

• Stress analysis 
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Chapter 2. Units 

The heart of science is measurement . 

— Erik Brynjolfsson 

All engineered systems require measurements for specifying the size, weight, speed, etc. of 
objects as well as characterizing their performance. Understanding the application of these units is 
the single most important objective of this textbook because it applies to all forms of engineering 
and everything that one does as an engineer. Understanding units is far more than simply being able 
to convert from feet to meters or vice versa; combining and converting units from different sources 
is a challenging topic. For example, if building insulation is specified in units of BTU inches per 
hour per square foot per degree Fahrenheit, how can that be converted to thermal conductivity in 
units of Watts per meter per degree C? Or can it be converted? Are the two units measuring the 
same thing or not? (For example, in a new engine laboratory facility that was being built for me, the 
natural gas flow was insufficient... so I told the contractor I needed a system capable of supplying a 
minimum of 50 cubic feet per minute (cfm) of natural gas at 5 pounds per square inch (psi). His 
response was “what’s the conversion between cfm and psi?” Of course the answer is that there is 
no conversion; cfm is a measure of flow rate and psi a measure of pressure.) Engineers have to 
struggle with these misconceptions every day. 

Engineers in the United States are burdened with two systems of units and measurements: 
(1) the English or USCS (US Customary System) © and (2) the metric or SI (Systeme International 
d’Unites) ©. Either system has a set of base unis , that is, units which are defined based on a 
standard measure such as a certain number of wavelengths of a particular light source. These base 
units include: 

• Length (meters, centimeter, feet, inches); 1 m = 100 cm = 3.281 ft = 39.37 in 

• Mass (lbm, slugs, kilograms); (1 kg = 2.205 lbm) (lbm = “pounds mass”) 

• Time (seconds; the standard abbreviation is “s” not “sec”) 

• Electric current (really electric charge is the base unit, and derived unit is current = 
charge/time) (1 coulomb = charge on 6.241506 x 10 18 electrons) (1 amp = 1 
coulomb/ second) 

Moles are often reported as a fundamental unit, but it is not; it is just a bookkeeping 
convenience to avoid carrying around factors of 10“’ everywhere. The choice of the number of 
particles in a mole of particles is completely arbitrary; by convention Avogadro’s number is defined 
by N a = 6.0221415 x 10 23 , the units being particles/mole (or one could say individuals of any kind, 
not limited just particles). 

Temperature is frequently interpreted as a base unit but again it is not, it is a derived unit , that is, 
one created from combinations of base units. Temperature is essentially a unit of energy divided by 
Boltzman’s constant. Specifically, the average kinetic energy of an ideal gas molecule in a 3- 
dimensional box is 1.5kT, where k is Boltzman’s constant = 1.380622 x 10 ’ J/K (really 
(Joules/molecule)/K). Thus, 1 Kelvin is the temperature at which the kinetic energy of an ideal gas 
(and only an ideal gas, not any other material) molecule is 1.5kT =2.0709 x 10’ ’ J. The ideal gas 

constant (91) with which are you are very familiar is simply Boltzman’s constant multiplied by 
Avogadro’s number, i.e. = kN A = 1.38 x 10’ 2 ’ J/moleculeK * 6.02 x 10 2 ’ molecules / mole = 
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8.314 J/moleK — 1.987 cal/moleK. There’s also another type of gas constant R — 91/(M, where (M 

= molecular mass of the gas; R depends on the type of gas whereas 91 is the “universal” gas constant, 
i.e., the same for any gas. Why does this discussion apply only for an ideal gas? By definition, ideal 
gas particles have only kinetic energy and negligble potential energy due to inter-molecular 
attraction; if there is potential energy, then we need to consider the total internal energy of the 
material (E, units of Joules) which is the sum of the microscopic kinetic and potential energies, in 
which case the temperature for any material (ideal gas or not) is defined as 



(Equation 1) 


where S is the entropy of the material (units J/K) and V is the volume. This intimidating-looking 
definition of temperature, while critical to understanding thermodynamics, will not be needed in this 
course. 


Derived units are units created from combinations of base units; there are an infinite number of 
possible derived units. Some of the more important/common/useful ones are: 

• Area = length 2 ; 640 acres = 1 mile 2 , or 1 acre = 43,560 ft 2 

• Volume = length’; 1 ft 3 = 7.481 gallons = 28,317 cm 3 ; also 1 liter = 1000 cm 3 = 61.02 in 3 

• Velocity = length/time 

• Acceleration = velocity/time = length/time 2 (standard gravitational acceleration on earth = 
g = 32.174 ft/s 2 = 9.806 m/s 2 ) 

• Force = mass * acceleration = mass*length/time 2 

o 1 kg m/s 2 = 1 Newton = 0.2248 pounds force (pounds force is usually abbreviated 
lbf) 

2 2 

• Energy = force x length = mass x length /time - 

o 1 kg m 2 /s 2 = 1 Joule (J) 

o 778 ft lbf = 1 British thermal unit (BTU) 

o 1055 J = 1 BTU 

o 1J = 0.7376 ft lbf 

o 1 calorie = 4.184J 

o 1 dietary calorie = 1000 calories 

• Power (energy/time = mass x length 2 /time’) 

o 1 kg m 2 /s 3 = 1 Watt 
o 746 W = 550 ft lbf/sec = 1 horsepower 

• Heat capacity = J/moleK or J/kgK or J/mole°C or J/kg°C (see note below) 

• Pressure = force/area 

o 1 N/m 2 = 1 Pascal 

o 101325 Pascal = 14.686 lbf/in 2 = 1 standard atmosphere 

• Volts = energy/charge = J/coulomb 

• Capacitance = amps / (volts/s) (1 farad = 1 coul 2 /J) 

• Inductance = volts / (amps/s) (1 Henry = 1 J s 2 / coul 2 ) 

• Resistance = volts/amps (1 ohm = 1 volt/amp = 1 Joule-s / coul 2 ) 
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^ ^ 2 / * 2 

Torque = force x lever arm length = mass x length"/time" — same as energy but one would 
usually report torque in N-m, not Joules, to avoid confusion. 

Radians, degrees, revolutions — these are all dimensionless quantities, but must be converted 
between each other, i.e. 1 revolution = 2 tc radians = 360 degrees. 


By far the biggest problem with USCS units is with mass and force. The problem is that pounds is 
both a unit of mass AND force. These are distinguished by lbm for pounds (mass) and lbf for 
pounds (force). We all know that W = mg where W = weight, m = mass, g = acceleration of gravity. 
So 

1 lbf = 1 lbm x g = 32.174 lbm ft/s 2 (Equation 2) 

Sounds ok, huh? But wait, now we have an extra factor of 32.174 floating around. Is it also true 
that 

1 lbf = 1 lbm ft/s 2 (Equation 3) 

which is analogous to the SI unit statement that 

1 Newton = 1 kg m/s 2 (Equation 4) 


No, 1 lbf cannot equal 1 lbm ft/s 2 because 1 lbf equals 32.174 lbm ft/sec 2 . So what unit of mass 
satisfies the relation 

1 lbf = 1 (mass unit)-ft/s 2 ? (Equation 5) 

This mass unit is called a “slug” believe it or not. By comparison of Equations (2) and (5), 

1 slug = 32.174 lbm = 14.59 kg (Equation 6) 

Often when doing USCS conversions, one uses a units conversion factor called 


32.174 lbm ft 
lbf? 


(Equation 7) 


One can multiply and divide any equation by g c — 1 as many times as necessary to get the units 
correct (an example of “why didn't somebody just say that?”) 


If this seems confusing, it is to me too. That’s why I recommend that even for problems in which 
the givens are in USCS units and where the answer is needed in USCS units, first convert everything 
to SI units, do the problem, then convert back to USCS units. I disagree with some authors who say 
an engineer should be fluent in both systems. The first example below uses the approach of 
converting to SI, do the problem, and convert back to USCS. The second example shows the use of 
USCS units employing g c : 


Example 1 

What is the weight (in lbf) of one gallon of air at 1 atm and 25°C? The molecular weight of air is 
28.97 g/mole = 0.02897 kg/mole. 


Ideal gas law: PV — n$tT 
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(P — pressure, V — volume, n — number of moles, — university gas constant, T — 
temperature) 

Mass of gas (m) = moles x mass/mole = n*M (*M = molecular mass) 

Weight of gas (W) = mg 

Combining these 3 relations: W = PVMg/DTT 




PVTftg 


latm 


101325N/m 

atm 


2 \ 



0.02897kg V 9.81m 
mole / 



8.314J 

moleK 


(25+273)K 


= 0.0440 


mole / \ s 


m 


2 


2 


(N)( m )(kg) 


m 

~2 


(Nm) kg 


m 

Zf 


= 0.0440 


= 0.0440 


mole 


= 0.0440 N °~ 22481bf = 0.00989 lbf = 0.01 lbf 

N 

Note that it’s easy to write down all the formulas and conversions. The tricky part is to 
check to see if you’ve actually gotten all the units right. In this case I converted everything 
to the SI system first, then converted back to USCS units at the very end — which is a pretty 
good strategy for most problems. 

Example 2 

A 3000 pound (3000 lbm) car is moving at a velocity of 88 ft/sec. What is its kinetic energy (KE) in 
ft lbf? What is its kinetic energy in Joules? 


KE = ^(mass)(velocity) 2 =^-(3000 lbm)^88 —j = 1.16 xlO 7 


• 2/2 ^ # 2 / . 2 * 

Now what can we do with lbm ft /sec - ??? The units are (mass)(length) - /(time) - , so it is a unit of 
energy, so at least that part is correct. Dividing by g c , we obtain 


KE = 1.16x10 


7 


lbm ft 2 1 ' 


2 


X 


g 


i in7 lbm ft 
1.16x10 -— 


2 \ 


c \ 


lbf sec 


2 \ 


32.174 lbm ft 


= 3.61 x 10 5 ft lbf 


KE = (3.61 x 10 5 ft lbf) 


1 J 


0.7376 ft lbf 


= 4.89 x 10 5 J 
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Note that if you used 3000 lbf rather than 3000 lbm in the expression for KE, you’d have the wrong 
units — ft lbf"/lbm, which is NOT a unit of energy (or anything else that I know of...) Also note 

* • •• ** 23 / 

that since g c = 1, we COULD multiply by g c rather than divide by g c ; the resulting units (lbm“ ft /lbf 
sec 4 ) is still a unit of energy, but not a very useful one! 


Many difficulties also arise with units of temperature. There are four temperature scales in 
“common” use: Fahrenheit, Rankine, Celsius (or Centigrade) and Kelvin. Note that one speaks of 
“degrees Fahrenheit” and “degrees Celsius” but just “Rankines” or “Kelvins” (without the 
“degrees”). 


T (in units of °F) — T (in units of R) - 459.67 
T (in units of °C) = T (in units of K) - 273.15 
1 K= 1.8 R 

T (in units of°C) = [T (in units of °F) — 32]/1.8, 

T (in units of °F) = 1.8[T (in units of °C)] + 32 
Water freezes at 32°F / 0°C, boils at 212°F / 100°C 


Special note (another example of “that’s so easy, why didn’t somebody just say that?”): when using units 
involving temperature (such as heat capacity, units J/kg°C, or thermal conductivity, units 
Watts/m°C), one can convert the temperature in these quantities these to/from USCS units (e.g. 
heat capacity in BTU/lbm°F or thermal conductivity in BTU/hr ft °F) simply by multiplying or 
dividing by 1.8. You don’t need to add or subtract 32. Why? Because these quantities are really 
derivatives with respect to temperature (heat capacity is the derivative of internal energy with respect 
to temperature) or refer to a temperature gradient (thermal conductivity is the rate of heat transfer 
per unit area by conduction divided by the temperature gradient, dT/dx). When one takes the 
derivative of the constant 32, you get zero. For example, if the temperature changes from 84°C to 
17°C over a distance of 0.5 meter, the temperature gradient is (84-17)/0.5 = 134°C/m. In 
Fahrenheit, the gradient is [(1.8*84 +32) - (1.8*17 + 32)]/0.5 = 241.2°F/m or 241.2/3.281 = 
73.5°F/ft. The important point is that the 32 cancels out when taking the difference. So for the 
purpose of converting between °F and °C in units like heat capacity and thermal conductivity, one can use 1°C = 
1.8°F. That doesn’t mean that one can just skip the + or — 32 whenever one is lazy. 

Also, one often sees thermal conductivity in units of W/m°C or W/mK. How does one convert 
between the two? Do you have to add or subtract 273? And how do you add or subtract 273 when 
the units of thermal conductivity are not degrees? Again, thermal conductivity is heat transfer per 
unit area per unit temperature gradient. This gradient could be expressed in the above example as 
(84°C-17°C)/0.5 m = 134°C/m, or in Kelvin units, [(84 + 273)K— (17 + 273)I<]/0.5 m = 134K/m 
and thus the 273 cancels out. So one can say that 1 W/m°C = 1 W/mK, or 1 J/kg°C = 1 J/kgK. 
And again, that doesn’t mean that one can just skip the + or — 273 (or 460, in USCS units) 
whenever one is lazy. 

Example 3 


BTU inch 


The thermal conductivity of a particular brand of ceramic insulating material is o.5 

ft z hour °F 

(I’m not kidding, these are the units commonly reported in commercial products!) What is the 

thermal conductivity in units of ^ atts ? 

meter °C 
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„ , BTU inch 1055 J ft 3.281 ft hour 1 Watt 1.8°F „ Watt 
ft 2 hour °F BTU 12 inch m 3600 s 1 J/s °C m°C 

Note that the thermal conductivity of air at room temperature is 0.026 Watt/m°C, i.e. about 3 times 
lower than the insulation. So why don’t we use air as an insulator? We’ll discuss that in Chapter 8. 
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Chapter 3. “Engineering scrutiny” 

“Be jour own worst critic , unlessyou prefer that someone else be your worst critic. ” 

-1 dunno, I just made it up. But, it doesn’t sound very original. 

Scrutinizing analytical formulas and results 

I often see analyses that I can tell within 5 seconds must be wrong. I have three tests, which should 
be done in the order listed, for checking and verifying results. These tests will weed out 95% of all 
mistakes. I call these the “smoke test,” “function test,” and “performance test,” by analogy with 
building electronic devices. 

1. Smoke test. In electronics, this corresponds to turning the power switch on and seeing if the 
device smokes or not. If it smokes, you know the device can’t possibly be working right (unless 
you intended for it to smoke.) In analytical engineering terms, this corresponds to checking 
the units. You have no idea how many results people report that can’t be correct because the 
units are wrong (i.e. the result was 6 kilograms, but they were trying to calculate the speed of 
something.) You will catch 90% of your mistakes if you just check the units. For example, 
if I just derived the ideal gas law for the first time and predicted PV = nOt/T you can quickly see 
that the units on the right-hand side of the equation are different from those on the left-hand 
side. There are several additional rules that must be followed: 


Anything inside a square root, cube root, etc. must have units that are a perfect square (e.g. 
m 2 /sec 2 ), cube, etc.) This does not mean that every term inside the square root must be a 
perfect square, only that the combination of all terms must be a perfect square. For 
example, the speed (v) of a frictionless freely falling object in a gravitational field is 


v = ^2gh , where g — acceleration of gravity (units length/time 2 ) and h is the height from 


which the object was dropped (units length). Neither g nor h have units that are a perfect 
square, but when multiplied together the units are (length/time 2 ) (length) = length 2 /time 2 , 
which is a perfect square, and when you take the square root, the units are 


v = yj length 2 /time 2 = length/time as required. 


Anything inside a log, exponent, trigonometric function, etc., must be dimensionless (I don’t 
know how to take the log of 6 kilograms). Again, the individual terms inside the function 
need not all be dimensionless, but the combination must be dimensionless. 

Any two quantities that are added together must have the same units (I can’t add 6 kilograms 
and 19 meters/second. Also, I can add 6 miles per hour and 19 meters per second, but I 
have to convert 6 miles per hour into meters per second, or convert 19 meters per second 
into miles per hour, before adding the terms together.) 


2. function test. In electronics, this corresponds to checking to see if the device does what I 
designed it to do, e.g. that the red light blinks when I flip switch on, the meter reading increases 
when I turn the knob to the right, the bell rings when I push the button, etc. — assuming that 
was what I intended that it do. In analytical terms this corresponds to determining if the result 
gives sensible predictions. Again, there are several rules that must be followed: 
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Determine if the sign (+ or -) of the result is reasonable. For example, if your prediction of 
the absolute temperature of something is —72 Kelvin, you should check your analysis again. 

Determine whether what happens to y as x goes up or down is reasonable or not. For 
example, in the ideal gas law, PV = n$tT: 

o At fixed volume (V) and number of moles of gas (n), as T increases then P increases 
— reasonable 

o At fixed temperature (T) and n, as V increases then P decreases — reasonable 
o Etc. 

Determine what happens in the limit where x goes to special values, e.g. zero, one or infinity 
as appropriate. For example, consider the equation for the temperature as a function of time 
T(t) of an object starting at temperature T } at time t = 0 having surface area A (units m 2 ), 
volume V (units m 3 ), density p (units kg/m') and specific heat C P (units J/kg°C) that is 
suddenly dunked into a fluid at temperature with heat transfer coefficient h (units 
Watts/m 2o C). It can be shown that in this case T(t) is given by 


T(t) = T„ + (T i -TJexip 


hA 


pVC 


t 


p ) 


(Equation 8) 


hA/pVCp has units of (Watts/m 2o C)(m 2 )/(kg/m'')(m'')(J/kg°C) = 1/s, so (hA/pVC P )t is 
dimensionless, thus the formula easily passes the smoke test. But does it make sense? At t — 
0, Tj = 0 as expected. What happens if you charge for a long time? The temperature can 
reach T*, but not overshoot it. In the limit / —> oo, the term exp(-(hA/pVC P )t) goes to zero, 
thus T —> T oo as expected. Other scrutiny checks: if h or A increases, heat can be transferred 

to the object more quickly, thus the time to approach decreases. Also, if p, V or C P 
increases, the “thermal inertia” (resistance to change in temperature) increases, so the time 
required to approach T^ increases. So the formula makes sense. 

• If your formula contains a difference of terms, determine what happens if those 2 terms are 
equal. For example, in the above formula, if = T*,, then the formula becomes simply T(t) 
= for all time. This makes sense because if the bar temperature and fluid temperature are 
the same, then there is no heat transfer to or from the bar and thus its temperature never 
changes. 

3. Performance test. In electronics, this corresponds to determining how fast, how accurate, etc. the 
device is. In analytical terms this corresponds to determining how accurate the result is. This 
means of course you have to compare it to something else that you trust, i.e. an experiment, a 
more sophisticated analysis, someone else’s published result (of course there is no guarantee that 
their result is correct just because it got published, but you need to check it anyway.) For 
example, if I derived the ideal gas law and predicted PV = 7nRT, it passes the smoke and 
function tests with no problem, but it fails the performance test miserably (by a factor of 7). But 
of course the problem is deciding which result to trust as being at least as accurate as your own 
result; this of course is something that cannot be determined in a rigorous way, it requires a 
judgment call based on your experience. 
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Scrutinizing computer solutions 

(This part is beyond what I expect you to know for AME 101 but I include itfor completenessj. 

Similar to analyses, I often see computational results that I can tell within 5 seconds must be wrong. 

It is notoriously easy to be lulled into a sense of confidence in computed results, because 
the computer always gives you some result, and that result always looks good when plotted 
in a 3D shaded color orthographic projection. The corresponding “smoke test,” “function test,” 
and “performance test,” are as follows: 

1. Smoke test. Start the computer program running, and see if it crashes or not. If it doesn’t crash, 
you’ve passed the smoke test, part (a). Part (b) of the smoke test is to determine if the computed 
result passes the global conservation test. The goal of any program is to satisfy mass, momentum, 
energy and atom conservation at every point in the computational domain subject to certain constitutive 
relations (eg., Newton’s law of viscosity T x = \sdujdy), Hooke’s Law a = Es) and equations of state 
(eg, the ideal gas law.) This is a hard problem, and it is even hard to verify that the solution is 
correct once it is obtained. But it is easy to determine whether or not global conservation is satisfied, 
that is, 


• Is mass conserved, that is, does the sum of all the mass fluxes at the inlets, minus the mass 
fluxes at the outlets, equal to the rate of change of mass of the system (=0 for steady 
problems)? 

• Is momentum conserved in each coordinate direction? 

• Is energy conserved? 

• Is each type of atom conserved? 

If not, you are 100% certain that your calculation is wrong. You would be amazed at how many 
results are never “sanity checked” in this way, and in fact fail the sanity check when, after months or 
years of effort and somehow the results never look right, someone finally gets around to checking 
these things, the calculations fail the test and you realize all that time and effort was wasted. 

2. Performance test. Comes before the function test in this case. For computational studies, a critical 
performance test is to compare your result to a known analytical result under simplified conditions. For 
example, if you’re computing flow in a pipe at high Reynolds numbers (where the flow is turbulent), 
with chemical reaction, temperature-dependent transport properties, variable density, etc., first 
check your result against the textbook solution that assumes constant density, constant 
transport properties, etc., by making all of the simplifying assumptions (in your model) that the 
analytical solution employs. If you don’t do this, you really have no way of knowing if your 
model is valid or not. You can also use previous computations by yourself or others for testing, 
but of course there is no absolute guarantee that those computations were correct. 

3. function test. Similar to function test for analyses. 

By the way, even if you’re just doing a quick calculation, I recommend not using a calculator. Enter 
the data into an Excel spreadsheet so that you can add/change/scrutinize/save calculations as 
needed. Sometimes I see an obviously invalid result and when I ask, “How did you get that result? 
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What numbers did you use?” the answer is “I put the numbers into the calculator and this was the 
result I got.” But how do you know you entered the numbers and formulas correctly? What if you 
need to re-do the calculation for a slightly different set of numbers? 

Examples of the use of units and scrutiny 

These examples, particularly the first one, also introduce the concept of “back of the envelope” 
estimates, a powerful engineering tool. 

Example 1. Drag force and power requirements for an automobile 

A car with good aerodynamics has a drag coefficient (C D ) of 0.2. The drag coefficient is defined as 
the ratio of the drag force (F D ) to the dynamic pressure of the flow = V 2 pv (where p is the fluid 
density and v the fluid velocity far from the object) multiplied by the cross-section area (A) of the 
object, i.e. 

F D ~\c D pv-A 

The density of air at standard conditions is 1.18 kg/m 3 . 

(a) Estimate the power required to overcome the aerodynamic drag of such a car at 60 miles per 
hour. 

Power = Force x velocity 

v = 60 miles/hour x (5280 ft/mile) x (m/3.28 ft) x (hour/60 min) x (min/60 sec) = 26.8 m/s 
Estimate cross-section area of car as2mx3m = 6m 

F d = 0.5 x 0.2 x 1.18 kg/m' x (26.8 m/s) 2 x 6 m 2 = 510 kg m/s 2 = 510 Newton 

Power = F D x v = 510 kg m/s 2 x 26.8 m/s = 1.37 x 10 4 kg m 2 /s' = 1.37 x 10 4 W = 18.3 
horsepower, which is reasonable 

(b) Estimate the gas mileage of such a car. The heating value of gasoline is 4.4 x 10 J/kg and its 
density is 750 kg/m'. 

Fuel mass flow required = power (Joules/sec) / heating value (Joules/kg) 

= 1.37 x 10 4 kg m 2 /s' / 4.4 x 10 J/kg = 3.10 x 10 4 kg/s 

Fuel volume flow required = mass flow / density 

= 3.10 x 10’ 4 kg/s / 750 kg/m' = 4.14 x 10~ m'/s x (3.28 ft/m) ' x 7.48 gal/ft' 

= 1.09 x 10’ 4 gal/sec 

Gas mileage = speed / fuel volume flow rate = 

[(60 miles/hour)/(1.09 x 10~ 4 gal/sec)] x (hour / 3600 sec) = 152.564627113 miles/gallon 


(Equation 9) 
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Why is this value of miles/gallon so high? 

o The main problem is that conversion of fuel energy to engine output shaft work is about 
25% efficient at highway cruise conditions, thus the gas mileage would be 152.564627113 x 
0.25 = 38.1411567782 mpg 

o Also, besides air drag, there are other losses in the transmission, driveline, tires — at best the 
drivetrain is 80% efficient — so now we’re down to 30.51292542 mpg 
o Also — other loads on engine — air conditioning, generator, ... 

What else is wrong? There are too many significant figures; at most 2 or 3 are acceptable. When we 
state 30.51292542 mpg, that means we think that the miles per gallon is closer to 30.51292542 mpg 
than 30.51292541 mpg or 30.51292543 mpg. Of course we can’t measure the miles per gallon to 
anywhere near this level of accuracy. 31 is probably ok, 30.5 is questionable and 30.51 is ridiculous. 
You will want to carry a few extra digits of precision through the calculations to avoid round-off 
errors, but then at the end, round off your calculation to a reasonable number of significant figures 
based on the uncertainty of the most uncertain parameter That is, if I know the drag coefficient only to the 
first digit, i.e. I know that it’s closer to 0.2 than 0.1 or 0.3, but not more precisely than that, there is 
no point in reporting the result to 3 significant figures. 

Example 2. Scrutiny of a new formula 

I calculated for the first time ever the rate of heat transfer (q) (in watts) as a function of time t from 
an aluminum bar of radius r, length L, thermal conductivity k (units Watts/m°C), thermal diffusivity 

a (units m 2 /s), heat transfer coefficient h (units Watts/m 2o C) and initial temperature T bar conducting 
and radiating to surroundings at temperature T^ as 

<7 = k ( T bar ~ T Je a " r - hr 2 (T bar - T. -1) (Equation 10) 


Using “engineering scrutiny,” what “obvious” mistakes can you find with this formula? What is the 
likely “correct” formula? 


1. The units are wrong in the first term (Watts/m, not Watts) 


2. The units are wrong in the second term inside the parenthesis (can’t add 1 and something 
with units of temperature) 


3. The first term on the right side of the equation goes to infinity as the time (t) goes to 
infinity — probably there should be a negative sign in the exponent so that the whole term 
goes to zero as time goes to infinity. 


4. The length of the bar (L) doesn’t appear anywhere 


5. The signs on (T bar — T^) are different in the two terms — but heat must ALWAYS be 
transferred from hot to cold, never the reverse, so the two terms cannot have different signs. 
One can, with equal validity, define heat transfer as being positive either to or from the bar, 
but with either definition, you can’t have heat transfer being positive in one term and 
negative in another term. 
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2 

6. Only the first term on the right side of the equation is multiplied by the e ( ~ alh ) factor, 
and thus will go to zero as t —> oo. So the other term would still be non-zero even when t —> 

which doesn’t make sense since the amount of heat transfer (q) has to go to zero as t —> 

2 

°o. So probably both terms should be multiplied by the e ( ~ alh } factor. 

Based on these considerations, the probable correct formula, which would pass all of the 
smoke and function tests is 


‘l-[U-(T k ,-T.) + hr 2 (T lmr -Tj]e 


-at/C 


Example 3. Thermoelectric generator 

The thermal efficiency (r|) = (electrical power out) / (thermal power in) of a thermoelectric power 
generation device (used in outer planetary spacecraft (Figure 2), powered by heat generated from 
radioisotope decay, typically plutonium-238) is given by 


71 = 


( T \ 


7i 


H / 


’ 1 a ~ 


•J 1 + ZT .. +TJT, 


H 


2 


(Equation 11) 


where T is the temperature, the subscripts L, H and a refer to cold-side (low temperature), hot-side 
(high temperature) and average respectively, and Z is the “thermoelectric figure of merit”: 

o 2 


z = 


pk 


(Equation 12) 


where S is the Seebeck coefficient of material (units Volts/K, indicates how many volts are produced for 

each degree of temperature change across the material), p is the electrical resistivity (units ohm m) 
(not to be confused with density!) and k is the material’s thermal conductivity (W/mK). 

(a) show that the units are valid (passes smoke test) 


Everything is obviously dimensionless except for ZT a , which must itself be dimensionless so that 
I can add it to 1. Note 


S’ 2 

z = —— r = 

pk a 



^ Volt 
\ K 

\2 

[ohm m) 

w\ 
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( J / cout 

2 
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1 


Js 


coul 
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niK 


K = 


_ coul 2 1 K 


K 


J 2 1 5(1/5) 1 


= 1 OK 


coul 


K 


(b) show that the equation makes physical sense (passes function test) 

o If the material Z = 0, it produces no electrical power thus the efficiency should be zero. If Z 
= 0 then 


?? = 


( T \ 

1-^ 


T, 


H ) 


Vi+o?:, -1 

■fi + 0:r„ +TJT„ 


1- 



VT-i 

( T \ 

1 1l 

0 

T 

1 H / 

Vi +t,/t h ~ 

1 

k T 

\ / 

1 +T l /T h 


= 0 


OK 
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o If T L — T h , then there is no temperature difference across the thermoelectric material, and 
thus no power can be generated. In this case 



OK 


o Even the best possible material (ZT a °o) cannot produce an efficiency greater than the 
theoretically best possible efficiency (called the Carnot cycle efficiency, see page 85) = 1 — 
T l /T h , for the same temperature range. As ZT a —> oo, 



OK 


Side note #1: a good thermoelectric material such as Bi 2 Te 3 has ZT a ~ 1 and works up to about 
200 °C before it starts to melt, thus 


?? = 


( T \ 

1 -^ 


T 


H ) 


vr+1-1 


= 0.203 1- 


VT+I + (25 + 273) /(200 + 273) 
25 + 273 


= 0.203 


I T \ 

1-^ 


T 


H / 


= 0.203)7 


Carnot 


200 + 273 


= 0.0750 = 7.50% 


By comparison, your car engine has an efficiency of about 25%. So practical thermoelectric 
materials are, in general, not very good sources of electrical power, but are extremely useful in some 
niche applications, particularly when either (1) it is essential to have a device with no moving parts 
or (2) a “free” source of thermal energy at relatively low temperature is available, e.g. the exhaust of 
an internal combustion engine. 

Side note #2: a good thermoelectric material has a high S, so produces a large voltage for a small 
temperature change, a low p so that the resistance of the material to the flow of electric current is 
low, and a low k so that the temperature across the material AT is high. The heat transfer rate (in 
Watts) q = kAAT/Ax (see Chapter 8) where A is the cross-sectional area of the material and Ax is its 
thickness. So for a given AT, a smaller k means less q is transferred across the material. One might 
think that less q is worse, but no. Consider this: 

The electrical power = IV = (V/R)V = V 2 /R = (SAT) 2 /(pAx/A) = S 2 AT 2 A/pAx. 

The thermal power = kAAT/Ax 

The ratio of electrical to thermal power is [S 2 AT 2 A/pAx]/ [kAAT/Ax] = (S 2 /pk)AT = ZAT, 
which is why Z is the “figure of merit” for thermoelectric generators.) 
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Figure 2. Radioisotope thermoelectric generator used for deep space missions. Note that 
the plutonium-238 radioisotope is called simply, “General Purpose Heat Source.” 


Example 4. Density of matter 


Estimate the density of a neutron. Does the result make sense? The density of a white dwarf star is 

12 / 3 * * 

about 2x10" kg/m ' — is this reasonable? 


The mass of a neutron is about one atomic mass unit (AMU), where a carbon-12 atom has a mass of 
12 AMU and a mole of carbon-12 atoms has a mass of 12 grams. Thus one neutron has a mass of 

l -< -> \ 


(1 AMU) 


1 C -12 atom 


1 mole C-12 


12 g C -12 


12 AMU A6.02xl0 23 atoms C-12 A mole C-12 


1 kg 

1000 g 


= 1.66 x l(T 27 kg 


A neutron has a radius (r) of about 0.8 femtometer — 0.8 x 10 15 meter. Treating the neutron as a 
sphere, the volume is 47 tr 3 / 3 , and the density (p) is the mass divided by the volume, thus 


mass 

volume 


1.66x10 27 kg _ 775xl0 nM 
—t—(o.8 x KV^m)"' m 


* * 3/3 * ** 

By comparison, water has a density of 10 1 kg/m 1 , so the density of a neutron is far higher (by a 
factor of 10 4 ) than that of atoms including their electrons. This is expected since the nucleus of an 
atom occupies only a small portion of the total space occupied by an atom — most of the atom is 
empty space where the electrons reside. Also, even the density of the white dwarf star is far less 
than the neutrons (by a factor of 10'), which shows that the electron structure is squashed by the 
mass of the star, but not nearly down to the neutron scale (protons have a mass and size similar to 
neutrons, so the same point applies to protons too.) 
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Chapter 4. Statistics 


“There are three kinds of lies: lies, damn lies, and statistics... ” 
- Origin unknown, popularized by Mark Twain. 

Mean and standard deviation 


When confronted with multiple measurements y l5 y 2 , y 3 , ... of the same experiment (e.g. students’ 
scores on an exam), one typically reports at least two properties of the ensemble of scores, namely the 
mean value and the standard deviation'. 


Average or mean value — (sum of values of all samples) / number of samples 


^1 + V2 + V3 + - + y n 1 y 

7/1 7/7 


y - 


y< 


(Equation 13) 


i=l 


Standard deviation — square root of sum of squares of difference between each sample and the 
mean value, also called root-mean-square deviation, often denoted by the Greek letter lower case O: 



(>i - y ) 2 + (y 2 - y Y + (y 3 - y Y + • ■• ■•+ (y„ - y ) 

n -1 


rr\2 


—\2 


rr\2 



(Equation 14) 


Warning: in some cases a factor of n, not (n-1), is used in the denominator of the definition of 
standard deviation. I actually prefer n, since it passes the function test better: 

o With n in the denominator, then when n = 1, y 1 = y , and 0 = 0 (that is, no sample deviates 
at all from the mean value.) 

o With n - 1 in the denominator, then when n = 1, again y 1 = y , but now O = 0/0 and thus 
standard deviation is undefined 

But the definition using n — 1 connects better with other forms of statistical analysis that we won’t 
discuss here, so it is by far the more common definition. 

Example: 

On one of Prof. Ronney’s exams, the students’ scores were 50, 33, 67 and 90. What is the 
mean and standard deviation of this data set? 


A/r — 50 + 33 + 67 + 90 / . . . . . r 

Mean — -= oU (a bit tower than the average 1 prefer) 


(50 - 60) 2 + (33 - 60)' + (67 - 60) 2 + (90 - 60) 

4--1 


Standard deviation — 


.2 


= 31.12 
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Note also that (standard deviation)/mean is 31.12/60 = 0.519, which is a large spread. More 
typically this number for my exams is 0.3 or so. In a recent class of mine, the grade 
distribution was as follows: 


Grade 

# of standard deviations above/below mean 

A+ 

> 1.17 a above mean (1.90, 1.81) 

A 

0.84 to 1.17 above mean 

A- 

0.60 to 0.67 above mean 

B+ 

0.60 above mean to 0.10 below mean 

B 

0.32 to 0.29 below mean 

B- 

0.85 to 0.68 below mean 

C+ 

1.20 to 1.07 below mean 

c 

1.67 to 1.63 below mean 

C- 

> 1.67 below mean (2.04) 


Stability of statistics 

If I want to know the mean or standard deviation of a property, how many samples do we need? 
For example, if I flip a coin only once, can I decide if the coin is “fair” or not, that is, does it come 
up heads 50% of the time? Obviously not. So obviously I need more than 1 sample. Is 2 enough, 1 
time to come up heads, and another tails? Obviously not, since the coin might wind up heads or 
tails 2 times in a row. Below are the plots of two realizations of the coin-flipping experiment, done 
electronically using Excel. If you have the Word version of this file, you can double-click the plot to 
see the spreadsheet itself (assuming you have Excel on your computer.) Note that the first time 
time the first coin toss wound up tails, so the plot started with 0% heads and the second time the 
first coin was heads, so the plot started with 100% heads. Eventually the data smooths out to about 
50% heads, but the approach is slow. For a truly random process, one can show that the uncertainty 
decreases as 1/Vn, where n is the number of samples. So to have half as much uncertainty as 10 
samples, you need 40 samples! 




Figure 3. Results of two coin-toss experiments. 
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Side note: if a “fair” coin lands heads 100 times in a row, what are the chances of it landing heads 
on the 101 st flip? 50% of course, since each flip of a fair coin is independent of the previous one. 


Least-squares fit to a set of data 


Suppose you have some experimental data in the form of (x l5 y t ), (x x , y 2 ), (x 3 , y 3 ), ... (x n , yj and you 
think that the data should fit a linear relationship, i.e. y = mx + b, but in plotting the data you see 
that the data points do not quite fit a straight line. How do you decide what is the “best fit” of the 
experimental data to a single value of the slope m and y-intercept b? In practice this is usually done 
by finding the minimum of the sum of the squares of the deviation of each of the data points (x l5 y a ), 
(x x , y 2 ), (x 3 , y 3 ), ... (x n , yj from the points on the straight line (x l5 mxj+b), (x 2 , mx 2 +b), (x 3 , mx 3 +b), 

... ((x n , mx n +b). In other words, the goal is to find the values of m and b that minimize the sum 

S = (y r (m Xl +b)) 2 + (y 2 -(mx 2 +b)) 2 + (y 3 -(mx 3 +b)) 2 + ... + (y n -(mx n +b)) 2 . 

So we take the partial derivative of S with respect to m and b and set each equal to zero to find the 
minimum. Note: this is the ONLY place in the lecture notes where substantial use of 
calculus is made, so if you have trouble with this concept, don’t worry, you won’t use it 
again in this course. A partial derivative (which is denoted by a curly “5” compared to the straight 
“d” of a total derivative) is a derivative of a function of two or more variables, treating all but one of 
the variables as constants. For example if S(x, y, z) = x 2 y 3 — z 4 , then dS/dx = 2xy 3 , dS/dy = 3x 2 y 2 
and dS/dz = -4z 3 . So taking the partial derivatives of S with respect to m and b separately and 
setting both equal to zero we have: 
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(Equation 15) 


dS d 


(>’, - (mx l + b)f + (y 2 - ( mx 2 + b)f + (y 3 - (mx, + b)f +... + (y n - (mx n + b)) 


2 


db db 
d 

db 

=> (-2y, + 2 mx x + 2b) +... + (-2y fJ + 2 mx n + 2b) = 0 


= 0 


p 2 - 2 y x mx x - 2y x b + in x] + 2mx x b + b 1 j +... + fy 2 - 2y n mx n - 2y n b + m 2 x 2 + 2 mxjo + b" j 


n 


n n 


n 


n 


2m^x ( + 2^1-2^y ( . = 0 


m 


2^, +bn = ^y, 


i =1 


i=l i =1 


i=l 


i =1 


(Equation 16) 
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These are two simultaneous linear equations for the unknowns m and b. Note that all the sums are 
known since you know all the and y. These equations can be written in a simpler form: 

Cm + Ab = D 

Am + nb = B (Equation 17) 

/= 1 /= 1 /= l /= l 


These two linear equations can be solved in the usual way to find m and b: 


1 

m = — 
A 


/ 


B - n 


\ 


AD-BC 
A 2 - nC 


;b = 


AD-BC 
A 2 - nC 


(Equation 18) 


Example 


What is the best linear fit to the relationship between the height (x) of the group of students shown 
below and their final exam scores (y)? Assuming this trend was valid outside the range of these 
students, how tall or short would a student have to be to obtain a test score of 100? At what height 
would the student’s test score be zero? What test score would an amoeba (height ~ 0) obtain? 


Student name 

Height (x) (inches) 

Test score (y) (out of 100) 

Juanita Hernandez 

68 

80 

Julie Jones 

70 

77 

Ashish Kumar 

74 

56 

Fei Wong 

78 

47 

Sitting Bear 

63 

91 


A = 68+70+74+78+63 = 353 
B = 80 + 77 + 56 + 47 + 91 = 351 
C = 68 2 + 70 2 + 74 2 + 78 2 + 63 2 = 25053 
D = 68*80 + 70*77 + 74*56 + 78*47 + 63*91 = 24373 


From which we can calculate m = -3.107, b = 289.5, i.e. 

Test score — -3.107 ^Height +289.5 

For a score of 100, 100 = -3.107*Height + 289.5 or Height = 61.01 inches = 5 feet 1.01 inches 
For a score of zero, 0 = -3.107*Height + 289.5 or Height = 93.20 inches = 7 feet 9.2 inches 
For a height of 0, score = 289.5 


21 





Figure 4. Least-squares fit to data on test score vs. height for a hypothetical class 

How does one determine how well or poorly the least-square fit actually fits the data? That is, how 
closely are the data points to the best-fit line? The standard measure is the so-called R "-value defined 
as one minus the sum of the squares of the deviations from the fit just determined (i.e. the sum of 
(^-(mXj+b))" divided by the sum of the squares of the difference between y { and the average value y 
(=70.2 for this case), i.e., 



(Equation 19) 


• • • ^ • 

For a perfect fit y { = mxj + b for all i, so the sum in the numerator is zero, thus R" = 1 is a perfect fit. 

The example shown above is pretty good, 


R 2 = 1 - 


= 0.9631 


80 -(- 3 . 107*68 + 298 . 5)) 2 + ( 77 -(- 3 . 107*70 + 298 . 5)) 2 + ( 56 -(- 3 . 107*74 + 298 . 5)) 2 + ( 47 -(- 3 . 107*78 + 298 . 5)) 2 + ( 91 -(- 3 . 107*63 + 298 . 5 )) j 

(80 - 70 . 2) 2 + (77 - 70 . 2) 2 + (56 - 70 . 2) 2 + (47 - 70 . 2) 2 + (91 - 70 . 2) 2 


and even fairly crummy fits (i.e. as seen visually on a plot, with many of the data points far removed 
from the line) can have R > 0.9. So R" has to be pretty close to 1 before it’s really a good-looking 
fit. 
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Chapter 5. Forces in structures 

“The Force can have a strong influence on the weak-minded ” 

- Ben Obi-wan Kenobi, explaining to Luke Skywalker how he made the famous “these aren’t the 
Droids you’re looking for” trick work. 

Main course in AME curriculum on this topic: AME 201 (Statics). 

Forces 

Forces acting on objects are vectors that are characterized by not only a magnitude (e.g. pounds force 
or Newtons) but also a direction. A force vector F (vectors are usually noted by a boldface letter) 
can be broken down into its components in the x, y and z directions in whatever coordinate system 
you’ve drawn: 

F = F x i + F y j + F z k Equation 20 

Where F x , F and F z are the magnitudes of the forces in the x, y and z directions and i, j and k are 
the unit vectors in the x, y and z directions (i.e. vectors whose directions are aligned with the x, y and z 
coordinates and whose magnitudes are exactly 1 (no units)). 

Forces can also be expressed in terms of the magnitude = (F x + F + F z ) and direction relative 
to the positive x-axis (= tan 4 (F y /F x ) in a 2-dimensional system). Note that the tan 4 (F y /F x ) function 
gives you an angle between +90° and -90° whereas sometimes the resulting force is between +90° 
and +180° or between -90° and -180°; in these cases you’ll have to examine the resulting force and 
add or subtract 180° from the force to get the right direction. 


Moments of forces 


Some types of structures can only exert forces along the line connecting the two ends of the 


A 


\d 




♦ 

Line of action 



Figure 5. Force, line of action and moment (= Fd) about a point A 
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structure, but cannot exert any force perpendicular to that line. These types of structures include 
ropes, ends with pins, and bearings. Other structural elements can also exert a force perpendicular 
to the line (Figure 5). This is called the moment of a force, which is the same thing as torque. Usually 
the term torque is reserved for the forces on rotating, not stationary, shafts, but there is no real 
difference between a moment and a torque. 

The distinguishing feature of the moment of a force is that it depends not only on the vector force 
itself (Fj) but also the distance (dj from that line of force to a reference point A. (I like to call this 
distance the moment arm) from the anchor point at which it acts. If you want to loosen a stuck bolt, 
you want to apply whatever force your arm is capable of providing over the longest possible d,. The 
line through the force is called the line of action. The moment arm is the distance (d; again) 
between the line of action and a line parallel to the line of action that passes through the anchor 
point. Then the moment of force (MJ is defined as 

Mi = Fjdj Equation 21 

where is the magnitude of the vector F. Note that the units of M () is force x length, e.g. ft lbf or 
N m. This is the same as the unit of energy, but the two have nothing in common — it’s just 
coincidence. So one could report a moment of force in units of Joules, but this is unacceptable 
practice — use N m, not J. 

Note that it is necessary to assign a sign to M ; . Typically we will define a clockwise moment as 
positive and counterclockwise as negative, but one is free to choose the opposite definition — as long 
as you’re consistent within an analysis. 

In order to have equilibrium of an object, the sum of all the forces AND the moments of the forces 
must be zero. In other words, there are two ways that a 2-dimensional object can translate (in the x 
and y directions) and one way that in can rotate (with the axis of rotation perpendicular to the x-y 
plane.) So there are 3 equations that must be satisfied in order to have equilibrium, namely: 

n n n 

2) F v,i = °;2 = 0; 21 M ‘ = 0 Equation 22 

1=1 /= 1 i=\ 

Note that the moment of forces must be zero regardless of the choice of the origin (i.e. not just at the 
center of mass). So one can take the origin to be wherever it is convenient (e.g. make the moment 
of one of the forces = 0.) Consider the very simple set of forces below: 
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A 



141.4 Ibf 


200 Ibf 
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0.5 ft 


0.707 ft 


0.5 ft 


/*<r 

0.707 ft X 


.i.-- 

D 


141.4 Ibf 


Figure 6. Force diagram showing different ways of computing moments 

Because of the symmetry, it is easy to see that this set of forces constitutes an equilibrium condition. 
When taking moments about point ‘B’ we have: 


2F X = +141.4 cos(45°) Ibf + 0 - 141.4 cos(45°) Ibf = 0 

2F y = +141.4 sin(45°) Ibf - 200 Ibf + 141.4 sin(45°) = 0 

2M b = -141.4 Ibf* 0.707 ft - 200 Ibf* 0 ft +141.4 Ibf* 0.707 ft = 0. 

But how do we know to take moments about point B? We don’t. But notice that if we take 
moments about point ‘A’ then the force balances remain the same and 

2M a = -141.4 Ibf* Oft -200 Ibf* 1 ft+ 141.4 Ibf* 1.414 ft = 0. 

The same applies if we take moments about point ‘C\ or a point along the line ABC, or even a point 
NOT along the line ABC. For example, taking moments about point C D’, 

2M d = -141.4 Ibf* (0.707 ft + 0.707 ft) + 200 Ibf* 0.5 ft +141.4 Ibf* 0.707 ft = 0 

The location about which to take the moment can be chosen to make the problem as simple as 
possible, e.g. to make some of the moments of forces = 0. 

Example of “why didn't the hook just say that The state of equilibrium merely requires that 3 
constraint equations are required. There is nothing in particular that requires there be 2 force and 1 
moment constraint equations. So one could have 1 force and 2 moment constraint equations: 

M /=1 >1 

where the coordinate direction x can be chosen to be in any direction, and moments are taken about 
2 separate points A and B. Or one could even have 3 moment equations: 


Equation 23 
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Equation 24 


i = 1 7=1 7=1 

Also, there is no reason to restrict the x and y coordinates to the horizontal and vertical directions. 
They can be (for example) parallel and perpendicular to an inclined surface if that appears in the 
problem. In fact, the x and y axes don’t even have to be perpendicular to each other, as long as they 
are not parallel to each other, in which case 2F X = 0 and 2F V = 0 would not be independent 
equations. 

This is all fine and well for a two-dimensional (planar) situation, what about ID or 3D? For ID 
there is only one direction that the object can move linearly and no way in which it can rotate. For 
3D, there are three directions it can move linearly and three axes about which it can rotate. Table 1 
summarizes these situations. 


# of spatial 
dimensions 

Maximum # of 
force balances 

Minimum # of 
moment balances 

Total # of unknown 
forces & moments 

i 

i 

0 

i 

2 

2 

1 

3 

3 

3 

3 

6 


Table 1. Number of force and moment balance equations required for static equilibrium as 
a function of the dimensionality of the system. (But note that, as just described, moment 
equations can be substituted for force balances.) 


Types of forces and moments 

A free body diagram is a diagram showing all the forces and moments of forces acting on an object. 

We distinguish between two types of objects: 

1. P'articles that have no spatial extent and thus have no moment arm (d). An example of this 
would be a satellite orbiting the earth because the spatial extent of the satellite is very small 
compared to the distance from the earth to the satellite or the radius of the earth. Particles 
do not have moments of forces and thus do not rotate in response to a force. 

2. Rigid bodies that have a finite dimension and thus has a moment arm (d) associated with each 
applied force. Rigid bodies have moments of forces and thus can rotate in response to a 
force. 

There are several types of forces that act on particles or rigid bodies: 

1. Rope , cable y etc . — Force (tension) must be along line of action; no moment (1 unknown 
force) 

T< -• 
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2. Rollers, frictionless surface — Force must be perpendicular to the surface; no moment (1 
unknown force). There cannot be a force parallel to the surface because the roller would 
start rolling! Also the force must be away from the surface towards the roller (in other 
words the roller must exert a force on the surface), otherwise the roller would lift off of the 
surface. 



3. Frictionless pin or hinge — Force has components both parallel and perpendicular to the 
line of action; no moment (2 unknown forces) (note that the coordinate system does not 
need to be parallel and perpendicular to either the gravity vector or the bar) 



4. Fixed support — Force has components both parallel and perpendicular to line of action 
plus a moment of force (2 unknown forces, 1 unknown moment). Note that for our simple 
statics problems with 3 degrees of freedom, if there is one fixed support then we already 
have 3 unknown quantities and the rest of our free body cannot have any unknown forces if 
we are to employ statics alone to determine the forces. In other words, if the free body has 
any additional unknown forces the system is statically indeterminate as will be discussed shortly. 



5. Contact friction — Force has components both parallel (F) and perpendicular (N) to surface, 
which are related by F = pN, where p is the coefficient of friction , which is usually assigned 
separate values for static (no sliding) (p s ) and dynamic (sliding) (p d ) friction, with the latter 
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being lower. (2 unknown forces coupled by the relation F = pN). p depends on both of the 
surfaces in contact. Most dry materials have friction coefficients between 0.3 and 0.6 but 
Teflon in contact with Teflon, for example, can have a coefficient as low as 0.04. Rubber 
(e.g. tires) in contact with other surfaces (e.g. asphalt) can yield friction coefficients of almost 



Actually the statement F = p s N for static friction is not correct at all, although that’s how it’s almost 
always written. Consider the figure on the right, above. If there is no applied force in the horizontal 
direction, there is no need for friction to counter that force and keep the block from sliding, so F = 
0. (If F 7^ 0, then the object would start moving even though there is no applied force!) Of course, 
if a force were applied (e.g. from right to left, in the —x direction) then the friction force at the 
interface between the block and the surface would counter the applied force with a force in the +x 

direction so that 2F X = 0. On the other hand, if a force were applied from left to right, in the +x 

direction) then the friction force at the interface between the block and the surface would counter 

the applied force with a force in the -x direction so that ZF X = 0. The expression F = p s N only 

applies to the maximum magnitude of the static friction force. In other words, a proper 

statement quantifying the friction force would be | F | < p s N, not F = p s N. If any larger force is 

applied then the block would start moving and then the dynamic friction force F = p d N is the 

applicable one — but even then this force must always be in the direction opposite the motion — so 

| F | = p d N is an appropriate statement. Another, more precise way of writing this would be 

- v _. . . . v 

F = fjL d N —, where v is the velocity of the block and vis the magnitude of this velocity, thus —is a 

v v 

unit vector in the direction of motion. 

Special note: while ropes, rollers and pins do not exert a moment at the point of contact, you can 
still sum up the moments acting on the free body at that point of contact. In other words, ZM A = 
0 can be used even if point A is a contact point with a rope, roller or pin joint, and all of the other 
moments about point A (magnitude of force x distance from A to the line of action of that force) 
are still non-zero. Keep in mind that A can be any point, within or outside of the free body. It does 
not need to be a point where a force is applied, although it is often convenient to use one of those 
points as shown in the examples below. 

Statically indeterminate system 

Of course, there is no guarantee that the number of force and moment balance equations will be 
equal to the number of unknowns. For example, in a 2D problem, a beam supported by one pinned 
end and one roller end has 3 unknown forces and 3 equations of static equilibrium. However, if 
both ends are pinned, there are 4 unknown forces but still only 3 equations of static equilibrium. 
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Such a system is called statically indeterminate and requires additional information beyond the equations 
of statics (e.g. material stresses and strains, discussed in the next chapter) to determine the forces. 


Analysis of statics problems 

A useful methodology for analyzing statics problems is as follows: 

1. Draw a free body diagram — a free body must be a rigid object, i.e. one that cannot bend in 
response to applied forces 

2. Draw all of the forces acting on the free body. Is the number of unknown forces equal to 
the total number of independent constraint equations shown in Table 1 (far right column)? 
If not, statics can’t help you. 

3. Decide on a coordinate system. If the primary direction of forces is parallel and 
perpendicular to an inclined plane, usually it’s most convenient to have the x and y 
coordinates parallel and perpendicular to the plane, as in the cart and sliding block examples 
below. 

4. Decide on a set of constraint equations. As mentioned above, this can be any combination 
of force and moment balances that add up to the number of degrees of freedom of the 
system (Table 1). 

5. Decide on the locations about which to perform moment constraint equations. Generally 
you should make this where the lines of action of two or more forces intersect because this 
will minimize the number of unknowns in your resulting equation. 

6. Write down the force and moment constraint equations. If you’ve made good choices in 
steps 2 — 5, the resulting equations will be “easy” to solve. 

7. Solve these “easy” equations. 


Example 1. Ropes 

Two tugboats, the Monitor and the Merrimac, are pulling a Peace Barge due west up Chesapeake 
Bay toward Washington DC. The Monitor’s tow rope is at an angle of 53 degrees north of due west 
with a tension of 4000 lbf. The Merrimac’s tow rope is at an angle of 34 degrees south of due west 
but their scale attached to the rope is broken so the tension is unknown to the crew. 
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Figure 7. Free body diagram of Monitor-Merrimac system 

a) What is the tension in the Merrimac’s tow rope? 

Define x as positive in the easterly direction, y as positive in the northerly direction. In order for 
the Barge to travel due west, the northerly pull by the Monitor and the Southerly pull by the 
Merrimac have to be equal, or in other words the resultant force in the y direction, R^, must be 
zero. The northerly pull by the Monitor is 4000 sin(53°) = 3195 lbf. In order for this to equal 
the southerly pull of the Merrimac, we require F Merrimac sin(34°) = 3195 lbf, thus F Merrimac = 5713 
lbf. 

b) What is the tension trying to break the Peace Barge (i.e. in the north-south direction)? 

This is just the north/south force just computed, 3195 lbf 

c) What is the force pulling the Peace Barge up Chesapeake Bay? 

The force exerted by the Monitor is 4000 cos(53°) = 2407 lbf. The force exerted by the 
Merrimac is 5713 cos(34°) = 4736 lbf. The resultant is R x = 7143 lbf. 

d) Express the force on the Merrimac in polar coordinates (resultant force and direction, with 0° 
being due east, as is customary) 

The magnitude of the force is 5713 lbf as just computed. The angle is -180° - 34° = -146°. 

Example 2. Rollers 

A car of weight W is being held by a cable with tension T on a ramp of angle 0 with respect to 
horizontal. The wheels are free to rotate, so there is no force exerted by the wheels in the direction 
parallel to the ramp surface. The center of gravity of the vehicle is a distance “c” above the ramp, a 
distance “a” behind the front wheels, and a distance “b” in front of the rear wheels. The cable is 
attached to the car a distance “d” above the ramp surface and is parallel to the ramp. 
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Figure 8. Free body diagram for car-on-ramp with cable example 


(a) What is the tension in the cable in terms of known quantities, i.e. the weight W, dimensions a, b, 
c, and d, and ramp angle 0? 

Define x as the direction parallel to the ramp surface and y perpendicular to the surface as 
shown. The forces in the x direction acting on the car are the cable tension T and 
component of the vehicle weight in the x direction = WsinB, thus 2F X = 0 yields 

Wsin0 - T = 0 => T = Wsin0 

(b) What are the forces where the wheels contact the ramp (F A and F B )? 

The forces in the y direction acting on the car are F y A , F B and component of the vehicle 
weight in the y direction = Wcos0. Taking moments about point A, that is 2M A = 0 (so that 
the moment equation does not contain F A which makes the algebra simpler), and defining 
moments as positive clockwise yields 

(Wsin0)(c) + (Wcos0)(a) - F y>B (a+b) - T(d) = 0 

Since we already know from part (a) that T = Wsin0, substitution yields 

(Wsin0)(c) + (Wcos0)(a) - F yB (a+b) - (Wsin0)(d) = 0 

Since this equation contains only one unknown force, namely F B , it can be solved directly to 
obtain 



acos(0) + (c- d)sin(6) 
a + b 
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Finally taking 2F y — 0 yields 


F y ,A + F y , B - WCOS0 = 0 

Which we can substitute into the previous equation to find F A 



W bcos(0) - (c - d)sin(6) 

a + b 


Note the function tests: 

1) For 0 = 0, T = 0 (no tension required to keep the car from rolling on a level road) 

2) As 0 increases, the tension T required to keep the car from rolling increases 

3) For 0 = 90°, T = W (all of the vehicle weight is on the cable) but note that F A and 
F B are non-zero (equal magnitudes, opposite signs) unless c = d, that is, the line of 
action of the cable tension goes through the car’s center of gravity. 

4) For 0 = 0, F A = (b/(a+b)) and F B = (a/(a+b)) (more weight on the wheels closer 
to the center of gravity.) 

5) Because of the — sign on the 2 nd term in the numerator of F A (-(c-d)sin(0)) and the + 
sign in the 2 nd term in the numerator of F B (+(c-d)sin(0)), as 0 increases, there is a 
transfer of weight from the front wheels to the rear wheels. 

Note also that F A < 0 for b/(c-d) < tan(0), at which point the front (upper) wheels lift off 
the ground, and that F y A < 0 for a/(d-c) > tan(0), at which point the back (lower) wheels lift 
off the ground. In either case, the analysis is invalid. (Be aware that c could be larger or 
smaller than d, so c-d could be a positive or negative quantity.) 


Example 3. Friction 

A 100 lbf acts on a 300 lbf block placed on an inclined plane with a 3:4 slope. The coefficients of 
friction between the block and the plane are p s = 0.25 and ju d = 0.20. 

a) Determine whether the block is in equilibrium 

b) If the block is not in equilibrium (i.e. it’s sliding), find the net force on the block 

c) If the block is not in equilibrium, find the acceleration of the block 
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Figure 9. Free body diagram for sliding-block example 

(a) To maintain equilibrium, we require that ZF X = 0 and 2F = 0. Choosing the x direction parallel 
to the surface and y perpendicular to it, 

2F y = N - (4/5)(300 lbf) = 0 => N = 240 lbf 
so the maximum possible friction force is 

Ffrktion, max = = 0.25 * 240 lbf = 60 lbf. 

The force needed to prevent the block from sliding is 

2F X = 100 lbf- (3/5)(300 lbf) + F needed = 0 
Fneeded = ^00 lbf + (3/5)(300 lbf) = 80 lbf 

Which is more than the maximum available friction force, so the block will slide down the plane. 

(b) The sliding friction is given by 

Ffricdon, max = K,N = 0.20 * 240 lbf = 48 lbf 

so the net force acting on the block in the x direction (not zero since the block is not at equilbrium) 
is 

2F X = 100 lbf- (3/5)(300 lbf) + 48 lbf = -32 lbf 

(c) F = ma => -32 lbf = 300 lbm * acceleration 

acceleration = (-32 lbf/300 lbm) ??? 

what does this mean? lbf/lbm has units of force/mass, so it is an acceleration. But how to convert 
to something useful like ft/sec 2 ? Multiply by 1 in the funny form of g c = 1 = 

32.174 lbm ft / lbf sec 2 , of course! 
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acceleration = (-32 lbf/300 lbm) (32.174 lbm ft / lbf sec 2 ) = -3.43 ft/sec 2 
or, since g earth = 32.174 ft/sec 2 , 

acceleration = (-3.43 ft/sec 2 )/(32.174 ft/sec 2 g earth ) = -0.107 g earth . 

The negative sign indicates the acceleration is in the —x direction, i.e. down the slope of course. 

A good function test is that the acceleration has to be less than 1 g earth , which is what you would get 
if you dropped the block vertically in a frictionless environment. Obviously a block sliding down a 
slope (not vertical) with friction and with an external force acting up the slope must have a smaller 
acceleration. 

Example 4. Rollers and friction 

A car of weight W is equipped with rubber tires with coefficient of static friction p s . Unlike the 
earlier example, there is no cable but the wheels are locked and thus the tires exert a friction force 
parallel to and in the plane of the ramp surface. As with the previous example, the car is on a ramp 
of angle 0 with respect to horizontal. The center of gravity of the vehicle is a distance “c” above the 
ramp, a distance “a” behind the front wheels, and a distance “b” in front of the rear wheels. 



Figure 10. Free body diagram for car-on-ramp with friction 

example 


(a) What is the minimum p s required to keep the car from sliding down the ramp? 

The unknowns are the resulting forces at the wheels (F y A and F B ) and the coefficient of 
friction p s . Taking ZF X = 0, ZF = 0 and 2M A = 0 yields, respectively, 
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-^ s F y ,A - ^ s F y ,B + Wsin0 - 0 


F y>A + F y>B - Wcos0 = 0 


(Wsin0)(c) + (Wcos0)(a) - F yB (a+b) = 0 
Which may be solved to obtain 



- W fcc ° s(e >- Csin(e) ;F, J , - W acos(e) + f sin(e) ; M , - “5^ , tan(e) 


a + b 


a + b 


cos(0) 


Note the function tests 


1) For 0 — 0, |u s — tan(0) — 0 (no friction required to keep the car from sliding on a 
level road) 

2) As 0 increases, the friction coefficient p s required to keep the car from sliding 
increases 

3) For 0 = 0, F y A = (b/(a+b)) and F B = (a/(a+b)) (more weight on the wheels closer 
to the center of gravity 

4) Because of the — sign on the 2 nd term in the numerator of F A (-c sin(0)) and the + 

sign in the 2 nd term in the numerator of F B (+c sin(0)), as 0 increases, there is a 
transfer of weight from the front wheels to the rear wheels. 

Note also that we could have also tried 2F = 0, ZM A and 2M B = 0: 

F y ,a + F yjB - WcosO = 0 

(Wsin0)(c) + (WcosO)(a) - F >yB (a+b) = 0 

(Wsin0)(c) - (WcosO) (b) + F vA (a+b) = 0 

In which case, the second equation could have been subtracted from the third to obtain: 

F y , A + F y , B - WcosO = 0 

which is the same as the first equation. So the three equations are not independent of each 
other, and we can’t solve the system. What’s wrong? The coefficient of friction fi s doesn't appear in 

the set of equations 2!F X = 0, 2M A and HM B = 0. We need to have each of the three unknowns 
F y ,A> F B and p s in at least one of the three equations. The set ZF X — 0, 2M A and 2M B — 0 
doesn’t satisfy that criterion. 

(b) At what angle will the car tip over backwards, assuming that it doesn’t start sliding down the 
ramp at a smaller angle due to low p s ? 

This will occur when F yA = 0, i.e. when sin(0)/cos(0) = tan(0) = b/c. This is reasonable 
because the tip-over angle should increase when c is made larger (center of gravity closer to 
the ground) or b made smaller (center of gravity shifted forward). Notice also that it doesn’t 
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depend on p s , that is, as long as it doesn’t slide due to low ju s , the tip-over angle only depends 
on the force balance. 

For what it’s worth, also note that the tip-over angle equals the sliding angle when tan(0) = 
p s = b/c. Since generally p s << 1, Except for a very top-heavy (large c) or rear-weight- 
shifted (small b) vehicles, the vehicle will slide down the ramp before it flips over backwards. 

Example 5. Pinned joint 



Figure 11. Free body diagram for pinned joint example 


A straight bar of negligible mass 12 inches long is pinned at its lower end (call it point A) and has a 
roller attached to its upper end (call it point B) as shown in the figure. The bar is at a 30° angle from 
horizontal. A weight of 100 lbf is hung 4 inches from the lower end (call it point C). 

a) What are the forces in the x and y directions on the pinned end? What is the force in the x 
direction on the roller end? 


The pinned end can sustain forces in both the x and y directions, but no moment. The roller 
end can sustain a force only in the x direction, and again no moment. Summing the forces in the 
y direction 


F y ,A + F y , B + F yC = 0 F YiA + 0 - 100 lbf = 0 ^ F y , A = +100 lbf. 

In other words, in the y direction the vertical force at point A must be +100 lbf since that is the 
only force available to counteract the 100 lbf weight. Next, taking moments about point A 
(since the lines of action of two of the unknown forces intersect at point A), 

2M a = 0 => -(4 in)(cos(30°))(-100 lbf) + (6 m)F xB = 0 => F xB = -57.7 lbf. 

Finally, for force balance in the x direction, 





F XjA = -F x ,b - F x ,c = -(-57.7 lbf) - 0 = +57.7 lbf 
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b) Would the forces change if the roller and pinned ends were reversed? 


In this case summing the forces in the x direction: 

F x ,a + F x ,b + F xC = 0 => 0 + F x>b + 0 = 0 => F xB = 0. 
For force balance in the y direction, 


F y>A + F y . B + F vC = 0 => F v-A + F v . b = 100 lbf 


Y>' 


y> 


Taking moments about point C just for variety (not the easiest way, since neither F A nor F B are 
known, we just know that F y A + F B =100 lbf), 

2M C = 0 => (4 in)(cos(30°))F yA - (8 in)(cos(30°))F y , B + (8 in)(sin( 30 °))F x , B = 0 

=> (4 in)(cos(30°))(100 lbf- F yB ) - (8 in)(cos(30°))F yjB + 0 = 0 

=> F y>B = +33.3 lbf => F yjA = +66.7 lbf 

which is quite different from case (a). 

c) What would happen if the lower end were fixed rather than pinned (upper end having the roller 
again)? 


In this case there are 4 unknown quantities (F x A , F y A , M A and F x B ) but only 3 equations (ZF X — 0, 

EF V = 0, 2M = 0) so the system is statically indeterminate. If one takes away the roller end 
entirely, then obviously F yA = 100 lbf, F xA = 0 and M A = +(100 lbf)(4 in)(cos(30°)) = 346.4 in 
lbf. 
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Chapter 6. Stresses, strains and material properties 

“It is not stress that kills us, it is our reaction to it’’ 

- Hans Selye, endocrinologist. 

Main course in AME curriculum on this topic: AME 204 (Strength of Materials). 

Stresses and strains 

As a follow-on to the discussion of statics we need to consider if a material is subject to a given 
tensile, compressive or side load, how much stress is imparted into the material, and will this stress 
cause it to break? As with many subjects in the class, you will learn about this in a lot more detail in 
future classes; here you’ll just get a taste of it. 

The normal stress (a) in a material is defined as 


a = F/A Equation 25 

where F is the force (either tension or compression) acting perpendicular to an imaginary plane 
surface passing through a piece of material and A is the cross section area. It is called “normal” not 
in the sense of being “typical” or “standard” but in the sense of being perpendicular or orthogonal 
to the cross-section of the material. Stress is defined as positive if the material is in tension (i.e. the 
material is being pulled apart) and negative if the material is in compression. Stress has units 
force/area, i.e. the same as pressure. The units are typically N/m 2 or lbf/in 2 . Often the unit of 
“kips” (kilopounds per square inch = 1000 lbf/in - ) is used to report stress. 

The strain (e) is the fractional amount of elongation or contraction in a material caused by a stress. 
For example, if under a given amount of tensile stress, a steel bar stretches from a length (L) of 1.00 
inch to 1.01 inch (a change in length, AL, of 0.01 inch) the strain = (1.01 — 1.00)/1.00 = 0.01. In 
other words, 


8 = AL/L Equation 26 

For most materials (other than gooey ones, i.e. Silly Putty™, Play-Doh™, ...) the amount of strain 
before failure is relatively small (i.e. less that 0.1, meaning that the material deforms less than 10% 
before failing.) 

An elastic material has a linear relationship between stress and strain, i.e. 


a — bE Equation 27 

where E is called the elastic modulus , i.e. the slope of the plot of O vs. 8 in the elastic region shown in 
Figure 12. Note that since 8 is dimensionless, E also has units of pressure. 

The strength of a material is generally reported in terms of the stress it withstands. For a sufficiently 
small stress, materials return to their original length or shape after the stress is removed. The 
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smallest stress for which the material does not return to its original length or shape after the stress is 
removed is called the yield stress (a yield ). Beyond this stress, generally the slope of the O vs. 8 plot 
becomes smaller. There is often an increase in slope as 8 is increased still further, up to a maximum 
a called the “ultimate stress”, beyond which a actually decreases as 8 increases, leading finally to 
fracture of the material. 

Note that we can write Equation 30 in the form F/A = (AL/L)E or F = (EA/L)(AL), which looks 
just like the force on a linear spring, F = kx, with k = EA/L. (One might wonder what happened to 
the — sign, that is, isn’t F = -kx? Stress is usually defined as positive in tension where as for the 
spring F is defined as positive in compression.) So E and the material dimensions A and L 
determine its “spring constant.” Figure 12 is typical of a ductile material such as steel that deforms 
significantly before failure. This is by no means the only shape that a vs. 8 curves may have. A 

brittle material such as a ceramic or concrete will fail without significant yielding, that is, the a vs. 8 
curve is nearly linear up to the failure point. This doesn’t mean that ceramics are necessary weak, in 
fact they may have higher E than ductile materials, but they are unforgiving to over-stressing (really, 
over-straining.) 

E and yield or ultimate strength have the same units (Pa or lbf/in 2 ) but there is no particular 
relationship between E and strength. Materials can be hard (high E) but break easily (low strength) 
or vice versa. Some examples of material properties are shown in Table 2. This table shows 
materials that are more or less isotropic , i.e. their properties are similar no matter what direction stress 
is applied relative to the material. Many engineering materials are anisotropic , i.e. they are not 
isotropic. A typical example of such materials is graphite-epoxy composites composed of fibers of 
graphitic carbon (which have very high tensile strength in the plane of the graphite sheets, and low 
strength in the direction perpendicular to this plane) that are bonded to an epoxy polymer, which 
has relatively low tensile strength but high compressive and shear strength. The result is a material 
that has very high strength for its weight. (The Boeing 787 uses composites for most of the 
structure; this has the advantage of high strength to weight ratio, no possibility of corrosion, and 
ease of forming into any desired shape.) 




proportionality limit 
elastic limit 
yield stress 


ultimate stress 


yielding 


strain hardening 


necking 


plastic behavior 


elastic 
behavi or 


fracture 


G 

stress 

/Pa 


a u 

Ultimate stress 


Fracture stress 

Oy 

Yield stress 

Op| 

Proportionality limit 


strain 

e 


A brittle material. This material is also strong because there is little 
strain for a high stress. The fracture of a brittle material is sudden with 
little or no plastic deformation.. Glass is brittle 


A strong material which is not ductile. Steel 
wires stretch very little, and break suddenly. 



A ductile material 

after the elastic region 
there is a strange 
section where 'necking' 
occurs - permanent 
deformation occurs In 
this ‘plastic region" 


A plastic material 
very small elestic 
region. 


s (strain) 


Figure 12. Typical stress-strain relationships: left: ductile material; right: comparison of various 
types of materials. Sources: http:/ / dolbow.cee.duke.edu/TENSILE/tutorial/node4.html , 
http://www.cyberp hysics.co.uk/topics/forces/young_modulus.htm 
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Material 

E (10 9 Pa) 

V 

Yield strength (in tension 
unless otherwise noted) 

(10 6 Pa) 

Ultimate 
strength 
(10 6 Pa) 

Aluminum, 6061-T6 

68.9 

0.32 

276 

310 

Steel, 4340-HR 

200 

0.30 

910 

1041 

Iron, pure 

200 

0.29 

30 

540 

Diamond 

700 - 1200 

0.10-0.29 

8680 - 16530 (compressive) 


High-density 

polyethylene 

0.18-1.6 


2.4-31.7 

10-50 

Alumina, A1 2 0 3 

370 

0.22 

3000 (compressive) 

300 

Solder (60% tin, 40% 
lead) 

30 

0.4 


53 

Silica aerogel 

0.001-0.01 

0.2 


0.016 


Table 2. Properties of some common materials (from http://www.matweb.com) 

Of course, in any design one must employ a material with a yield strength greater than the actual 
stress in the material; the ratio of the yield stress of the material to the actual predicted stress in the 
material is call the factor of safety. 

Some factoids about materials 


How strong are these materials? How does this compare with the strength of the attractive forces 

between the atoms (F atoms )? That is, can we estimate the strength of the material a ~ F at 0 ms/- A atoms , 
where A atoms is the cross-section area of the atoms? How could we estimate F atoms and A atoms based 
on macroscopically measurable properties? Let’s start with the size of one atom. Let’s consider a 
typical material like aluminum. Its molecular weight is 27 g/mole, and its density is (by coincidence) 

/ 3 73 * * 

2.7 g/cm 1 , and 1 mole = 6.02 x 10“' atoms. So the volume occupied by each atom is 


cm 3 27 g mole 

2.7 g mole 6.02 x 10 23 atom 


= 1.66 x 10 -23 


cm 3 
atom 


^ ^ * * ^ 31/3 8 10 

Thus each atom occupies a roughly cubic space of (1.66 x 10’“’) = 2.55 x 10 1 cm = 2.55 x 10’ m, 

* 10 2 20 2 * * 

or a cross-section area of (2.55 x 10’ m)“ = 6.51 x 10 “ m . What is the attractive force between the 

atoms? The heat of formation of Al(gas) from Al(s) is 330 kj/mole (see http://webbook.nist.gov) . 
This is the energy needed to separate A1 atoms in the solid phase from each other to make a gas. 
On a per-atom basis this is (330,000 J/mole)(mole/6.02 x 10 2 ’ atoms) = 5.48 x 10 19 J. Then, since 
Energy = force x distance, we can roughly estimate the attractive force as energy/distance or 

F atoms = (5.48 x 10’ 19 J)/(2.55 x 10’ 10 m) = 2.15 x 10’ 9 N. 


Then finally, the force per unit area is 

Fatoms = (2.15 X 10’ 9 N)/(6.51 X 10’ 20 m 2 ) = 3.30 x 10 10 Pa = 33 GPa 


Note that this is comparable to the elastic modulus (E), not the tensile or shear strength , which is about 1000 times 
smaller. Why is the strength so much smaller than the elastic modulus? For a perfect crystal with no 
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defects, the above estimate would be appropriate. But real materials have defects in their crystalline 
structure. The strength of materials is determined mostly by the microstructural properties like the 
number of defects, the size of the “grains” (individual crystals), and the response of the defects to 
strain. This is why small amounts of additives (like adding carbon to iron to make steel) to a 
material and the details of how the material is processed (e.g. heat treating, rolling, etc.) affect its 
strength so much, but do not significantly affect other properties such as E, V (see below), density, 
etc. You’ll learn much more about this in MASC 310. 

As materials deform under tension, they become longer of course, but they also become narrower. 
The ratio between the change in diameter (d) of a cylindrical sample and change in length (L) is 

called Poisson's ratio (v), i.e. 

V = -(Ad/d)/(AL/L) Equation 28 

The minus sign is there because under tension AL > 0 (sample lengthening) but Ad < 0 (sample is 
narrowing). Note that the volume (V) of the cylindrical sample is L*7td 2 /4 before applying the stress, 
and (L+AL)*7c(d+Ad) 2 /4 after applying the stress. So 

(V+AV)/V = 1 + (L+AL)(d+Ad) 2 / Ld 2 

or, keeping only terms with one Ad or AL (not (Ad) 2 , (Ad)(AL), (AL) 2 etc.) 

AV/V = 2(Ad/d) + (AL/L) = -2(AL/L) [-(Ad/d)/(AL/L)] + (AL/L) 

= (AL/L)(1 - 2v) Equation 29 

function test. From the above equation, it is apparent that for a material to have no change in volume 
under stress, one would need 1 - 2v = 0 or V = 0.5. In reality most materials have V ~ 0.3, which 
means that their volume increases under tensile load (AL/L > 0). Certainly one would not expect V 
> 0.5, for this would imply the volume decreases under tensile load, and increases under 
compressive load (AL/L < 0) — not very likely! 

Pixample 

A V 2 inch diameter steel bar increases in length from 10 cm to 10.4 cm under an applied force of 
10,000 lbf. 

(a) What is the stress in the bar? 

Stress = force/area = 10,000 lbf / (tc( 0.5 in) 2 /4) = 50920 lbf/in 2 . 

(b) What is the strain in the bar? 

Strain = AL/L = (10.4 cm — 10 cm)/(10 cm) = 0.04 

(c) What is the change in diameter of the bar? 
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Note that the volume of the bar isn’t constant; to answer this question you’ll have to use 
Poisson’s ratio.) 


v = -(Ad/ d)/ (AL/L); if v ~ 0.3 for steel as in Table 2, then 
Ad = -vd(AL/L) = -(0.3)(0.5 inch)((0.4 cm)/(10 cm)) = 0.006 inch 

Shear forces 


Tension and compression are forces that act in the direction perpendicular to a particular imaginary 
plane cut through a piece of material. The force that acts parallel to a particular imaginary plane cut 
through the material is called the shear force (V) (why V? I dunno...). The shear stress (x) is the shear 
force per unit area, i.e. 

X = V/A. Equation 30 

To describe all of the stresses (normal and shear) in a material, we need something more than just a 
vector (which sufficed to describe forces.) For example, a two-dimensional object in the x-y plane 
has two components of tension or compression (call them O x and O) and two shear stresses (one 
each in the x and y directions; call them x and x ; usually these are just called x and x ). This is 
4 entities acting in the x-y plane; the composite effect of these 4 stresses cannot be represented in 
terms of a single scalar or vector. In three dimensions, the situation is even more problematic. As 
shown in Figure 13, a three-dimensional object will have three components of tension or 
compression (one each in the x, y and z directions) and six components of shear (in the x-y plane, in 
the x and y directions; in the y-z plane, in the y and z directions; and in the x-z plane, in the x and z 
directions). So stress is actually not a single scalar value or even a single vector but a 3 x 3 matrix 
called the Cauchy stress tensor: 







Equation 31 


Keep in mind that for each force shown in Figure 13, there is an equal and oppositely-directed force 



Figure 13. Diagram of normal (a) and shear (x) stresses in an imaginary, infinitesimally 

small cube of material 
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on the opposing side of the imaginary cube. This looks fairly complicated, and perhaps it is, but one 
saving grace is that, in order for the moments about a very small cube of material to sum to zero, 

one must have T xy _ x = T , T zx _ z = T xz _ x and T zy _ z = T yz _ y , in other word the matrix is symmetric. So there 
are only 6, not 9, independent stresses, or in two dimensions there are only 3 independent stresses, 
namely G x , a and T = T . Note the analogy between the number of independent stresses (3 in two 
dimensions, 6 in three dimensions) and the number of degrees of freedom of a free body (see Table 
1) which determines the number of constraint equations needed to enforce static equilibrium; this is 
not a coincidence at all but a necessary outcome of the number of degrees of freedom of an object. 

Principal stresses 

Note that normal stress is defined as the stress in the direction perpendicular to an imaginary plane 
and shear stress is defined as the stress in the direction parallel to that same imaginary plane — but 
how should that plane be chosen? For a simple shape like a cylinder it seems natural to define a 
plane parallel to the ends of the cylinder, but what about oddly shaped objects? Will the stresses be 
different depending on how one chooses the coordinate system? Will an object fail or not fail under 
stress depending on how one chooses the coordinate system? What is important to learn from this 
sub-section is that the magnitudes of both normal stress and shear stress are entirely 
dependent on the choice of coordinate system. 

It can be shown that the coordinate system (x, y, z) of Figure 13 can be rotated such that all of the 

off-diagonal terms (i.e. all the shear stresses x) are zero; these coordinates are called the principal 
directions and the corresponding stresses the principal stresses. Now we’re down to 3 independent 
stresses in this coordinate system. Furthermore, in the principal directions two of the three 
coordinates yield the maximum and minimum stress attainable from any rotation of the coordinates. 
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But proving this or using these results is beyond the scope of this course; wait for AME 204. In this 
course we will consider only the simpler two-dimensional case (Figure 14). If I know in some 

coordinate system (x,y) the normal stresses G x and G and the shear stress X x , then by rotating the 
coordinate system by an angle 0 P , the principal stresses (call them cq and G 2 , corresponding to the 
maximum (most positive or least negative) and minimum (most negative or least positive) stresses in 
the material) are obtained; their values are given by 



+ Oy 

2 



°X + Oy 
2 



n 1 
6 P = —tan 

! 2 


2r 


xy 


°x-^ y/ 


Equation 32 


Note that cq and 0 2 could be both positive, both negative, or one of each depending on the values 
of O x , a and X x . If both are positive then the larger one is the only one you need to worry about in 

terms of material failure, and failure could occur in tension only. If both cq and 0 2 are negative, 
then the more negative one is the only one you need to worry about in terms of material failure, and 
failure could occur in compression only. If cq is positive and C 2 is negative or vice versa, then you 
need to worry about both in terms of material failure, the positive one in tension and the negative 
one in compression. 

Also note that there are some good function tests you can perform on the formula for principal 
stresses: 


If x —0, then cq — O x and G 2 — a , that is, the normal stresses are the principal stresses since 
there is no shear in this coordinate systems. 

If o and x are both zero, that is, if there is only one normal stress and no shear stress, then 
cq — O x and 0 2 =0. 

If o = 0 but x i 1 0, then the principal stresses are 



which makes sense because adding the shear increases cq to a value larger than its value if 
there were no shear (cq), and decreases o 2 to a value smaller than its value if there were no 
shear (0). In other words, the addition of shear increases both the minimum and maximum 
normal stresses. 


Also, by rotating the coordinate system by a different angle 0 S , the maximum shear stress (t max ) is 
obtained; its value is given by 


T = 

max 


/ \ 2 

a -a 

x y 


2 


+ V = 


°i V 


2 


d<- = — tan 
s 2 


-1 


/ \ 

a -a 

x y 


2t 


xy / 


= e P ± 45 c 


Equation 33 
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In this coordinate system, the normal stresses are the same and equal to (a x + a)/2, i.e., the average 
of O x and O . 

Note that if in the above equations the stress in the x direction is non-zero (cf x /= 0) but the stress in 
the y direction is zero (o = 0) and the shear in the x-y plane is zero (x = 0), then the principal 
stresses are cq = o x , o 2 = 0 and T max = oj 2. Note also that just because I’m only pulling on the 
material in one direction (say, in tension) that doesn’t mean that there is no shear stress in the 
material; it’s all a matter of my choice of coordinates. The important conclusion is that a material 
under any type of stress has both normal and shear stresses; to determine the conditions for 
failure, it is not sufficient just to calculate the stresses in one particular coordinate system. 
One must determine the maximum normal and shear stresses in the material based on the above 

equations for cq, cq and T max and choose an appropriate dimensions and materials that can withstand 
such stresses. An analogy of sorts is with Alfred Hitchcock movies — typically the main character is 
an ordinary person doing some ordinary task, then something happens to him/her that causes 
him/her to become involved in some terrifying event. The message of his movies is typically, “you 
think you’re not involved... but you ARE.” The same thing applies to stresses: “you calculate 
normal stress and you think you’re not involved with shear stress... but you ARE.” Note that 
according to Eq. 33, the only situation where the material has no shear stress at all (T max = 0) is when 
a x = a y and t xy = O.c 

E xample 

A horizontal steel bar V 2 inch in diameter is pulled with a tension of 10,000 lbf, then a load is 
hung on it in such a way that the shear force is 5,000 lbf. 

(a) What is the maximum normal stress in the bar? 

Normal stress along length of bar = G x = Force/Area 
= (10,000 lbf)/(71(0.5 in) 2 /4) = 50,930 lbf/in 2 
Shear stress = x xy = (Shear force)/Area = (5,000 lbf)/( tc( 0.5 in) 2 /4) = 25,465 lbf/in 2 

From the above equation for principal stresses we find that cq, o 2 are 

r y = 0, x xy = 25,465 lbf/in 2 
,-10,548 lbf/in 

So the maximum normal stress is +61,477 lbf/in 2 (+ sign indicating tension, - sign would 
indicate compression). Also note that the transformation from the coordinate system x, y to 
the principal directions requires an angle of rotation given by 
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' 2(25,465 Ibf I in 2 )' 
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22.5° 


(b) What is the maximum shear stress in the bar? 





+ 25,465 2 


±36,013 lbf/in 


Pressure vessels 

In a cylindrical vessel containing a pressure P, there are 2 stresses to be considered: the hoop stress 
trying to pull the cylinder apart radially, and the longitudinal stress or axial stress trying to pull it apart 
axially. Referring to Figure 15, note that the total force trying to pull the cylinder apart radially is PA, 
where the area A = 2rL, where r is the cylinder radius and L its length. The total wall cross section 
area resisting this force is 2xL, where T is the wall thickness (not to be confused with T the shear 
stress used above; wall thickness has units of length, shear stress units of pressure). Thus, 


Hoop stress (o^) — (total force)/(area of wall resisting force) 

= (P * 2rL) / (2 tL) = Pr/x Equation 34 

Similarly for the longitudinal stress, the total force trying to pull the cylinder apart axially is PA = 
2 

P(7ir“) and the total wall cross section area resisting this force is 271 * r * X, thus 


Longitudinal stress (Oj) — (total force)/ (area of wall resisting force) 

= P(7ir 2 )/(27crx) = Pr/2x Equation 35 

Note that the hoop stress is twice the longitudinal stress. This is why an overcooked hot dog usually 
cracks along the longitudinal direction first (i.e., its skin fails from hoop stress, generated by internal 
steam pressure). Also note that both hoop and longitudinal stress are both positive, i.e. in tension if 
the pressure inside the vessel is higher than that outside the vessel. Of course, if the pressure 
outside is higher (i.e. a vacuum chamber or submarine) then both will be hoop and longitudinal 
stress are both negative, i.e. in compression. 



Figure 15. Diagram of pressure vessel showing hoop stress (a h , left) and longitudinal stress 
(o 15 right). Figures from http://www.efunda.com/ . 
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So for the pressure vessel we have hoop stress (call it G x , where x is the radial direction) — Pr/x and 

longitudinal stress (call it G, where y is the axial direction) = Pr/2x. In this coordinate system there 
is no shear stress. Thus the principal stresses are 


O l 9 a 2 


°x + °y 


2 



Pr/r + Pr/2r 
2 



Pr Pr 
r ’2r 


So for this case the principal stresses G t and G 2 are just the calculated stresses in the x and y 
directions; in fact this will happen any time x = 0. On the other hand, the maximum shear stress 
for the pressure vessel is 


max 




Pr/r -Pr/2r\ 2 

2 ) 




Equation 36 


Thus, unless the yield strength in shear was less than Y 4 the yield strength in tension or compression, 
the material would fail in tension or compression before it failed in shear. 

.Example (to be continued below ..J 

An iron pipe 1 foot in diameter and 50 feet long has a wall thickness of 1/2 inch. The material 

properties are: elastic modulus (E) = 30 x 10 6 lbf/in 2 , yield stress (O yield ) = 30 x 10' lbf/in 2 in 
tension , yield stress = -30 x 10 4 lbf/in 2 in compression and yield stress = 10 x 10 3 lbf/in 2 in 
shear . If the ends of the iron pipe are sealed and the pipe is used as a cylindrical pressure vessel, 
with the high pressure inside: 

(a) At what pressure (in units of lbf/in 2 ) will the iron yield? 


Maximum stress — hoop stress (G h ) — +Pr/x (tension) 

Pressure at yield = G y x/r = (30 x 10' lbf/in 2 )(0.5 in )/(6 in) = 2,500 lbf/in 2 



For what yield stress in shear (in units of lbf/in 9 ) will the iron yield in shear rather than in 
tension for this pressure? 


T_ = ± 


max 


Pr 

4r 


p r (iSOOIbf /in 2 )(6in) 


'yield 


= 7500 lbf/in 


4t 4(0.5 m) 

Since the actual yield stress in shear is 10,000 lbf/in", the pipe will not yield in shear at this 
pressure and thus it will yield in tension instead as calculated in part (a). 


Bending of beams 

One of the most common problems in structural mechanics is the compute the stresses in a beam 
subject to a load, perpendicular to the axis of the beam, distributed over the length of the beam. 
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Figure 16. Left: Force balance on a differential element of beam of length dx, showing the 
applied load per unit length w, shear force V and bending moment M. From 
http:/7www.efunda.com/ . Right: Schematic of compressive and tensile stresses in beam 
caused by bending moment. From http://strengthandstiffness.com/6_beams/page_6b.htm . 


The load is typically reported as a force per unit length along the beam (w), with units N/m or 
(more likely) lbf/ft. As shown in 

Figure 16 (left), this load causes a shear force in the beam (V) = wL, where L is the distance from 

the end of the beam. Or, if w is not constant, we can say that dV = w dx and V = f Q W dx. Then 

the moment about one end of the beam is given by (defining counterclockwise moments are 
positive, which is standard in structural mechanics) dM = -Vdx, thus 


dM 

dx 




d 2 M 
dx 2 


Equation 37 


Knowing w(x) for the beam and the boundary conditions at the ends of the beam, one can 
determine the moment M at any location along the beam. For example, with a constant load w per 
unit length and which is pinned at one end (able to withstand a force in both x and y directions, but 
unable to cause a moment) and has a roller at the other end (so able with withstand a force in the y 
direction only, so there are 3 unknown forces and 3 degrees of freedom, i.e. a statically determinate 
system) we have 
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Figure 17. Distribution of bending moments along a beam for uniform or point loading 

d"M/dx = -w = constant 

with the boundary conditions 

M = 0 at x = 0 and x = L 

for which the solution is 

dM/dx = -wx + q; 

M = -wx 2 /2 + qx + c 2 ; 

M = 0 at x = 0 => c 2 = 0; 

M = 0 at x = L =^> -wL 2 /2 + q(L) + 0 = 0 => q = wL/2 

=> M(x) = (wx/2)(L — x) (uniformly loaded beam) Equation 38 

Note that the maximum of M is at x = L/2 with a value of wL 2 /8. 

Who cares about these moments? Well, what we do care about is the stress in the beam. The beam 
must resist this moment by the stresses in the material. In order to do that, the beam has to be (for 
downward loading, i.e. weight on the beam) in tension on the bottom and compression on top. In 
other words 
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Forces (loads) on beam => Bending moments in beam => Stresses in beam material 


It’s important to understand that the bending moments cause far more stress in the beam 
than the stresses caused by the direct application of the force w(x). It’s beyond the scope of 
this course to derive the relationship between bending moment M(x) and stresses (you’ll learn about 
this in AME 204), but for a slender beam (one for which its length L is much larger than its height in 

the y direction) by far the largest stress in the beam is the normal stress in the x direction (cy) 
resulting from this bending moment and this stress given by 


& x (x,y)=- 


M(x)y 

I 


Equation 39 


where y is the vertical distance from the “neutral axis” of the beam (where cr x = 0) (the “neutral axis” 
is half-way through the beam for a symmetrical cross-section), M(x) is the moment just computed, 
and I is the moment of inertia of the beam cross-section (units are length 4 ). The moment of inertia 
about an axis A-A’ is defined in Figure 18. In comparison to O x (x,y) given by this equation, the 
normal stress in the y direction (a y ) and the shear stress T in the (x,y) coordinate system are much 
smaller and so can be neglected. (Of course we still have to find the maximum shear stress T max 
obtained by rotating the coordinate system, as discussed below.) 


Note that for a given total cross-section area (thus total weight of beam) one can have large I (thus 
lower stress o) by having more material at larger distances from the axis A-A’. This is the reason for 
using I-shaped beam sections. Formulas for I for common shapes include: 


Circular cross-section of diameter d: I = Ted 4 /64 

Thin-wall hollow tube of diameter d and wall thickness T: I = Tcd’x/8 

Rectangular cross-section I = ab Vl2 (a = width of beam; b = height of beam) 

I beam of width a, height b, thickness of central section T w and thickness of top and bottom 

atf (a-r w )(b-2r h f 


sections T h : I — / = 


12 


12 


Note that the I-beam formula satisfies the function tests: when T w — a or T h — b/2, the I-beam is 
just a “filled” rectangle with I = ab Vl2, and if T w = a and T h = b/2, the beam has no material thus I 



Figure 18. Definition of moment of inertia (I) about axis (dotted line) A-A 5 
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0. 


Now let’s use this for a simple case of a uniformly-loaded or point-loaded rectangular cross-section 
beam (Figure 19) of height b and thickness a. For this case, the moment of inertia I = ab 3 /12, thus 
G max (at the center of the beam (x = L/2), at the top or bottom, where y = +b/2 at the top of the 
beam and —b/2 at the bottom of the beam) is given by 


M y 

.y __ max J max 

max 


/ 


( ^/ 8 )(± < , / 2 ) , t075 ^ 

ab~ /12 ab 


(uniform loadl 


nr 


AT 


/ ^ ^ # 4 \ / 


1 1 


Equation 40 


where the — sign refers to the compression at the top of the beam and the + sign refers to the 
tension at the bottom of the beam (recall that the sign convention for stresses is that compression is 
negative and tension is positive.) Does this result make sense? 

• Smoke test: The units (for uniform load) are wL 2 /ab 2 = (Force/Length)(Length) 2 /(Length’) 

= Force/Length - , which is stress — OK. 

• Function test #1: a longer (larger L) or more heavily loaded (larger w or P) beam should have 

more stress - OK 

• Function test #2: a thicker (larger a) or taller (larger b) should have less stress — OK 


Note that as mentioned below Eq. 33, for a situation in which there is normal stress O x in only one 
direction (a = 0) and no shear stress in that x-y coordinate system (t = 0), there is still a shear 
stress T max = oj 2 according to Eq. 33. So the beam would fail in shear if the yield strength in shear 
is less than half of the smaller of the yield strength in tension or compression. 

Also note that for a given total load (in units of force) = wL (uniform load) or P (point load) and 
beam length L, the stress is proportional to 1 /ab 2 , whereas the weight is proportional to volume = 
abL. Thus to minimize the stress for a given weight of beam, one wants to minimize the ratio 
(abL/ab 2 ) = L/b, meaning that (since L is already fixed) we want to maximize b (and thus minimize 
a) — in other words, a tall skinny beam cross-section works better than a short fat one. That’s 
another reason for using the I-beam shape. Another way of thinking of this is that since O max is 
proportional to l/ab“, one gets more benefit from increasing b than increasing a. Increasing a 
results in a proportional decrease in stress and a proportional increase in weight, so the strength-to- 
weight ratio doesn’t change. However, increasing b resulting in a more-than-proportional (l/b 2 ) 
decrease in stress, thus the strength-to-weight ratio increases. 
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Figure 19. Schematic of example beam-loading problem 


By examination of Fig. 19, one might notice that there is a compressive stress acting on the beam by 
virtue of the loading. In Fig. 19, this “direct” compressive stress would be Force/Area = wL/aL = 
w/a. Is this a lot or a little? By comparison the compressive stress caused by the bending moment 
is (Eq. 41) 0.75wL/ab. The ratio of the bending-induced to direct compression is then 
(0.75wL 2 /ab 2 )/(w/a) = 0.75(L/b) 2 . Since we have already assumed L >> b (a long, slender beam), 
the bending-induced compression far exceeds the “direct” compression. Similarly, for the shear 
force V caused by the “direct” loading, from Eq. 37 dV/dx = -w, thus V = -wx = -wL/2 at the 
center of the beam. This results in a shear stress = force/area = -(wL/2)/ab. The ratio of the shear 
stress caused by the bending moment = oj 2 = (0.75wL 2 /ab 2 ) = 0.375wL 2 /ab 2 to that caused by the 
“direct” loading is then 0.75(L/b). So again, since L >> b, the stress caused by the bending 
moment is much more than that due to direct loading, which explains why we can usually ignore the 
direct loading when determining the point at which a beam will fail. 


Another property of some interest is the maximum deflection A (i.e. the sag in the middle of the 
beam) due to the applied load. Its value is given by 


A = 


5 wU 
384 El 


(uniform load); A = 


PL 

48 El 


(point load) 


Equation 41 


Note that for the same total applied load (wL — P), the maximum deflection is (1/48)/(5/384) — 1.6 
times larger for the point load than the uniform load. 

Another common stress analysis problem is a circular disk (e.g., the top or bottom of a cylindrical 
pressure vessel) of radius r and thickness x with pressure P on one side. The maximum stress 
(including the transformation to principal stresses) is given by 

Pr 2 

a max = ±c ~y Equation 42 

x 
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where c = 1.24 if the edges of the disk are free to pivot (e.g. like a drum head, which is not very 
realistic) or c = 0.696 if the edges are rigidly clamped and unable to pivot (which would be the case 
if the disk were welded or bolted on to the end of the cylindrical part of the pressure vessel.) 

Example 

(a) If the iron pipe from the pressure vessel example above has no pressure inside but instead is 
used as a beam with one pinned end and one roller end instead, what is the maximum point 
load (P, units of force, not to be confused with pressure P) that could be applied at the 
middle of the beam before the iron yields? 


(PLIA)(-dl2) PL 
I Ttd^r / 8 itd 2 r 

Jid. 2 xo yield _ ?t(12in) 2 (0.5iri)(30xl0 i lbf /in 2 ) 
L (50 ftx\2in / ft) 


11,304 Ibf 


• • • • • • ^ 
(Note: the material is weaker in tension, so the failure will occur when a = +30 x 10' 

lbf/in 2 in tension at the bottom of the beam (y = -d/2) rather than a = -30 x 10 4 lbf/in 2 
in compression at the top of the beam (y = +d/2)). 

However... we need to check for failure due to shear stress also. At the top or bottom 
of the beam, for a given value of the normal stress in the x direction (along the length of 

the beam) of G x , assuming no normal stress in the y direction (a = 0) and no shear 
stress (x = 0), there is a shear stress given as (Eq. 33): 


r 


max 





2 


As discussed below Eq. 33, in this particular case the maximum shear stress is half as 
much as the normal stress O v . Thus 

A 

_ I Vnax.Vax \ (PL ! Ar)(-d !l) PL 

yie,d 2 1 2 jtd 3 t/8 2jtd 2 T 

p _ 2jtd 2 r(r yield ) _ 2^(12m) 2 (0.5m)(10xl0 3 /(7//m 2 ) 536 

L (50 ft x 1 2in / ft) ’ 1 


So the beam yields at a lower point load P in shear than in tension (7,536 lbf vs. 11,304 lbf.) 
(b) What is the maximum deflection of the beam? 



PL 3 
48 El 


(U,300lbf)[(50ft)(l2in/ft)] 3 
48(30 x \(f lbf / in 2 )(jt(\ 2m) 3 (0.5m)/8) 


= 5.00m 
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Notice that the 50 foot long beam deflects/bends a maximum (i.e. at the failure load) of 5 
inches, i.e. only 5 / (50 x 12) = 0.0083 = 0.83%. 

(c) Assuming that shear stress is not an issue (as it was shown to be at the end of part (a)), how 
much compressive force (pre-stressing) should be applied to the pipe to maximize the point 
load that could be applied? What would this maximum point load be? 

At the maximum load condition there is G x = -30 x 10' lbf/in (compression) at the top 
of the beam and O x = +30 x 10' lbf/in (tension) at the bottom of the beam. This is sort 
of a waste, since the top of the beam could take a lot more compressive stress before it 
failed. So by pre-compressing the beam, one could even things out. So to have both the 
top and bottom of the beam at their maximum stress, we would have a stress of -30 x 
10 4 = PC — S at the top of the beam (- sign indicating compression, PC is the pre¬ 
compression, S is the stress due to the applied load — on top, + on the bottom) and +30 
x 10 = PC + S at the bottom of the beam. So combining this two equations, we have 

-30 x 10 4 = PC - S 
30 x 10 3 = PC + S 

=> -27 x 10 4 = 2 PC => PC = -13.5 x 10 4 ; S = 16.5 x 10 4 

We already showed that without pre-compressing, a point load of 11,300 lbf produces a 
stress of +30 x 10' lbf/in“ in the beam. With pre-compressing, we can withstand ±16.5 
x 10 4 of stress in the beam due to the loading, that is, 5.5 times more stress. Since the 
relationship between the applied load P and the stress is linear, we can conclude that P = 
5.5 (11,300) = 62,200 lbf. 

(d) If the pipe has welded disk end caps of the same material and thickness as the pipe, at what 
pressure P (again this is pressure P, not to be confused with point load P) would the end cap 
fail? 


a =±c 

max 


PP 

~2 


cr_T z (30x 10'lbf / in 2 )(0.5in) ,,..2 


P = ± 


max 


cr 


= + 


0.696(6m) ; 


= 299 lbf / in 


where again the tensile (not compressive) strength is chosen because it the smaller value. 


Buckling of columns 

Another way in which a structural element under compression can fail is by buckling. This is not 
strictly a failure of the material, but effectively eliminates the load-carrying capability of the structure. 
The compressive force (F) at which buckling occurs in a column of length L is given by 

2 2 

Fbuckling = njTEI/L Equation 43 
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where E is the elastic modulus discussed above, I is the moment of inertia of the cross-section of 
the column in the plane perpendicular to the direction of the applied force (which is parallel 
to the long direction of the column) and n is a constant that depends on the way in which the 
column ends are or are not held: 

• Both ends pinned, i.e., free to pivot: n = 1 

• Both ends clamped, i.e., unable to pivot: n = 4 

• One end pinned, one end clamped: n = 2 

It should be noted that this buckling formula is valid only for a “slender” column (where the length 
L is much greater than the width of the column cross section) and it assumes that the column cross- 
section does not change (in other words, it would not account for the crumpling of an aluminum 
beverage can, where the buckling occurs due to a change in the cross-section of the can. 

When computing I for a rectangular cross-section of a buckling column, which is the “a” dimension 
and which is the “b” dimension? Since the column can buckle either way, you have to use the lesser 
I, i.e. where a is the smaller dimension, which says that to avoid buckling, you don’t want tall skinny 
beam cross-sections, you want round or square ones. Note that this conflicts with the desired cross- 
section to minimize stress due to bending moments, i.e. it was just mentioned above that for best 
strength to weight ratio, you want tall skinny I-beams. Thus the optimal I-beam cross-section will 
be a compromise between the two shapes. 

'Example 

What is the buckling load of a polyethylene plastic drinking straw (E ~ 10 9 Pa), W’ in diameter and 
1/32” wall thickness, 6” long, with both ends free to pivot? 

^buckling — EI/L; 

n = l; 

E ~ 10 9 Pa; 

I = 7td 3 x/8 = tc( 0.25) 3 (1 /32)/8 = 1.92 x l(h 4 in 4 = 7.99 x 10 -11 m 4 
L = 6 in = 0.152 m 

F buckling = 1 tt 2 (10 9 Pa)( 7.99 x 10' 11 m 4 )/( 0.152 m) 2 = 34 N = 7.7 lbf. 

In practice the buckling load would be less because this analysis assumes the load is exactly along the 
axis of the column, whereas in reality there would be some sideways (shear) load. 


55 



Chapter 7. Fluid mechanics 


“The goal in life is to be solid, whereas the way that life works is totally fluid, so you can never actually achieve that 
goal Damien Hirst (British artist). 

Main course in AME curriculum on this topic: AME 309 (Dynamics of Fluids). 

Fluid statics 


Fluid mechanics is basically 2F = d(mv/dt) (Newton’s 2 nd Law, the sum of the forces is equal to the 
rate of change of momentum) applied to a fluid. What distinguishes a fluid from a solid is that a 
solid deforms only a finite amount due to an applied shear stress (unless it breaks), whereas the fluid 
continues to deform as long as the shear stress is applied. This makes fluid mechanics a lot more 
complicated (at least to me) than solid mechanics. 

Hydrostatic pressure 

Let’s look first at fluid statics , i.e. when 2!F = 0. If a fluid is not moving at all, as in a glass of water, 
then the fluid is static , that is, zero velocity everywhere, and has only a hydrostatic pressure. Imagine a 
column of water of height z, cross-section area A and density p (units M/L 3 , i.e. mass of fluid per 
unit volume.) 

Table 3 gives the density of several common liquids and gases. The weight of the water is the mass 
x g = density x volume xg = pzAg. This weight is distributed over an area A, so the force per unit 
area (the hydrostatic pressure) is pzAg/A = pgz. This is added to whatever pressure P(0) exists at z 
= 0. So the hydrostatic pressure P(z) is 

P(z) = P(0) - pgz Equation 44 

where z is defined as positive upward, i.e. decreasing depth. This result assumes that the density (p) 
is constant. This is reasonable for water and practically all liquids, even at pressures of thousands of 

atm. It’s also ok for gases if z is not too large, i.e. such that pgz << P G . 


Fluid 

Density 

(P, kg/m 3 ) 

Dynamic viscosity 
(H, kg/m sec) 

Kinematic viscosity 
(v = ju/p, m 2 /sec) 

Water 

997.1 

8.94 x UT 4 

8.97 x 10~ 7 

Air 

1.18 

1.77 x UT 5 

1.50 x 10 5 

Motor oil 

917 

0.260 

2.84 x 10~ 4 

Mercury 

13500 

1.53 x UT 3 

1.13 x 10~ 7 


Table 3. Properties of some common fluids at ambient temperature and pressure. 
Buoyancy 
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According to Archimedes’ principle, an object of volume V placed in a liquid of density p will exert 
a buoyant force equal the weight of the fluid displaced = p f gV. The net force on the object is the 

difference between this Archimedean (buoyant) force and the weight of the object = p o gV, where p () 
is the average density of the object (just total mass/total volume). Thus the net force F acting on the 
object is 


F = (p f - P„)gV 


Equation 45, 


where the sign convention is such that the force is positive (directed upward) when the object 
density is less than the fluid density (i.e., the object floats upward). ¥ unction test if the density of the 
object and the fluid are the same, the object is “neutrally buoyant,” and there is no net force on the 
object (F = 0). 


Example 


a) The deepest part of the ocean is a spot called “Challenger Deep” in the Marianas Trench in 
the western Pacific Ocean. The depth is 35,838 feet. The density of seawater is 1026 kg/m 3 . 
What is the hydrostatic pressure (in atmospheres) at this depth? Remember, at sea level, the 
pressure is 1 atm and increases as the depth increases. 


Ocean depth: 2 —35838 ft —10923 m; seawater density p — 1026 kg/m’ 


P = P( 0) - pgz = latm - (l026kg / m 3 j(9.81m/s ,2 )(-10923m) 


1 atm 

101325 N/m 2 


= I086atm 



The density of air at sea level is 1.18 kg/m’. If the air density were constant (not a function 
of elevation), at what elevation would the pressure be zero? 


P = P( 0) - p air gz = 0 => P{ 0) = p air gz thus z = - 


m 

Pair£ 


101325 N/m 2 
lASkg/m 3 x9.81m/A 2 


= 8732m 


c) Until 2012, the only vessel ever to carry people to Challenger Deep was the bathyscaphe 
Trieste in 1960. It used gasoline (p = 739 kg/m’) for flotation since no air tank could be 
made light enough to sustain an 1100 atm pressure difference and still provide positive 
buoyancy (since the gasoline is essentially incompressible it could be contained in a thin- 
walled tank that did not need to sustain a pressure difference between the gasoline and the 
surrounding seawater.) The Trieste used 22,500 gallons of gasoline for flotation. How much 
buoyant force could this much gasoline produce? 

In this case the surrounding fluid is seawater and the “object” is the gasoline itself, thus 
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v ’ \ m m / sec 


22,500 gal 


fi 


7.46 gal 


/ \ 

3\ 

m 


13.281 ft) 

/ 


F = 2.404x10 


5 kg m ^ A riA ..ir\5*T 1 Ibf r /f i . . i rv4 


= 2.404 x 10 


sec 


4.4487V 


= 5.41 xlO*Ibf = 27.0 tons 


Equations of fluid motion 


Bernoulli’s equation 

One of the most common problems in fluid flows is to determine the relationship between velocity, 
pressure and elevation of a flowing fluid in a pipe or other duct. To do this, we enforce conservation of 
energy on the fluid, i.e. the energy contained by the fluid at one point in the flow is the same as any 
other, but the energy may be transformed from one form to another. Moreover, it is more 
convenient to work with power (rate of change of energy) rather than energy itself. If the flow is 
steady so that no energy is accumulating or dissipating within the pipe then the power (sum of all 
forms) must be constant. This course in general is not intended to provide derivations of formulas 
you will study in much greater detail in later courses, but it is worthwhile to do so for Bernoulli’s 
equation just as an example of the value and power of units, and the concept of conservation (e.g., 
of energy) applied to a fixed volume (often, a fixed mass is analyzed rather than a fixed volume). 

There are 3 types of power that must be considered, and their sum conserved. In words, the 
conservation of energy can be stated as: 

(Power needed to push fluid into the tube inlet - power extracted at the tube outlet) 

+ (kinetic power of the fluid flowing into tube - kinetic power of the fluid flowing out of tube) 

+ (power associated with change of gravitational potential energy of fluid) 

= 0 


Let’s compute the individual terms then add them up. 

1. Power needed to push fluid into the tube inlet or power extracted at the tube outlet: 


Power = force x velocity = [ ^°~ C ^ {velocity x area ) 


area ) 

= pressure x (volume/time) = pressure x 


mass / time Pm 


(Equation 46) 


mass/volume p 


where m is the mass flow rate (units kg/sec), discussed in more detail in the next sub-section. 

2. Kinetic power of the fluid: 


„ d (1 l dm 2 1-2 

Power = — — mv =-v =—mv 

dt \2 2 dt 2 


(Equation 47) 


3. Power associated with gravitational potential energy of fluid 
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(Equation 48) 


Power = 



dm 

dt 


gz = rhgz 


Combine : sum of the powers at inlet (call it station 1) = sum of powers at outlet (call it station 2). 
Assume mass flow rates are equal at inlet and outlet (if they’re not, the flow can’t be steady because 
mass will be accumulating or being lost from the pipe) 

Pm, 1.2 • 9 1. 2 

+ ~ m i V l + = JLJL + ~ m 2 V 2 + m 2^2 

A 2 Pi 2 

1 2 1 , 

+ -AT + m = ^2 + — P V 2 + Pg*2 


(Equation 49) 


This is Bernoulli’s equation, which is merely a statement of conservation of energy for an 

incompressible (p = constant), inviscid (no viscosity), steady, one-dimensional flow between 

locations 1 and 2. Recall that the term pv 2 /2 is called the dynamic pressure , i.e. the increase in pressure 
that would occur if the fluid were decelerated (at constant z) from velocity U to a velocity of zero. 
¥ unction test it was shown that for a static fluid (v 1 = v 2 = 0), P(z) = P(0) - pgz, or P(z) + pgz = P(0), 
which is the same as Bernoulli’s equation for z 2 = 0. 

If there are more than one inlets or outlets, we still have to conserve energy, thus the sum of the 
Bernoulli terms must be the same at the inlet and outlet. For example, if there are two inlets (say la 
and lb) and two outlets (say 2a and 2b) then 


1 


1 


1 


1 


Pa+-p V \a+PSZ la + P lb + ~ PVlb + PgZ lb = Pla + - P V 2a + PgZla + P 2b + ~ P V 2b + PgZ 2 b 


Note that Bernoulli’s equation assumes that the density (p) is constant. At first glance this might 
suggest that it cannot be used for air or other gases, which are compressible. Actually, Bernoulli’s 
equation can be used for gases if the Mach number (ratio of velocity to sound speed) is significantly 
less than 1. This applies to most of our common flow situations, e.g. for all practical purposes the 
air flowing over a car can be considered to have constant density, as discussed later. 

Besides the assumption of steady ID flow, constant density, what other significant limitations does 
Bernoulli’s equation have? The most important is that friction losses (viscosity) are not 
considered. Can we just add another term to Bernoulli’s equation to account for viscosity? No, 
because viscosity is dissipative and causes a loss in the total power (sum of the three terms). Where 
does the power go? Into thermal energy of the fluid (i.e. it gets hotter.) Moreover, the amount of 
power lost is path dependent , i.e. a longer or narrower tube will have more loss, whereas the above 
three terms don’t depend on the length or diameter of the tube. Of course, viscosity can be 
incorporated into fluid flow analysis, but it’s much more difficult and can’t be done with a simple 
equation like Bernoulli’s that depends only on the initial and final states (1 and 2 in our notation) 
and not at all on the path between states 1 and 2. 
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Conservation of mass 


How does one determine the velocity v? For steady flow in a pipe, channel or duct, the mass flow 
rate m (in kg/s) has to be the same everywhere in the flow system. This mass flow rate is the 
product of the fluid density (p), velocity (v) and the cross-section area (A) of the pipe or duct 
through which the fluid flows, i.e., 

m = PjVjAj = p 2 v 2 A 2 Equation 50. 


In the case of Bernoulli’s equation, we have already assumed that the density is constant, so for this 
case (i.e. liquids as well as gases at low Mach number) we can simplify this to 

Mi = M2 

note that UA has units of (length/time)(length) 2 = length 3 /time = volume/time, i.e. the volumetric 
flow rate , usually given the symbol Q. Thus for an incompressible fluid = Q 2 . If there are more 
than 1 inlets or outlets then the sum of the mass flows at the inlets must equal those at the outlets, 
i.e. 


P\a V \a\a + P\b V \b\b Pla V 2a^2a + P2b V 2b^2b 
or for incompressible flow (p = constant) 

+ V lb^lb = V 2 Afl + V 2b^2b OT Q\a + Q\b = Qla + ^2/? 


Example 

Water flows from a faucet at elevation 2 = 0 with supply pressure (P t ) of 30 lbf/in 2 = 
207,000 Pa (above atmospheric), area 5 cm", to the roof of a house with 2 = 5 m and 
through a no 22 le with area 1 cm , into ambient air with a pressure (P 2 ) of 0 lbf/in - above 
atmospheric. What is the velocity of the water leaving the no 22 le? Assume that viscous 
effects are negligible. 

Neither v t nor v 2 are known, but P 1? P 2 , 2 1 and 2 2 are all known, so we have 2 equations 
(Bernoulli and mass conservation) for the two unknowns. 

First apply mass conservation: v 1 = v 2 A 2 /A 1 = v 2 (l cm 2 )/(5 cm 2 ) = 0.2 v 2 . 

Then use Bernoulli and solve for v 2 : 
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| P (■vl - vf ) = (Pi - p 2 ) + pg (z, - z 2 ) 
v 2 2 - (o.2v 2 )“ = 2 [(/> - P 2 )/p + g(z! - z 2 )] 
v 2 =^(2/0.96)[(P l -P 2 )/p + <? (z 1 -z 2 )] 

= ^(2/0.96)[(207,000iV/m 2 -0)/(1000^/m 3 ) + (9.81m/5 2 )(0-5m)] = 18.l7«/5 
Note also that velocity at the supply faucet v 1 = 0.2 v 2 = (0.2)(18.1 m/s) = 3.63 m/s. 


Viscous effects 


Definition of viscosity 

Bernoulli’s equation pertains only when there is no viscosity (i.e. the flow is inviscid). Fluids resist 
motion, or more specifically resist a velocity gradient , through viscosity (p), defined by the relation 
(called Newton’s Law of Viscosity) 

d\ _ 

T v _ = fl —- Equation 51 

dy 


where T xz is the shear stress in the x-z plane, v x is the component of velocity in the x direction, and 
dv/dy is the velocity gradient in the y direction. This type of viscosity is called the dynamic viscosity. 
Since T xz has units of force/area = (ML/T 2 )/L 2 , u has units of L/T and y has units of L, the 
viscosity p has units of M/LT, e.g. (kg/ms). This unit has no particular name, but 1 g/cm s = 0.1 
kg/m s = 1 Poise. The unit centipoise — 0.01 Poise = 0.001 kg/ms is frequently used because the 
dynamic viscosity of water at ambient temperature is almost exactly 1 centipoise. 

Another type of viscosity is the kinematic viscosity , which is just the dynamic viscosity divided by 
density: 

V = — Equation 52 

P 

which has units of (M/LT)/(M/L 3 ) = L 2 /T, e.g. m 2 /s. Again, this unit has no particular name but 1 
cm 2 /s = 10’ 4 m 2 /s = 1 Stoke, and again 1 centistoke = 0.01 Stoke, which is very nearly the kinematic 
viscosity of water at ambient temperature. The units of kinematic viscosity, L 2 /T, are the same as 
that of diffusion coefficients, e.g. the property that describes how fast a drop of ink will spread out 
in a beaker of water, thus my favorite interpretation of V is that is it the momentum diffusivity. This 
describes how quickly or slowly the momentum of the fluid is exchanged with the solid object 
passing through the fluid (or fluid passing through the solid object, as in the case of flow in a pipe.) 
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No-slip boundary condition 

At the boundary between a fluid and a solid object, the velocity of the fluid and the solid must be 
the same. This is called the no-slip condition. This means that any time a solid is moving through a 
fluid (e.g. an airplane flying through the air) or a fluid is moving through a solid (e.g. flow through a 
pipe) there will be a velocity gradient (i.e. du/dy in the above equation) because there is a difference 
between the fluid velocity far from the boundary and the fluid velocity at the boundary, and this 
velocity difference occurs over some finite distance. This velocity gradient is the source of the 
viscous drag — the velocity gradient creates a shear stress on the fluid (see definition of viscosity 
above) that resists the motion of the fluid. It should be noted\ however ; that while the action of viscosity always 
causes drag, not all drag is due to viscosity. 

Reynolds number 


How important is viscosity in a given flow? That depends on the dimensionless quantity called the 

Reynolds number (Re): 


~ pvL vL 

Re = —— = — 

p 


(Equation 53) 


v 


where L is a characteristic length scale of the flow which has to be specified. (Note that the symbols 
“v” (Roman letter Vee’) for velocity and “v” (Greek letter c nu’) for kinematic viscosity are similar, be 
careful!) In the case of a wing, L is usually chosen to be the length of the wing in the streamwise 
direction (which is called the cord of the wing.) For flow in a pipe, L would be the pipe inside 
diameter. For flow around a cylinder or sphere, L would be the diameter of the cylinder or sphere. 
Also, the fluid velocity v changes as the fluid approaches the object, so v is chosen to be the value 
far away from the object. For flow inside pipes, v is the average velocity of the fluid, i.e. the volume 
flow rate (gallons per minute, m f sec) divided by the cross-section area of the tube. 

The standard catechism of fluid mechanics states that “Reynolds number is the ratio of inertial 
forces to viscous forces” but this is nonsense. First of all there is no such thing as “inertial forces” 

in mechanics. Second, pvL is not a unit of force, nor is p. Here’s my interpretation of Re. Re-write 
Re as 



pv 2 Dynamic pressure 

-/—-/- 

pv / L Shear stress due to viscosity 


or 



L /v Viscous diffusion time scale 
- /—-/ - 

L / v Flow time scale 


The first interpretation notes that pv 2 /2 is the dynamic pressure noted above in the context of 
Bernoulli’s equation and the velocity gradient dv/dy is proportional to v/L, thus the shear stress T ~ 
pv/L. The second interpretation (which is my personal favorite) notes that the time scale for any 
type of diffusion process is L 2 /D where D is the diffusion coefficient for that process (the viscous 
diffusivity V in this case), and the time scale for the fluid to move a distance L is simply L/v. So the 
second interpretation states that the Reynolds number is the ratio of the time for the momentum (or 
lack of momentum, as in a stationary wall with a fluid moving past it) to diffuse across a distance L 
to the time for the fluid to move a distance L. 


62 



Why is Reynolds number useful? Besides determining how important viscosity is in a given flow, it 
allows one to employ scaling. For example, suppose you want to determine the drag coefficient C D 
(another dimensionless number, to be discussed shortly) on a 5 meter long car at 30 m/sec (about 
67 mi/hr) and you don’t have a wind tunnel large enough to put a real car it in, but you have a water 
channel that is big enough for a 1/5 scale (1 meter long) model of a car. The kinematic viscosity of 
air at ambient temperature and pressure is about 1.5 x 10" 5 m 2 /s and that of water is 1.0 x 10’ 6 m 2 /s. 
Then by choosing the velocity of water in the water channel to get the same Reynolds number, you 
can obtain a valid measurement of C D : 


^air^air 


v L 

water water 


(30 m/s)(5 m) _ (v water )(l m) 


V 

air 


V 


1.5 x 10 5 rn / s 1.0 x 10 ° m z / s 


-6 


water 


V water = 10 m/s 


So a 5 meter long model in air moving at 30 m/s will have the same behavior as a 1 meter long 
model in water moving at 10 m/s. There may be other advantages to using water, e.g. the use of 
fluorescent dye molecules that make it easier to visualize the flow using a sheet of laser light. 

Navier-Stokes equations 

As previously stated, fluid mechanics is just F = dp/dt = d(mv)/dt applied to a fluid (here p = mv is 
mass x velocity, i.e. the momentum of the fluid). Note that the force F, momentum p and velocity v 
are all vectors, hence the boldface notation. The set of equations that describe F = d(mv)/dt 
applied to a fluid, including viscosity effects, is called the Navier-Stokes equations , shown here in 2 
dimensions, for an incompressible fluid (p = constant) that follows Newton’s law of viscosity (x = 
lidv/dy as described above): 



dv.. dv 
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dx 


dy 


dv dv y 
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(x momentum) 
(y momentum) 


dx 


+ 




(mass conservation) 


(Equation 54) 


Here v x and v v are the components of fluid velocity vector v in the x and y directions, respectively. 
Of course, at every point (x, y) in the flow, the velocity components v x and v may be different. The 
left-hand side of the first two equations is basically the d(mv)/dt = m(dv/dt) + v(dm/dt) terms, i.e. 
the rate of change of momentum of the fluid. In particular, the p(dv x /<3t) terms are just m(dv/dt) = 
mass x acceleration and the pv(<3v x /<9x) terms are just v(dm/dt), i.e. the increase or decrease in 
momentum within an infinitesimal volume due to an increase or decrease of mv within the volume. 
The right-hand side of the first two equations is the forces acting on the fluid due to pressure, 
gravity, and viscosity. The third equation is required to conserve the mass of fluid, i.e. it basically 
says that for a fluid of constant density, the rate of volume flow into an infinitesimal volume must 
equal rate of volume flow out of that volume. 
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These equations are very difficult to solve for all but the simplest situations. Certainly we will not 
try to do it in this course. For the purposes of this course, the key points to note about the Navier- 
Stokes equations are: 

1. The first two equations are just expressions of F = d(mv)/dt applied to a fluid 

2. There are two equations because momentum is a vector and thus there are x and y 
components; for a three-dimensional system another equation for the 2 component of 
momentum would be required 

3. The terms on the right hand side of the first two equations are just the forces F broken 
down into their components in the x and y directions 

4. The terms on the left hand side of the first two equations are just the rate of change of 
momentum d(mv)/dt broken down into their components in the x and y directions. 

5. There are three equations for the three unknowns v x , v and pressure P. (Unlike velocity, P 
is a scalar so it doesn’t have x and y components.) 

6. The equations are linear except for the v x (5v x /5x), v y (dvjd y), v x (3v /3x) and v (5v /5y) 
terms. This nonlinearity is very significant because 

a. It makes fluid mechanics difficult — the nonlinear terms are responsible for very 
complicated phenomena such as flow instabilities, turbulence and shockwaves. 

b. It makes fluid flow fundamentally different than linear systems. For example, if I 
have one solution to the Navier-Stokes equation, call it v 1 (x,y), Pi(x,y) and a second 
solution v 2 (x,y), P 2 (x,y), it is generally NOT the case that v 1 (x,y) + v 2 (x,y), P t (x,y) + 
P 2 (x,y) is also a solution. As an example of a linear system, consider traveling waves 
on a string. A rightward-traveling wave and a leftward-traveling wave can pass 
through each other without any change in the waves after the passage. However, a 
rightward-traveling flow structure (say, a spinning vortex) and a leftward-traveling 
flow structure will interact with each other in such a way that each will be 
permanently changed by the interaction. 

Laminar and turbulent flow 

When Re is low, which means that viscous effects are relatively important, the flow will be steady 
and smooth, which is called “laminar flow.” At higher Re, viscosity is not strong enough to 
suppress the instabilities (due to the nonlinear terms in the Navier-Stokes equations) and the flow 
becomes turbulent. While you have an intuitive feel of what turbulence is, a precise definition of what 
is or is not turbulent is not a simple matter and we will not attempt to define it rigorously here. 

The Reynolds number at the transition from laminar to turbulent flow depends on the type of flow, 
for example: 

• Flow in circular pipes: Re d = vd/v ~ 2,200 (v = mean velocity of flow in the pipe; d = 
inside diameter of pipe) 

• Flow along a flat plate: Re L = vL/v ~ 500,000 (v = velocity of flow far from the plate; L = 
distance from the “leading edge” of the plate.) 

Since value of Re at the transition from laminar to turbulent flow can be vary widely depending on 
the type of flow, the actual value of Re is not meaningful in comparing different flows. That is, one 
cannot say if 10,000 is a large or small value of Re until one also specifies the type of flow; it is a 
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relatively high value (well into the turbulent regime) for pipe flow, but a low value (well within the 
laminar regime) for flow along a flat plate. 

Lift, drag and fluid resistance 

Lift and drag coefficients 

Any object moving through a fluid will experience a force (F D ) in the direction opposing the motion. 
This force is called drag. Recall the definition of drag coefficient (page 13): 

F„ - i C D pv 2 A 

where C D is the drag coefficient, p is the fluid density, v the fluid velocity far from the object and A 
is the cross-section area of the object. 

An object moving through a fluid may also experience a force in the direction perpendicular to the 
direction of fluid motion. This force is called lift and is defined in a way similar to drag: 

1 9 

F l = —C L pv A (Equation 55) 

2 

where C L is the lift coefficient. While all objects moving through a fluid experience drag, only some 
will experience lift. The main goal of aircraft wing design is to maximize the lift to drag ratio, i.e. 
C L /C d . A glider may have a lift to drag ratio of 50, whereas commercial passenger aircraft wings are 
in the range 15 — 20 (which is about the same as an albatross). 


Flow around spheres and cylinders 

In the case of laminar flow at very low Re, and only in this case, the drag coefficient C D on a 
sphere is equal to 24/Re (for laminar flow at low Re only - got it???). Combining this result with 
the definition of drag coefficient and definition of Re, we obtain 

F drag = 37C(ivd (laminar flow around spheres) (Equation 56) 


If the sphere is moving due to gravity alone, the buoyant force F buoyant - (p fluid -p spheIe )gV - 

(Pfluid-Psphere)g(47t/3)r 3 = (p flu i d -p sphere )g(7t/6)d 3 . Note that the buoyant force does not depend on v, 
but the drag force does. Thus a dropped sphere will initially accelerate until its velocity is just that 
required for the drag force to equal the buoyant force, at which point there is no acceleration, and 
the velocity has reached a constant value called the terminal velocity. For the case of the sphere, 
equating F drag and F buoyant we obtain: 

v termina i = gd 2 (p fluid - Psphere)/18pt (laminar flow around spheres) (Equation 57) 


where the + sign is consistent with the fact that if p fluid > p sphere , 

v). 


the sphere moves upward (positive 
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The above terminal velocity is only valid for laminar flow. For turbulent flow, there is no simple 
analytical relationship between Re and C D , so one must resort to experiments or detailed (and 
difficult) computer simulations of the Navier-Stokes equations. Figure 20 shows a comparison of 
the actual C D vs. Re with that predicted by the low-Re laminar-flow model. It can be seen that the 
relation for laminar flow is reasonable up to about Re ~ 3 but at higher Re, the flow is NOT laminar 
and thus the laminar flow result C D = 24/Re does not apply. As one would expect, C D is higher 
with turbulent flow. Note also that at Re ~ 3 x 10 s there is a sudden decrease in C D . 

To provide a physical explanation of the drag coefficient plot, note that, apart from the small dip 
near Re ~ 3 x 10 5 , at high Re, C D is close to 1 and doesn’t change much with Re. This is because at 
high Re, the high momentum (relative to viscous effects) of fluid causes the flow behind the sphere 
separates and there is a region behind the sphere with v ~ 0. Thus, most of the dynamic pressure 
( =1 / 2 pv“) of the flow is lost, thus according to Bernoulli’s equation (which does not strictly apply 
because the flow is neither steady nor inviscid, but is still useful for estimation purposes) the 

pressure on the downstream side of the sphere is higher than that on the upstream side by Vipv . 
Hence, the net force on the sphere due to this separation-induced drag is F D = Vipv A, and thus the 

drag coefficient = F D /(V 2 pv 2 A) = (V 2 pv 2 A)/(V 2 pv 2 A) = 1. Thus, for any blunt object at high Re, C D 
is usually close to 1 (another example of “that’s easy to understand, why didn’t somebody just state 
that?”) At low Re (less than about 10 for the sphere), the drag is predominantly viscous and thus C D 
is higher. 
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A similar C D vs. Re plot for cylinders in cross-flow (i.e. with the flow in the direction perpendicular 
to the axis of the cylinder) is also shown in Figure 20, but for cylinders there is no simple analytical 
relationship analogous to the C D = 24/Re result for laminar flow over spheres. Note again the 
sudden decrease in C D at almost the same Re as for spheres. Note that a 3 cm golf ball hit at 70 m/s 
in air has a Reynolds number of (70 m/s)(0.03 m)/(1.5 x 10° m 2 /s) = 1.4 x 10 s . Dimpling the golf 
ball decreases the transition Re somewhat, and thus enables a lower C D (actually the dimpling also 
increases the lift due to the backspin on the ball, but that’s beyond our scope.) 

Example 

(a) Howard and Samantha go skydiving. Howard weighs 175 lbf with all his gear, and has a 
cross-sectional area when free falling of 8 ft. If Howard’s terminal velocity is 150 miles per 
hour, what is his drag coefficient? 

1 2 

Drag force F D = —C D pv~A 



2 F 

D _ 


o i 4.448A 

2x175/0/- 

wf 

pv 2 4 A 

so 

oc 

H 

. 

i-H 

^ 150m/5280/1 hr m 


m 1 

\ hr mi 3600 sec 3.281 ft t 


2 




3.281/1 


0.395 


(b) Samantha has an unusual skydiving style. She free-falls lying perfectly straight horizontally, 
and her shape can be treated as roughly that of a cylinder 5 ft long and 2 ft in diameter. She 
weighs 125 lbf with all her gear. Assuming that her drag coefficient can be modeled as that 
of a circular cylinder, what is her terminal velocity? To do this problem you will have to 

1) Guess a terminal velocity 

2) Compute her Reynolds number 

3) Look up her drag coefficient in Figure 20. 

4) Compute her drag force 

5) Does her drag force equal her weight? If not, adjust your guess of terminal velocity 
and go back to step 2. 


1) “Guess” v — 113 mph — 50.51 m/s 



Reynolds # Re - * . 

v 1.5x10 ~ 5 m 2 /s 


2.053xl0 6 


3) From Figure 20, C D ~ 0.4 

4) Compute drag force: 
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F d = -C D pv 2 A = —(0.4) (1.18%/ ot 3 )(50.51ot / sf (2 ft)(5 ft)(m / 3.28 1 ft) 2 
2 2 

= 559.47V = I25lbf 

5) Does her drag force equal her weight? 

Yes, 125 lbf = 125 lbf so mission accomplished. 

Flow through pipes 

For flow through pipes, drag coefficient is not used because what we’re really interested in is not the 

drag force but rather the pressure drop (AP), and thus a slightly different quantity called the friction 
factor (f) is used to quantify the effect of viscosity on the flow in the pipe: 

A P 
pv 2 L 
2 d 

where v is the average velocity of the fluid flowing thought the pipe, p the fluid density, L is the 
length of the pipe and d its diameter. For laminar flow only in pipes, f = 64/Re d where Re d = 

pvd/ju is the Reynolds number based on pipe diameter d, not pipe length L!, thus 

AP = (64/Re d )(pv 2 /2)(L/d) = (64g/pvd)(pv 2 /2)(L/d) = 32juvL/d 2 (Equation 59). 

Sometimes it’s more convenient to deal with volume flow rate (Q) rather than velocity (v). Q is the 
velocity multiplied by the cross-section area of the pipe, thus Q = v7id 2 /4. Thus we can write one 
last relation: 

AP = (128/7i)juQL/d 4 (laminar flow only!) (Equation 60). 

Note the significance of this result: if you double the flow rate Q or the length of the pipe L, the 
pressure drop doubles (makes sense.) Also, for a given flow rate Q, if you double the diameter of 
the tube, the pressure drop decreases by a factor of 16! So use a bit bigger pipe in your plumbing 
design! 

The results leading to the last 2 equations assumed f = 64/Re d and thus are valid only for laminar 
flow. For turbulent flow, the friction factor depends not only on Re d but also the roughness of the 

pipe wall, which is characterized by a roughness factor — 8/d, where 8 is a measure of the roughness (i.e. 
height of the bumps on the wall) and d is (as always) the pipe diameter. The combined effects of 
roughness and Re d are presented in terms of the Moody chart (Figure 21). 



(Equation 58), 
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Figure 21. “Moody Chart” showing the effect of Re d (that is, Reynolds number based on pipe 
diameter not length) and surface roughness s/d on the friction factor (f). 


Note that laminar flow prevails up to Re d = 2,200 (this value is essentially independent of the pipe 
roughness factor), then for higher Re d , C D increases suddenly but in a way that depends on the pipe 
roughness — as one would expect, rougher pipes have higher C D . It’s remarkable (to me, anyway) 
that at high Re d a tiny amount of roughness has a huge effect on f. For example, at Re d = 10 8 , f 
increases by a factor of 3 as one changes from a perfectly smooth pipe (e/d = 0) to e/d = 0.001. In 
other words, a roughness of one part in 1000 increases the pressure drop by a factor of 3. Size does 
matter! 


Alternatively, if you don’t like using the Moody diagram, the following empirical formula for 
turbulent flow can be used (for laminar flow, use f = 64/Re d as mentioned above): 





(turbulent flow) 


(Equation 61) 


but note that this formula has f on both sides of the equation, and you can’t simplify it any further, 
so for a given Re d and e/d, you have to guess a value of f and see if the right and left hand sides of 
the equations are equal, and adjust your guess of f until the two sides are equal (this is called a 
transcendental equation , one that cannot be solved in closed form). 
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Example: A 50 foot long garden hose has an inside diameter of 5/8” and a roughness (s) of 1/32”. 
Water flows through the pipe at a velocity of 3 ft/s. 

a) What is the flow rate in gallons per minute? 

Flow rate = velocity x cross-sectional area 

= 3 ft/sec x (tc/ 4) (0.625 inch) 2 x (ft / 12 inch) 2 x (7.48 gallon / ft/) x (60 sec / min) 

= 2.87 gallon / min 

b) What is the pressure drop in lbf/in 2 ? 


Re d = vd/v 

= [(3 ft/sec)(m/3.281 ft)] [(0.625 inch)(ft/12 inch)(m/3.281 ft)] / 1.0 x 10' 6 m 2 /sec 


— 14,515 > 2200 turbulent 

8/d = (1/32)/(5/8) = (1/32)/(20/32) = 1/20 = 0.05 


Use 


1 


77 


= -2 log 


Eld 2.51 

+ 


3.7 Re„J/ 


and “guess” f — 0.0731: 



1 


Vo.0731 


= 3.699;-2 log 


eld 2.51 


3.7 Re, 7/ 


= -2 log 


0.05 


2.51 


3.7 14,515V0.0731 


3.698 


the equation is satisfied, so f — 0.0731. Then 

r A P ^ A p _ rPv 2 L _ 0Q731 (1000fcg/m 3 )[(3/?/sec)(m/3.28iyi)7 (50 ft)(12inch I ft) 

1 pv 2 zA 1 2 d 2 0.625 inch 

2 d 

= 29,3351V/m 2 = 4.26/6/ /in 2 


Do you think I did these calculations by hand? No way! I used an Excel sheet (double click to 
open). The cells shaded in blue are the things you change, and the other cells are calculated values, 
except for the “Friction factor (guess)” cell, which you have to adjust until the left-hand side (LHS) 
and the right-hand side (RHS) of the equation for the friction factor are equal, and thus the “fraction 
error” goes to zero. (You can also use Excel’s “goal seek” feature to do this adjustment 
automatically.) 
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Pipe dia (in) 

Pipe len (ft) 

Roughness 

(in) 

Viscosity 

(m A 2/s) 

rho 

(kg/m A 3) 

0.625 

50 

0.03125 

1.00E-06 

1000 

Pipe dia (m) 

Pipe len (m) 

epsilon 



0.01587423 

15.23925632 

0.05 



Velocity 

(ft/s) 

Velocity (m/s) 

Re 



3 

0.914355379 

14515 



Friction 

factor 

(guess) 

LHS 

RHS 

fraction 

error 


0.07309956 

3.698645 

3.698296 

9.4274E-05 


Delta P 
(N/m A 2) 

Pressure drop 
(lbf/in A 2) 




29335.00 

4.26 





Compressible flow 


All of the above discussion of fluid mechanics relates to cases with constant density (p), which is 
certainly reasonable for liquids (e.g. water) under most conditions and even air if the velocity (v) is 
“small enough”. How small is small enough? We have to compare U to something else that also 
has units of velocity. That “something else” turns out to be the speed of sound (c). The ratio of 
these is the Mach number (M), i.e. 


M — v/c Equation 62. 

For an ideal gas, the sound speed c is given by the formula 

c = (yRT) 1/2 


where 


• y is the specific heat ratio of the gas (—1.4 for air at ambient temperature, but may be as low as 
1 for gas molecule with many atoms, and as high as 5/3 for a monatomic gas like helium) 

• R is the gas constant for the specific gas of interest = 3i/M, where 9^ is the universal gas 
constant = 8.314 J/mole K and M is the molecular mass of the gas (in kg/mole, = 0.02897 
kg/mole for air, thus R = 287 J/kgK for air.) 

• T is the gas temperature (in K of course) 

How does Mach number affect density (p), temperature (T) and/or pressure (P)? That depends on 
the process the gas experiences as it accelerates or decelerates. A detailed discussion of compressible 
gas dynamics is way beyond the scope of this course, but T11 give you the results for the simplest 
case of one-dimensional steady flow of an ideal gas in a duct of changing area A with constant 
specific heats (C P and C v ) between locations 1 and 2 assuming no heat transfer, no friction and no 
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shock waves (we call this special case “isentropic flow,” meaning no change in the entropy of the 
gas) as well as no potential energy (elevation) change : 






= rJi + 


y - 1 

2 




Equation 63 


These equations are plotted in Figure 22 below. 


While a lot of simplifying assumptions were made (note all the underlined words above), this 
“isentropic flow” is still useful as the simplest model of flow in nozzles of jet and rocket engines, as 
well as intakes in jet engines. Note that as Mach number increases (for example, during expansion 
in a nozzle), pressure, density and temperature all decrease. However, to obtain transition from 
subsonic (M < 1) to supersonic (M > 1) flow, the area must pass through a minimum, i.e. a throat, 
which occurs at M = 1. Thus, rocket nozzles must have an hourglass shape in order to accelerate 
the exhaust to high Mach numbers and therefore produce the maximum possible thrust. Often the 
areas are referenced to the minimum area at the throat (A) where M = 1, in which case 



Equation 64 
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Figure 22. Plot of pressure, density, temperature and area as a function of Mach 

number for one-dimensional, isentropic flow of an idea gas with y= 1.4. The “t” 
subscript indicates the pressure, density or temperature when M = 0. (Double-click 
to open Excel spreadsheet). 


One important but often overlooked point about the above equation: the Mach number M is the 
local (at station 1 or 2) flow velocity U divided by the local sound speed (which depends on the local 
temperature.) So you can’t divide the local velocity by the sound speed at ambient temperature to 
get the Mach number! That is, M t = v 1 /(yRT 1 ) 1/2 and M 2 = v 2 /(yRT 2 ) 1/2 but you can’t say M 2 = 
v 2 /(YRT0 1/2 ! 


How to scrutinize this result? The units are clearly ok since y and M are dimensionless. Also, if A 1 
= A 2 then Mj = M 2 , that is, nothing changes. But here’s a great function test: in the limit of small M 
(small compressibility effects), the results should reduce to Bernoulli’s equation. The second 
relation involves pressure (P) and velocity (v), so looks a lot like Bernoulli. Recall the binomial 
expansion theorem which says that for m << 1,(1 + m) n ~ 1 + mn, thus 
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This is not quite Bernoulli’s equation, which requires p = p { = p 2 = constant. But from the first 
part of Equation 62, in the limit of small M 1 and M 2 , p t = p 2 . Notice that the exponent on the 

density terms, 1 /(y-1), is smaller than that on the pressure terms, y/(Y"l)> an d thus density can be 
assumed constant even when pressure is not. A formal derivation requires carrying out higher order 
terms (i.e. M 4 terms) in the binomial expansion so I’ll skip that... 

Example 

The (now decommissioned) SR-71 aircraft flew at Mach 3 at an altitude of 80,000 feet. Assuming 
isentropic flow, what is the temperature and pressure on the leading edges of the wings where the 
flow (in the frame of reference of the aircraft) has decelerated from M 1 = 3 to M 2 = 0? 

From http://www.digitaldutch.com/atmoscalc/calculator.htm : at an altitude of 80,000 ft, the 
standard atmospheric conditions are = 0.0273 atm (that is, 0.0273 sea-level atmospheres!) and Tj 
= 221 K. 
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T 2 =619A: = 654F 


So the dynamic pressure loading on the wings (1 atm, or 14.7 lbf/in 9 ) and the temperature (654°F) 
are really high, and things get worse in a hurry as increases - notice the terms appear 

everywhere, and there is an additional y/ (y-1) ~ 3.5 exponent on the pressure equation — so dynamic 
pressure increases like M in high-speed flight! 


If you’re dying to know more about compressible flow, check out 
http://ronney.usc.edu/AME436Sl 3/Lecture 11 . 

Example 


Can the air flow over a car be considered an incompressible flow (p ~ constant)? In decelerating the 
air flow from say 75 mph to 0 mph, how much does the density change assuming isentropic flow? 
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1 + 
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2 




= 1.0047 


In other words, the density changes by less that 0.5% in decelerating from 75 mph to 0 mph, so for 
all but the most stringent accuracy requirements, the density change can be neglected. In other 
words, for most practical purposes the air flow over a vehicle at 75 mph can be treated as 
incompressible , i.e. in the same way as if the flow were water at the same Reynolds number. 
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Chapter 8. Thermal and Energy Systems 


“Passion is energy. Feel the power that comes from focusing on what excitesyou .” — Oprah Winfrey 

Main courses in AME curriculum on this topic: AME 310 (Thermodynamics); AME 331 
(Heat Transfer). 

Conservation of energy - First Law of Thermodynamics 

Statement of the First Law 

The First Law of Thermodynamics states that “you can’t win,” meaning that energy is conserved, i.e., 
the energy contained in an isolated system (one that does not exchange energy with its surroundings) 
cannot change. Of course, energy can be converted from one form to another, which is the whole 
point of energy engineering. 

Quantitatively, the 1st Law of Thermodynamics for a control mass, i.e. a fixed mass of material (but 
generally changing volume, for example the gas in a piston/cylinder) can be stated as follows: 

AE = AQ - AW Equation 65 

where 


E = total energy contained by the mass - a property of the mass (in Joules, BTUs, etc.) 

Q = heat transfer to the mass (in Joules, BTUs, etc.) 

W = work transfer to or from the mass (see below) (in Joules, BTUs, etc.) 

A = Change from state 1 to state 2; thus the 1 st law could also be written 

1-2 — Wi-2 

i.e., the change in energy contained by the mass is equal to the heat transferred to the mass minus 
the work transferred out of the mass. Work transfer is generally defined as positive if out of the 
control mass, in which case - sign applies, i.e. AE = AQ - AW; if work is defined as positive into 

system then AE = AQ + AW. A process in which no heat transfer occurs is called an adiabatic 
process. 

Note that there’s nothing profound about the above equation, it’s just “energy bookkeeping.” It 
merely states that the change in the energy E contained by a substance is equal to the energy transfer 
to the substance (via heat transfer Q) minus the energy transfer from the substance (via work 
transfer W). 

In the above equation, 1 is the initial state or condition of the system (temperature, pressure, volume, 
etc.) and 2 is the final state. 1—>2 is the process or series of states leading from the initial state 1 to the 
final state 2. 


E 2 -E 1 = Q 
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What is the difference between heat and work? Why do we need to consider them separately? 

1. Heat transfer is disorganized energy transfer on the microscopic (molecular or atomic) 
scale and has entropy transfer associated with it. (What is entropy? We’ll talk about this in 
the context of the 2 nd Law of Thermodynamics, but basically it’s a measure of the level of 
disorganization of the system.) 

2. Work transfer is organized energy transfer which may be at either the microscopic scale or 
macroscopic scale and has no entropy transfer associated with it. 

The total energy of the substance (E) consists of 

• Macroscopic kinetic energy (KE = Vzmv 2 ) (m = mass, U= velocity) 

• Macroscopic potential energy (PE = mgz) (g = acceleration of gravity, z = elevation) 

• Microscopic internal energy (U) (which consists of both kinetic (thermal) and potential 
(chemical bonding) energy, but we lump them together since we can’t see it them separately, 
only their effect at macroscopic scales.) Generally this is written not as E but as mu, where 
m = mass and u = internal energy per unit mass (units Joules/kg, BTU/lbm, etc.). 

Thus, the total energy contained by a piece of material is given by 
E = KE + PE + U = Vzmv 2 + mgz + mu 

(See the “energy family tree” ( Figure 23)). 



Figure 23. Energy family tree. 


There are several ways to transfer heat to/from a system, that is, by conduction, convection and/or 
radiation as will be discussed later. Another way is to have a chemical reaction, for example 
combustion, occur within the mass. Stricdy speaking, this is not heat transfer, it is a change in the 
potential energy part of internal energy of the mass (usually changed into the kinetic part of the 
internal energy, i.e. the substance gets hotter). But chemical energy release due to combustion is 
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often modeled as heat transfer from an external source. How much? That depends on 2 things: the 
mass of fuel being burned (m) and the heating value of the fuel, denoted as Q R : 

AQ (due to combustion) = mQ R Equation 66 

The units of Q R are Joules/kg, but be careful — this is per kg of fuel, not per kg of fuel+air mixture! 
The most fuel I can add without wasting fuel is when I have just enough oxygen (from the air) to 
burn all of the carbons to make C0 2 and all the hydrogens to make H 2 0. This is called the 
stoichiometric mixture, and the ratio of (fuel mass)/(fuel mass + air mass) at stoichiometric is about 
0.064 for typical hydrocarbons — which means the mixture is mostly air. If I have a higher mass 
ratio, then I’ll be adding fuel that I can’t burn because I don’t have enough oxygen, which is a waste 
of fuel and generates pollutants such as CO (carbon monoxide) (poisonous) and unburned 
hydrocarbons that helps create ozone (0 3 ) (bad stuff!) in the atmosphere. 

Some typical values of Q R (in J/kg) are given in Table 4. 


Fuel 

Heating value, Q R 

a/kg) 

(fuel mass)/(fuel mass + air mass) 
at stoichiometric 

Gasoline 

44 x 10 6 

0.0642 

Methane 

50 x 10 6 

0.0550 

Methanol 

20 x 10 6 

0.104 

Ethanol 

27 x 10 6 

0.0915 

Coal 

34 x 10 6 

0.0802 

Paper 

17 x 10 6 

0.122 

Fruit Loops 

16 x 10 6 

Probably about the same as paper 

Hydrogen 

120 x10 6 

0.0283 

U 235 fission 

82,000,000 x 10 6 

n/a 


Table 4. Heating values of some common fuels 

By comparison, the energy content of the lithium-ion batteries used in your cell phone and laptop 
computer (that is, by discharging the battery to extract its electrical energy, not burning the battery!) 
is about 0.8 x 10 6 J/kg — more than 50 times less than hydrocarbon fuels. This is why most of 
us don’t drive a battery-powered car, and none of us fly in battery-powered aircraft! Also, starches 
and sugars are all about 110 cal/ounce (just look at the nutritional information on the side of the 
box of cereal or any other dry food), which translates to 16 x 10 6 J/kg. Note also that nuclear 
energy sources such as uranium-235 fission have millions of times more energy per unit mass than 
fuels, which explains their value for bombs, submarine propulsion, etc. 

Describing a thermodynamic system 

In order to characterize thermodynamic systems, we need to describe the behavior of the material as 
its temperature, pressure, volume, etc. changes. Hence the following terminology has been 
developed: 
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• Property — a quantitative description of a piece of matter (e.g. pressure, temperature, density, 
volume, mass, internal energy, enthalpy, entropy, as well as the composition) 

• State — a list of the properties of a piece of matter (e.g. liquid water at 25°C and 1 atm) 

• Process — a sequence of states (e.g. heat transfer to water at a constant pressure of 1 atm, with 
temperature increasing from 25°C to 35°C) 

• Cycle — a sequence of processes that returns to the original state (e.g. compress liquid water at 
25 °C from 1 atm to 100 atm, transfer heat at constant pressure until the temperature is 400°C, 
expand the water (actually steam) until the pressure is back to 1 atm, then transfer heat from the 
water (actually steam) until the temperature is 25 °C again) 

Conservation of energy for a control mass or control volume 

With the aforementioned forms of energy only (internal, kinetic, potential) along with the two types 

of energy transfers (heat and work) the first law for a control mass can be written as 

(mu, + 1 / 2 mv 2 + mgz 2 ) " (mu t + 1 / 2 mv 1 + mgz/) = Qi_ 2 —W 12 or 

(u 2 + Vzv 2 + gz 2 ) - (u t + 1 / / 2V 1 2 + gz t ) = (Q 1-2 — WjJ/m Equation 67 


where the subscript “1” indicates some initial state and “2 
occurred. 


?? 


a final state after some process has 


This form of the first law is useful for a “control mass” i.e. a fixed mass of material than may undergo 
volume changes. But what if I have a system in which mass flows in and out, for example a water 
turbine or a jet engine? In that case it’s usually more convenient to work with a “control volume ,” i.e. a 
fixed volume in space whose size and shape does not change but does allow for the possibility of 
mass flow in and out. In this case, it’s generally more convenient to write the first law as a rate 
equation: 


dE 

dt 



-W + m in 




\ 


+ 8Zm 


-m 


OUT 



Equation 68 


where 


E — energy contained by the control volume 

m in fa out ~ tnass flows into / out of control mass (kg / s) 

Q,W — rates of heat transfer in and work transfer out (Watts) 

h in , h out = enthalpy of material at inlet/outlet per unit mass (Joules/kg) = u + P/p, where u = 
internal energy per unit mass as above, P = pressure and p = density. 

Note that in this case, instead of using the subscripts “1” and “2” to denote the states before and 
after the process as was done with the control mass form of the first law, in the control volume 
form the subscripts “in” and “out” are used. This is because for the control volume one must 
distinguish between an inlet, and thus causes energy to be added to the control volume, from an 
outlet, which causes energy to leave the control volume. The difference is only a sign (+ or -), but 
the sign is important! On the other hand, for the control mass form of the first law, the subscripts 


79 



“1” and “2” are interchangeable, that is, 1 could be at the beginning of the process and 2 and the 
end or vice versa, and the conservation of energy is still enforced. When the Second Law of 
Thermodynamics is considered, however, there is a very definite requirement as to which state, 1 or 
2, happened first. Only for a very special type of process, called a reversible process , could either 1 or 2 
be the initial state. 

For many materials over not too large a temperature range, 

h 2 — hj ~ C P (T 2 — Tp) Equation 69, 

where C P is the specific heat or. heat capacity at constant pressure of the material (units J/kg°C) and T is the 
temperature. Also, the internal energy per unit mass (u) can be represented in a similar way: 

u 2 — u i ~ C v (T 2 — T t ) Equation 70, 

where C v is the specific heat or heat capacity at constant volume of the material (units again J/kg°C). The 
ratio of C P to C v is called the specific heat ratio (y) that we’ve already used: 

y = C P / C v Equation 71. 

For liquids and solids, y is very nearly 1 so this distinction between C P and C v is of no practical 
importance. But it is hugely important for gases; even though for gases 1 < y < 5/3, many of the 
formulas (e.g. at the end of Chapter 7) involve a factor of (y -1). Important point: the terms 
heat” or “heat capacity” at “constant pressure” or “constant volume” are 
misleading for two reasons: 

1. Substances do not contain heat, they contain energy; heat is a mode of energy transfer 
between a substance and its surroundings 

2. C P and C v can be and often are used even in processes for which pressure or volume is not 
constant 

Despite these atrocities of terminology, we are forced to use the terms because they are engrained in 
the vernacular of science and engineering. “Specific enthalpy” and “Specific internal energy” would 
be much better terms for C p and C v , respectively (in which case we would probably call them C h and 
Q instead.) 

Note that the control volume form of the First Law looks a lot like the control mass form; about the 
only thing one couldn’t figure out by inspection is the substitution of the enthalpy h for internal 

energy i. Note also that since h = u + P/p, thus the difference between h and u is just P/p. This is 
due to something called the “flow work” required to push the material into the control volume and 
the work obtained when extracting it from the exit of the control volume. This flow work term 
doesn’t apply to the control mass since (by definition) there is no mass entering or leaving the 
control mass! 

A very useful and important special case of the control volume form of the First Law is the steady 
state , steady flow case, where all properties (E, h m , h out u ln , u out , z in , z out ) and all fluxes (m in ,m oun Q,W) 


terribly 
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are constant (not changing over time) and moreover m in = fn out (otherwise the mass contained 
within the control volume would change over time). In this case the First Law is written as 


0 = Q-W + m 


(K- h ou,) + 




2 + (sz in -gz out ) 


Equation 72 


Processes 

As previously mentioned, a “process” is a sequence of states. Normally in simple thermodynamic 
analyses, one assumes that one of the properties of the substance (temperature, pressure, volume, 
entropy, internal energy, etc.) is constant during a given process. It is beyond our scope to identify 
which property is most nearly constant during a given type of process except to mention a few 
specific cases. For heat addition due to combustion in a piston-type engine, the process is nearly 
constant volume. For steady-flow heat addition due to combustion in a gas turbine, or heat transfer 
to water/steam in a boiler, the process is nearly constant pressure. Compression of a substance is 
usually idealized as being “reversible,” meaning that the compression can be reversed (i.e. the 
substance can be expanded) until the substance returns to its original state and the same amount of 
work transferred into the system can be transferred out during expansion. This is similar to an ideal 
spring. Moreover, compression processes are usually idealized as being adiabatic (without heat 
transfer). It can be shown that a reversible and adiabatic process results in no change in the entropy 
(discussed later) of the substance. Furthermore, it can be shown that isentropic compression or 

expansion of an ideal gas with constant specific heat ratio (y) follows the relation 

PjVj 7 = P 2 V y 2 (ideal gas, constant y, reversible adiabatic process) 

P = pressure, V = volume 


Y 

This relationship (PV = constant) is called the isentropic compression law (but note the restrictions on 
when it applies! Ideal gas, constant specific heat ratio, reversible and adiabatic process!) 

Examples of energy analysis using the 1 st Law 

Example #1 — gas compression 

For isentropic compression of 480 cm' 1 (= 480 x 10° m~) of air in a cylinder of car engine initially at 
300K and 1 atm (= 101325 Pa) by a volume ratio of 8, neglecting kinetic and potential energy in the 
gas, 


a) What is the pressure and temperature of the air after compression? 

b) What is the work required? 

c) If there are 8 cylinders and the engine rotates at 3000 RPM, what power is required to do 
this compression? 

a) For an ideal gas with constant y undergoing an isentropic process , 

PjVj 7 = P 2 V 2 7 where P = pressure, V = volume 
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Thus P 2 /P\ = V\ Y IV 2 or P 2 = P l {V l IV 2 ) — (1 atm)(8) 1-4 — 18.4 atm. 

Since this is an ideal gas we can also say PjVJ = mRTj and P 2 V 2 = mRT 2 ; combining these 
relations with Py Y = P 2 VJ we obtain T 2 =T l (V l IV 2 ) 

V, = 480 cm 3 , V 2 = 480/8 = 60 cm 3 , T t = 300K => T 2 = 300K (480/60) 0 * 4 = 689.2 K 

Note that the pressure after compression (P 2 ) seems to fail the function test — how can the 
pressure ratio be more than the volume ratio, that is, how can the post-compression pressure 
be more than 8 atm? It’s because we put work into the gas, as evidenced by the temperature 
rise. If the temperature were constant during the compression, then from the ideal gas law 
PV = mRT with m, R and T all constant, we would have P 2 /P\ = Vj/V 2 = 8 (or less, if there 
were leaks in the cylinder, which is the whole reason for doing a compression test — to check 
for leaks.) This confused me as a fledgling auto mechanic in high school; even then I was 
doing function tests, and I couldn’t understand why when I did a compression test on my 
engine, the pressure ratio was higher than the volume ratio. I was thinking the process was 
isothermal, not isentropic (I didn’t know that in high school...) 

b) Treat the mass of gas in the cylinder as a control mass 

E 2 - E t = Q U2 - W, 2 

In an isentropic ( reversible, adiabatic ) process there is no heat transfer (which is what adiabatic 
means), only work transfer, so Q 12 = 0. Also, KE = PE = 0, so 

E = U + KE + PE = U, thus 

W w = E, - E 2 = U, - U 2 = m(uj - u 2 ) = mC y (T 1 - T>). 

From the ideal gas law, m = P/V^/RTp for air at 300K, R = 287 J/kg°C = 287 J/kgK (see 
http://www.efunda.com/materials/common matl/show gas.cfm?MariName=AirOC) for 
example), thus 

m = P^/RT 1 ! = (101325 Pa)( 480 x 10 6 m 3 )/(287 J/kgK)(300K) = 5.65 x 10' 4 kg 
For air at 300K, C v ~ 720 Jkg°C = 720 J/kgK (same website) thus 

W t _ 2 = (5.65 x 10' 4 kg)(720 J/kgK)(300K - 689.2K) = -158.3 J 
(Function test: work is negative because it’s work going into the gas) 

c) All automotive engines are 4-stroke engines and have only one compression stroke for every 
2 revolutions of the engine, so there are only 1500 compression strokes per minute per cylinder, 
thus 

Power = work/time = (-158.3J/comp)(1500 comp/min cyl) (8 cyl) (min/60 sec) 

= -3.17 x 10 4 J/s = -3.17 x 10 4 Watts x (1 hp / 746 Watts) = -42.4 horsepower 

(Performance test: this sounds like a lot of power, but keep in mind (a) this is a fairly large 
engine, 480 cm/ cyl x 8 cyl = 3840 cm' 1 = 3.84 liter; (b) the engine rotation rate is fairly high, 
3000 rev/min, whereas your typical highway cruise is closer to 2000 rev/min and (c) this 
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assumes the air is coming in at 1 atm, which means wide open throttle, i.e. “pedal to the 
metal.” At this condition the engine would produce well over 100 net horsepower.) 


Wouldn’t it be better not to compress the air, and get 42 more horsepower? No way! If you 
don’t compress, you don’t get any power at all. You don’t have an engine any more, just a 
complicated heater. As it turns out, and you’ll learn if you take AME 436, that you get more 
work out of expanding the hot gas than the work input required to compress the cold gas by 
the same volume ratio. This is the only reason that internal combustion engines work. In 
fact, the higher the volume compression ratio, the more work you get out for a given 
amount of heat input and thus the higher efficiency you get. This is discussed further below. 

Example #2 — potential energy 


The Upper Fall of Yosemite Falls in Yosemite National Park is a sheer plunge of 435 meters. 
Assuming no air drag (yeah, right), no heat transfer to/from the air or rocks (yeah, right) and no 
work extracted from the falling water (hey, this is a National Park, no hydroelectric plants allowed!) 

a) What is the velocity of the water just before it hits the rocks at the bottom of the upper falls? 
Assume that the water is nearly at rest at the top of the falls. 

b) After churning around in the rocks (with no further change in elevation) until the velocity is 
very small compared to that just before hitting the rocks, how much warmer is the water? 

a) Draw a control volume where in = top of falls and out = bottom of falls, just above the 
rocks. Then from the steady state , steady flow form of the First Law, with no heat transfer, no 
work transfer, and (until after the water hits the rocks) no change in enthalpy: 




-gZ in with V; 



9.81 m/s, (z out -z in ) 


-435 m 



92.4 m/s 



Draw a control volume where in — bottom of falls, just above the rocks and out — water 
downstream after churning around in the rocks until the velocity is very small. 






, with v in 


92.4 m/s, v out 


and (for water) C P — 4184J/kg°C 



C P (T in -T 0 J, 


T out - T in = [(0 m/s) 2 - (92.4 m/s) 2 ]/2 /4184J/kg°C = 1.02°C = 1.84°F 


Example #3 — kinetic energy 


A jet engine on an aircraft flying at 500 mi/hr has an inlet air mass flow of 10 kg/s, an inlet air 
temperature of 250K, a fuel mass flow of 0.3 kg/s, and an exhaust temperature of 900K. All flows 
are steady. What is the velocity of the jet exhaust? Assume C P = 1400 J/kg°C for fuel, air and 
exhaust. The heating value of jet fuel is the same as gasoline. Neglect elevation change, and neglect 
any work extracted from the engine (e.g. to drive an electrical generator.) 


83 



In this case we have two inlets, one for fuel (v lnfuel = 0) and one for air (v lnair = 500 mi/hr), but only 
one outlet (for the jet exhaust) (v out = ???). The heat input is the mass flow rate of fuel multiplied by 
the heating value of the fuel. Thus we can apply the control-volume form of the first law as follows: 
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= Q-W + m. 
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By using the 2 nd Law (which we haven’t covered yet) one can determine the exhaust temperature 
without having to specify it as we have done here. Note that the exhaust velocity is much higher 
than the flight velocity, which is required if we want any net thrust! In practical jet engines, however, 
only a portion of the fuel is burned and most of the air ingested into the engine goes through a giant 
fan without being burned. The power needed to drive this fan comes from taking work out of the 
combusted stream through several stages of turbines. This type of engine is called a turbofan and is 
much more efficient that the turbojet we just analyzed. The reasons that the turbofan is more 
efficient is discussed in great detail in AME 436. 


Second Law of thermodynamics 

The First Law of Thermodynamics told us that energy is conserved, i.e., the energy contained in an 
isolated system (one that does not exchange energy with its surroundings) cannot change. But this 
isn’t the whole story, because it does not place any restrictions on the direction of a process. For 
example, one can readily fill a (constant-volume) combustion chamber with a mixture of methane 
and air at 300K, ignite the mixture with a spark, and observe a flame burn the mixture to form 
carbon dioxide, water and nitrogen at 2000K. Clearly this does not violate the First Law. But when 
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was the last time you saw carbon dioxide, water and nitrogen at 2000K in a chamber spontaneously 
cool off to form methane and air at 300K? Clearly this does not violate the First Law either, since 
energy is conserved in either the forward or reverse direction, but you have never seen the reverse 
process and you never will. 

So clearly we need a Second Law of Thermodynamics that places restrictions on the direction of 
processes. The Second Law invokes a property of substances called entropy , which is the measure of 
the “disorganization” or “randomness” of a substance. The hotter or less dense a substance is, the 
less information we have about where the individual molecules are, and thus the higher its entropy 
will be. The Second Law can be stated simply as 

The entropy of an isolated system always increases or remains the same 

meaning that the entropy never decreases. The methane — air mixture at 300K has a lower entropy 
than the carbon dioxide, water and nitrogen mixture at 2000K, so only the usual combustion process 
is physically possible, never the reverse. (Of course I could take that carbon dioxide, water and 
nitrogen at 2000K, cool it off to 300K, break the molecules apart, rearrange them to form methane 
and air, but to do this I would need to increase the entropy of the surroundings by more than the 
entropy change of combustion, so there would be a net increase in the entropy of the universe.) 

A detailed discussion of entropy and the Second Law is far beyond the scope of this course, so only 
two important consequences of the Second Law will be mentioned here. The first such consequence is 
that 


It is impossible to create a device that has no effect other than the transfer of heat 
from a lower temperature to a higher temperature. 

If this were not true, then it would be possible for an object initially at uniform temperature to 
spontaneously become hotter on one side and colder on the other — which is obviously a more 
organized (lower entropy) state than the original, uniform-temperature object. This statement is 
sometimes stated as “heat is always transferred from hot to cold, never the reverse” — which is only 
a requirement if there is no other effect. Obviously a refrigerator transfers heat from a lower 
temperature (your food and drink) to a higher temperature (the air inside your kitchen) but it has 
other effects too — namely there is a work input to the process. The second consequence of the 2 nd 
law that will concern us is the limitations on the efficiency of heat engines and refrigerators, 
discussed below. 


Engines cycles and efficiency 

The First Law said that “you can’t win”, meaning that one cannot, for example, transfer 1 Joule of 
energy into a device in the form of heat transfer and get more than 1 Joule of work transfer out of 
said device (despite the fact that every day the news media reports that someone somewhere in the 
world has done exactly that.) The Second Law is an even more insidious and depressing because it 
says in effect that “you can’t break even, ” meaning that for 1 Joule of energy into a device in the 
form of heat transfer, one cannot even get as much as 1 Joule of work transfer out of said device. 
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How much work output can one obtain if not 100% of the heat transfer? To determine this, first 
we’ll state that “it can be shown” that the entropy production in any heat transfer process is given by 
Q/T, where Q is the amount of heat transfer and T is the temperature (must be absolute 
temperature, K not °C!) at which the heat transfer occurs. It is important to recall that there is no 
entropy transfer associated with work transfer, so there are no W/T terms to consider. The 
same thing can be said about a change in the macroscopic kinetic energy and/or potential energy of 
the substance; there is no entropy change associated with macroscopic kinetic or potential energy 
changes. So if we transfer Q H units of transfer heat into an engine (or any device) at temperature T H 
and transfer Q r units of transfer heat out of that system at temperature T L , the net entropy 
production must be positive, i.e. 

Entropy production = Q H /T H - Q L /T L > 0 => T L /T H > Q L /Q H 

Where the negative sign in front of Q L appears because it is out of the system. Thus for the best 
possible heat engine, Q H /T H — Q L /T L = 0 or 

t l /t h = q l /q h 

The efficiency (r|) of the engine is the ratio of work output (W) to heat input (Q H ) (see page 15), and 
by the first law W = Q H - Q L , thus 

*1 = (Qh - Qi)/Qh = 1 - Ql/ Qh (for any engine) 
and for the best possible engine with zero net entropy production 

r\ = 1 - T l /T h (for the best possible engine) Equation 73 

This best possible device is a theoretical abstraction called a Carnot cycle engine . While it’s not obvious, 
it can be shown that any real engine must have a lower efficiency than a Carnot engine, i.e. 

T) < 1 — T l /T h (real engines) Equation 74. 

Note also that “engine” does not necessarily mean something with pistons and cylinders; it refers to 
any device that generates work transfer (shaft work, electrical work, etc.) using heat transfer as the 
energy source {i.e. a heat engine in the vernacular of thermodynamics). 

Figure 24 shows some consequences of the First and Second Laws as applied to heat engines (the 
black box; what’s inside is irrelevant). 

No practical engine operates on a cycle similar to that of the Carnot cycle in which heat is added at 
constant temperature T H (sounds weird, adding heat at constant temperature, but it can be done in 
theory...) and rejecting heat at another constant temperature T L . Some other idealized 
thermodynamic cycles that more nearly approximate real cycles include 

• Otto Cycle , which is a model for spark-ignition reciprocating-piston engines like those in most 
automobiles. The cycle assumes an ideal gas with the following processes: 
o Isentropic compression (see page 62) by a volume ratio of r 
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o Heat addition at the minimum volume with no change in volume during the heat 
addition, and 

o Isentropic expansion by a volume ratio of 1/r to get back to the initial volume 
The ideal, theoretical efficiency of this cycle is given by 

rj = 1 —4r (ideal Otto cycle) Equation 75 

r Y 


which satisfies the function tests 

o T| = 0 when r = 1 (no compression, no net work done; you have a heater, not an 
engine, in this case) 

o r| —> 1 as r —> co (efficiency can never exceed 1) 

Crayton Cycle , which is a model for gas turbine engines. The cycle assumes an ideal gas with 
the following processes: 

o Isentropic compression by a pressure ratio of r 
o Heat addition at constant pressure, and 

o Isentropic expansion by a pressure ratio of 1/r to get back to the initial pressure 
The ideal, theoretical efficiency of this cycle is given by 

1 1 
*1-1- r (y-l)/y 

which satisfies the same function tests as the Otto cycle. 

Kankine Cycle , which is a model for steam turbine engines. The cycle assumes 

o Isentropic compression of liquid (water or whatever fluid) by a pressure ratio of r 
o Heat addition at constant pressure until the fluid is in the vapor (gas) state 
o Isentropic expansion of the steam by a pressure ratio of 1/r to get back to the initial 
pressure 

There is no simple expression for the efficiency of the Rankine cycle because it does not 
assume a fluid with a simple equation of state like an ideal gas. Note that there is no 
difference between the Brayton and Rankine cycles except for the type of fluid used; both 
assume isentropic compression, constant-pressure heat addition and isentropic expansion 
back to the starting pressure. 


(ideal Brayton cycle) 


Equation 7 6 
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COP Watts of heat 
transfer 
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small, i.e. COP < TJ(T H - T L ) 


Figure 24. Some consequences of the First and Second Farrs of thermodynamics as applied to heat engines and heat 

pumps. 
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The principle of increasing entropy can also be applied to pure heat transfer. For a device with no 
work input or output (W = 0), in order to satisfy energy conservation, Q H = Q L . Since the 2 nd law 
requires Q H /T H - Q L /T L > 0, with Q H = Q L , 


1/T h — 1/T l > 0 or T h > T l (heat transfer with no work transfer) 

which ensures heat transfer can only occur from a higher temperature to a lower temperature. Note 
that this not say that heat transfer can never occur from a lower temperature to a higher temperature, 
but that it cannot occur when there is no work transfer. 

Following along this same line, one can also consider the case opposite of heat engines, namely 
refrigerators that obviously do enable heat transfer from low temperature to high temperature. If 
there is work transfer W into the device, then Q H = Q r + W and in this case heat can flow from a 
lower temperature to a higher temperature if Q r is sufficiently small. In other words, the work 
transfer W (which causes no entropy production or loss) decreases the ratio of Q r /Q H so that it can 
be less than T L /T H so that Q h /Th — Q L /T L > 0 can be satisfied as required to have entropy 
production > 0. 

What is the best possible performance of a refrigerator? In this case the concept of “efficiency” 
doesn’t apply since work is an input rather than an output, but one can define a different figure-of- 
merit called the coefficient of peformance (COP) = Q L /W. For the best possible (Carnot cycle in 
reverse) refrigerator, the COP would be 

cop - what you get _ Ql _ Qh - w _ q h l _ 1 t _ i 1 _ t l 

What you pay for W W W f] Canwr , _ T l T h -T l 

T h 

and thus for any “real” refrigerator 

T 

COP <--— (refrigerator) Equation 77. 

Note that, as a function test, as T H approaches T L , the COP approaches infinity, since no work input 
(W = 0) is required to transfer heat across zero temperature difference. 

Note also that a heat pump used to heat homes is the same device as a refrigerator or air conditioner 
in that heat is transferred from a lower temperature to a higher temperature at the expense of some 
work input. The only difference is that in a refrigerator, Q L is the desired heat transfer ( from an 
object, to make it colder than ambient temperature) and Q H is waste (heat transfer to the 
surroundings at higher than ambient temperature, which is why your cats love to sleep behind your 
refrigerator), whereas with a heat pump Q H is the desired heat transfer (to your living room) and Q L 
is waste (heat transfer from the cold outside environment). Because of this, the definition of COP 
for a heat pump is different from that of a refrigerator. For the best possible (Carnot cycle in 
reverse) heat pump, the COP would be 
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and thus for any “real” heat pump 


COP < 


T 


H 


t h -t l 


(refrigerator) 


Equation 78. 


The advantage of a heat pump compared to simple electrical heating is that several units of heat 
transfer (Q H ) can be obtained for one unit of work transfer (W), whereas with simple electrical 
heating, only one unit of heat transfer is obtained per unit of work transfer. Of course, as the 
difference between T H and T L increases, COP decreases and thus the advantage of the heat pump 
decreases. Moreover, heat pumps, being functionally equivalent to air conditioners or refrigerators, 
are far more complicated that simple electrical heaters. For this reason heat pumps are not in 
widespread use, even in locations where air conditioners are installed and you could use the same 
device run forwards and backwards for both heating and cooling. (Keep in mind that heat pumps 
need electrical power whereas heating can also be done with natural gas which costs about Y 4 as 
much for the same energy.) 

Heat transfer 


The First Law of thermodynamics places restrictions on how energy can be converted from one 
form to another, and the Second Law places restrictions on the direction which processes may occur, 
but neither one says anything about how fast such processes occur. Here we’ll just look at the rates of 
heat transfer, which is only one piece of the puzzle. We’ve already talked about rate processes in 
terms of fluid mechanics, and you know something about dynamics (F = ma applied to a solid 
body.) If there are chemical reactions, we would need to compute their rates also, but that’s beyond 
the scope of this course. 

Heat transfer may occur by one or more of three forms: conduction , convection or radiation , which we’ll 
discuss separately below. 

Conduction 


Conduction heat transfer occurs in an immobile material (/.&, not a moving fluid) due to the vibrations 
of the molecules within the material. The more rapidly vibrating (warmer) portion of the material 
induces faster vibrations in the initially cooler part of the material and thus enabling randomly 
directed kinetic energy to pass through the material. The rate of said heat transfer is described by 
Fourier's Fan/. 


0 = -kA 


AT 

Ax 


Equation 79 


where Q is the rate of heat transfer (in Watts or some other unit of power), k is the thermal 
conductivity of the material, A is the cross-section area of the material exposed to heat transfer (i.e. the 
area in the direction perpendicular to the direction of the temperature gradient), AT = T H - T L is the 
temperature difference across the material (T H -= hot side temperature, T L = cold side temperature) 
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and Ax is the thickness of the material. AT can be specified in either absolute (K or R) or relative 
(°C or °F) units since the addition factor (273 from °C to K, or 460 from °F to R) will cancel out in 

the T h - T l term. Since Q is in units of watts, A is meters 2 , AT is degrees C or K and Ax is meters, 
the units of k must be W/m°C or equivalently W/mK. Note that the minus sign ensures that heat 
transfer is positive when AT is negative — in other words, heat must flow from high temperature to 
low temperature as required by the Second Law of Thermodynamics. Some typical thermal 
conductivities are given in Table 5. Note that all of these values are approximate because k depends 
on the temperature and composition of the material. In particular, pure metals have much higher 
conductivities than alloys (mixtures of metals) because in the case of metals the mobile electrons can 
transport thermal energy much faster than vibration within the solid structure itself can. 


Material 

Thermal conductivity (W/m°C) 

Air 

0.026 

Glass wool insulation 

0.04 

Plastics 

0.2 - 0.4 

Water 

0.6 

Glass 

1 

Silicon dioxide ceramic 

2 

Steel 

20-40 

Aluminum 

200 

Copper 

400 

Diamond 

2000 


Table 5. Thermal conductivities of some common materials at room temperature 


TLxample . The walls of a house are filled with glass wool insulation. The wall is 4 inches thick. The 
temperature on the inside wall is 70°F and outside it’s a nippy 0°F. The wall is 10 feet high and 20 
feet wide. 

a) What is the rate of heat transfer through this wall? 


0 = -kA 


AT 

Ax 


= 284VF 




3.281 ft) 


2 \ 

/ 


70°F -0°F 

' 1° C 

4 in 

V 1.8 °F 


'12//I 3.28lft N 
v ft m , 


\ 


/ 



If electrical heating is used at a cost of 10 cents per kilowatt-hour, what is the cost to heat 
this house for 1 day (just on the basis of heat loss through this one wall)? 


$0.10 kW 
kWhrlOOOW 


(284VF)(24/zr) = $0.68 


which doesn’t seem like much. But the real cost of heating is due to heat loss through windows, 
which are much thinner and made of materials with much higher k than glass wool insulation. 
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Convection 


Convection is heat transfer due to fluid flow across a surface. It is in general much faster than 
conduction (immobile material, no fluid flow) because in the case of convection, there is a 
continuous supply of hot fluid to deliver thermal energy to the cold surface, or of cold fluid to 
remove thermal energy from the hot surface. The rate of heat transfer by convection is given by 


Q = hAAT = hA(T sulface - T fluid ) 


Equation 80 


where h is the convective heat transfer coefficient (units Watts/m 2o C), A is the area of the surface exposed 
to convective heat transfer, T surface is the surface temperature and T fluid is the fluid temperature far 
away from the suface. In general it is very difficult to compute h because it depends on both fluid flow 
and conductive heat transfer (from the surface to the fluid adjacent to the surface) and thus involves 
Fourier’s Law coupled to the Navier-Stokes equations. For the purposes of this course I’ll just give 
you a value of h when needed. Some typical values are 10 W/m“°C for buoyant convection in air 
(when there is no forced flow, just the rising or falling of air due to a temperature difference 
between a surface and the surrounding air), 100 — 1,000 W/m 2o C for turbulent flow of water over a 
surface, and up to 10,000 W/m 2o C for heat transfer to boiling water. 


Example. Due to a tornado the insulation blew off the wall of the house in the previous example, so 
now there is a 70°F wall exposed directly to 0°F air. What is the rate of heat transfer now assuming 
buoyant convection only (i.e. the tornado has passed and the wind is calm) with h = 10 W/m 2o C? 


Q = hA(T suiface 





3.281 ft) 



= 1225W 


which is 25 times more than the insulated wall! 


Radiation 

Radiation heat transfer is heat transfer due to electromagnetic radiation between objects. Radiation 
what makes a fire feel warm even when you’re 10 feet away, i.e. too far for conduction or 
convection to be significant. The rate of heat transfer between two surfaces at temperatures T H and 
T r is given by 



Equation 81 


where O is the Stefan-Rolt^mann constant = 5.67 x 10’ 8 Watts/m 2 K 4 , 8 is the emissivity of the surface (a 
dimensionless number between 0 and 1; closer to 1 for opaque, non-reflecting surfaces and closer to 
0 for highly polished, reflecting surfaces). Note that temperatures must be specified in an absolute scale {i.e. 
Kelvins, not Celsius) since there is a () 4 term and thus the 273 or 460 conversion factor does not 
simply subtract out as it did with in the case of conduction or convection. In other words, for 
conduction and convection, 
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CT„+273) - (T l +273) = T h - T L 
so that either °C or K are acceptable, whereas for radiation 
(T h +273) 4 - (T l +273) 4 ^ T H 4 - T L 4 


so only K is acceptable. 


Example. For the unfortunate house above with a bare wall exposed to ambient air, what is the rate 
of heat transfer by radiation? The wall emissivity (b) is 0.5. 


Q = osA(T 4 -T 4 ) 

k -8 W 


= 5.67x10 


2tz4 


m z K 


(0.5) 


(10 ft x 20 ft) 


/ \ 2 \/ 

m 


3.28 1 ft 


(((70 + 460 )R) 4 - ((0 + 460)/?) 4 )( ^ 


= ni2w 


which is less than convection in this case. But note that since conduction and convection increase 
linearly with temperature, radiation increases with temperature to the fourth power ; thus at sufficiently high 
temperature, heat transfer by radiation will always exceed that due to conduction and convection. 
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Chapter 9. Written and oral communication 


“The single biggest problem in communication is the illusion that it has taken place .” — George Bernard Shaw. 
The Golden Rule of written technical communication can be stated as follows: 


If l didn’t already know this subject and were reading this paper (or listening to this 
presentation) for the first time y would I understand it? 

Written papers or reports 

1. The usual organization of a paper or report is 

A. Heading: title, authors, affiliations 

B. Abstract: explains what was done and what the main conclusions are. Must be short (a 
few hundred words at most, depending on the journal requirements), no matter how 
long and how complicated the paper is. 

C. Introduction: 

i. Explain what your problem is and why it is important. 

ii. State what is known about the subject. 

Hi. Complain about what is lacking in the current state of knowledge 

iv. Txplain what you will do that is better (may be in a separate Objectives section). 

D. Method: experimental apparatus, numerical model, whatever 

E. Results: what you found and how it compares with previous works 

F. Conclusions: what you learned 

G. Future work (optional) 

H. Acknowledgements (optional) — organizations that funded the work and/or people who 
helped but wasn’t included on the author list (e.g. people who gave advice but didn’t 
participate in the work itself, technical support, computer programmers.) 

Almost all novice writers (and many experienced writers) fail in two ways when organizing a paper. 
First, they fail to state clearly their objectives (what they are trying to learn) and their message or 
conclusions (what they found). Any piece of information that does not help to support the 
message doesn’t belong in the paper (unless it helps to show what isn’t certain about the 
conclusions.) For each paragraph, each picture, etc., ask yourself 2 questions: 

I. What is the message I am trying to get across? and 
2. Does this picture or text do that? 

In particular, students like to report on everything they did that went wrong before they got to the 
final results, just to show that they worked hard even if they didn’t accomplish much. This may be 
ok for a lab report in a class, but for a technical paper, nobody wants to know that the first 7 
voltmeters you tried didn’t work because someone spilled coffee on them. (Worst line ever: “We 
didn’t finish the project, but we have all the parts ...” Which means you have some parts that 
might be useful; you never know if you have all of them until after the project is done.) 

2. Every symbol in the text 

A. is defined in a Nomenclature section (preferred) or defined at its first appearance in the 
text (often you have to do it this way because of space limitations, but it’s annoying to 


94 



have to scan through a long document to find out what a is.) Hopefully you know by 
now how important it is to define your symbols. 

3. Every equation that is set apart from the text 

A. has a number 

B. has all of its symbols defined if not already or defined in a Nomenclature section (if used) 

4. Every word 

A. is spell-checked 

B. is defined the first time it is used if it is a “buzz word” or acronym. Example: “Many 
engineers use an excessive number of Three Letter Acronyms (TLAs).” 

5. Every figure 

A. is assigned a number — “Figure 1,” “Figure 2,” etc. 

B. is referred to in the text — no “orphan” figures! 

C. is referred to as “Figure x” if it appears at the beginning of a sentence, otherwise it is 
called “Fig. x” 

D. has a sensible scale on each axis (i.e. 0, 1,2, 3; not -0.37, 0.15, 0.67) 

E. When showing multiple plots of similar results, use the same scales. For example, if 
showing the burning velocities of methane-air and propane-air mixtures as a function of 
fuel concentration on separate plots, use the same scale for burning velocity on each 
unless they have drastically different ranges. 

F. Uses a logarithmic scale if a large numerical range of data (more than one decade) is 
covered (otherwise all the data having low numerical values are squashed together) 

G. has the units defined on each axis 

H. has a caption (in addition to its figure number) 

I. has all relevant conditions (pressure, temperature, whatever is important) stated on the 
plot or in the caption 

J. has all plot symbols (squares, circles, filled or open, ...) and curves (solid, dotted, dashed,..) 

defined either in a legend box within the figure (preferred method) or in the caption 

K. Does not have a lot of “white space” 

L. if it is a picture, it has a scale on the picture or has a statement in the caption such as 
“field of view is xxx cm by yyy cm” or has some object of easily identifiable scale (a 
person, coin, etc.) in the picture 

M. must be readable - caution on pictures!!! 
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Equivalence ratio 

Figure 25a. Terrible figure. 

Figure 25a shows a terrible figure with many common mistakes. What’s wrong? 

• The scales on each axis are terrible — weird numbers, not 1, 2, 3, ... 

• Units are not defined on the vertical axis (Seconds? Nanoseconds? Millenia?) 

• The plot symbols are defined using meaningless notation (“Condition 17” means 
nothing to the reader.) 

• There is a tremendous amount of “white space” 

• Most of the data squashed together because a linear scale was used - the scale has to be 
large enough to cover the large values of rise time in “Test —117”, which goes up to 300, 
but most of the data is in the 10 — 50 range 

• There are tick marks inside, outside, all over the place (I prefer tick marks on the inside 
only). Also, the major and minor tick marks are the same length so it’s hard to 
distinguish between them. 

• The plot symbols are too small to see 

• The numbers are too small to read 

• All of the grid lines make it hard to read the data and legend. (I don’t like grids at all, 
they clutter the figure — if someone really wants to pick points off your graph, they can 
draw their own grid lines or ask you to email the data file to you.) 

• There are ugly looking jagged lines connecting the data points (rather than a smoothed 
curve) 
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Three data sets have lines connecting them, whereas the fourth does not (is there 
something different about the fourth data set that makes it ineligible for connecting 
lines?) 

The axes and tick marks are too thin 


big part of the problem is that most people just let their plotting program make bad plots using all the default 
settings , and somehow try to rationalise that still a good plot. Figure 25 b shows a reasonable figure presenting 
exactly the same data as Figure 25a. 



6. Every reference cited in the text 

A. appears in the reference section 

B. is a plain number (i.e. 11, 12, 13; not 11, 11a, lib) or follow the Harvard system (e.g. 
Smith and Jones, 1953) depending on the instructions to authors 

C. if a number, may be superscript or in [brackets] or (parenthesis) depending on the 
instructions to authors 

7. Every reference in the reference section 

A. is called out in the text; it is not acceptable to simply have a list of references at the end 
of a document without referring them in the text so that the reader knows what 
information was used from that reference 

B. has the journal name or book title (journal titles may or may not be abbreviated 
depending on the instructions to authors) 

C. has the page number (may be just the first page of the article or inclusive pages 
depending on the instructions to authors) 

D. has the journal issue number (if a journal article) 

E. has the publisher (if a book) 

F. has the year of publication 

G. may or may not have the title of the article depending on the instructions to authors 
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Bottom line: ask yourself, if I were reading this paper for the first time, and I were not 
already aware of the results, would I understand this paper??? 


Oral presentations 

Most of the rules listed above for written presentations apply to oral presentations. In particular: 

1. Organization 

A. Title page title, authors, affiliations, acknowledgements (optional) 

(no abstract) 

B. Introduction: 

i. Explain what your problem is and why it is important. 

ii. State what is known about the subject. 

Hi. Complain about what is lacking in the current state of knowledge 

iv. Cxplain what you will do that is better (may be in a separate Objectives section). 

C. Method: experimental apparatus, numerical model, whatever 

D. Results: what you found and how it compares with previous works 

E. Conclusions: what you learned 

F. Future work (optional) 

2. Every symbol, buzz word, acronym, etc. must be defined the first time it is used (no 
Nomenclature section) 

3. Equations aren’t numbered 

4. All of the rules for figures and pictures still apply 

5. References may be mentioned, especially if there are key works that your work builds 
upon or refutes, but are not numbered 

There are also special rules for presentations: 

6. Use a laptop-based powerpoint presentation. This makes it much easier to combine/split 
previous presentations, add color, animations, sound effects, etc. But the most valuable 
aspect is probably that it allows you make last-minute changes. Also it is useful because then 
you can email the presentation to interested people, or post it on your website. 

7. Do not use 5-point font! Reduce the amount of material presented and use big fonts! 
Make sure everything is legible. A. good rule of thumb is that if the slide is printed on standard 8.5” x 
11 ” paper, you should be able to put the page on the floor and read everything on the page while standing up 
and looking down at the page. 

8. Use color. The human eye is much more sensitive to variations in color than shades of 
gray. Key words can be given emphasis using colors. Plots in journals usually have to be in 
black only, but in a presentation you can use colors to make it easier for the audience to 
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distinguish between different data sets (as in Figure 25b above) that it is to distinguish 
between data sets by looking at the different symbols. 

9. Include movies. Why limit yourself to static presentations when you have the power of a 
computer? A picture is worth a thousand words, and a movie is worth a thousand pictures. 

10. Address the audience. Say things like, ’’this plot shows you the effect of x on y...” 
rather than ’’this plot shows the effect of x on y...” 

11. Keep reminding the audience of your nomenclature. That is, if you show an equation 

I' _ 2 

F — me 

don’t say ’’this equation shows that eee equals emm cee squared,” (the audience can already 
see that). Instead say, "this equation shows that the energy of a substance is equal to its 
mass times the speed of light squared” (the audience has forgotten your definitions of E, m 
and c that you gave 12 slides back). 

Bottom line: ask yourself, if I were in the audience listening to this presentation for the first 
time y would I understand it??? 
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Appendix A. Suggested course syllabus 

Below is a “template” for a semester-long course syllabus, patterned after my own 


AME 101 — Introduction to Mechanical Engineering and Graphics - Fall 2010 

Lecture: Tuesday and Thursday 8:00 - 9:20 OR 9:30 - 10:50 am, ZHS 252 

Labs: Tuesday OR Thursday, 12:30 — 1:50 pm, SAL 126 (Tuesday only), SAL 127 (both days) 

Final exam: 8:00 AM class: Tuesday, Dec. 14, 4:30 pm — 6:30 pm 

9:30 AM class: Thursday, Dec. 9, 11:00 am — 1:00 pm 
Note: I will attempt to coordinate a single final on Dec. 14. 

Web page: Accessible to registered students through Blackboard at http:/7blackboard.usc.edu . A 
direct link to the syllabus and other information is http:/7ronney.usc.edu/AME101F10/ 

Instructor: Paul Ronney 

Office: Olin Hall 430J 
Phone: 213-740-0490 
Email: ronney@usc.edu 

Office hours: Wednesdays 1:00 — 3:30 pm; other times by appointment. 

Teaching Assistants: 

Shalini Reddy (sreddy@usc.edu) ; Office hours Thursdays 2:00 pm — 5:00 pm 
Qianyu (“Cherry”) Du ( qianyuli@usc.edu) ; Office hours Fridays 9:00 am — 12:00 noon 

Grader: To be announced 

Texts: 

• Lecture notes 

• Handouts in laboratory sessions 

• (OPTIONAL) An Introduction to Mechanical Engineering by Jonathan Wickert, 2 nd edition, 
Thomson-Engineering Publishing, 2005, ISBN-10: 0534552978; ISBN-13: 978-0534552978 

• (OPTIONAL) Introduction to Solid Modeling Using SolidWorks 2010, by William Howard and 
Joseph Musto, McGraw-Hill, 2009, ISBN-10: 0073375438; ISBN-13: 978-0073375434 

Grading: __ 


Homework 

20% 

Design projects & competitions 

15% 

Laboratory (graphics) 

30% 

Midterm exams (2) 

10% each 

Final exam 

15% 


• Breakdown of laboratory grade 

o 4 homeworks (60% of lab grade) 
o 1 mini-project (40% of lab grade) 

• NO LATE HOMEWORK WILL BE ACCEPTED, PERIOD, NO EXCEPTIONS in 

either lecture or lab. 


















Graphics Laboratory 


The Graphics Laboratory aspect of AME 101 will introduce you to a powerful Computer Aided 
Design (CAD) tool, SolidWorks, which is widely used in industry today. As an introductory course, 
it is not intended to make you an expert with this software; however, you will acquire a basic 
knowledge of CAD skills extensively used in mechanical engineering today. 

This is a hands-on, learn-by-doing class and all instruction will require active use of software, 
SolidWorks 2009-2010, which is available in all the ISD-managed computer laboratories. This 
software is also available for installation on your own machines should you desire to work at home. 
Detailed instructions for home version installation are posted on the Blackboard 
(https://blackboard.usc.edu/) . 

In lieu of a formal course textbook, a short presentation will be posted on Blackboard each week 
before the class. The presentation will cover the material for the week and conclude with one or 
more tutorials and/or exercises. The tutorials and exercises are designed to show you how certain 
tasks may be accomplished and allow you to practice either on your own or in the lab sessions where 
help will be available. Your mastery of the material will depend entirely on how much you work 
with the software. 

The lab homeworks must be submitted to the grader with the following attachments: 

• Template “cover sheet” (see below) including name, date and images of the problems 

• SolidWorks files of the assignment 

AME 101 Graphics Lab Homework Template 
LAB HOMEWORK ASSIGNMENT #1 


Name: _ Lab section:_ Date: 

• Problem 1 
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Tentative schedule 

“Plans are nothing... planning is everything” — Dwight D. Eisenhower 


Week 

Monday 

Date 

Lecture subject 

Graphics lab 
subject 

Tues. 

lecture 

Thurs. 

lecture 

Assignment 

due 

1 

8/23 

Introduction, units 

Introduction 

PDR 

PDR 


2 

8/30 

Units 

Sketch basics 1 

PDR 

PDR 


3 

9/6 

Engineering scrutiny 

Sketch basics 2 

PDR 

PDR 

Ll 

4 

9/13 

Excel for engineers; 
statistics 

Sketch basics 3 

PDR 

PDR 

G1 

5 

9/20 

Forces and moments 

on structures 

Feature basics 1 

PDR 

PDR 

L2 

6 

9/27 

Forces and moments 

on structures 

Feature basics 2 

PDR 

PDR 


7 

10/4 

Forces and moments 

on structures 

Feature basics 3 

PDR 

PI 

L3 

8 

10/11 

Materials and 

stresses 

Stress/ Simulation 

PDR 

Qi 

G2 

9 

10/18 

Materials and 

stresses 

Stress/ Simulation 

PDR 

PDR 

R1 

10 

10/25 

Fluid flows 

Assembly 

PDR 

PDR 

L4 

11 

11/1 

Fluid flows 

Assembly 

PDR 

PDR 

G3 

12 

11/8 

Energy and thermal 
systems 

Motion 

PDR 

PDR 

L5 

13 

11/15 

Energy and thermal 
systems 

Motion 

PDR 

Q2 


14 

11/22 

Energy and thermal 
systems 

Drawing 

P2 

XXX 

G4 

15 

11/29 

Energy and thermal 
systems 

Drawing 

PDR 

PDR 

L6 





XXX 

Final (12/9 
or 12/14) 

GP, R2 


Legend for schedule 

PDR 

SL 

Qn 

Ln 

Gn 

GP 

Pn 

Rn 

XXX 


PDR lectures 
Substitute lecturer 
Midterm exam n 
Lecture homework n due 
Graphics lab homework n due 
Graphics lab project due 
Design project n contest 
Design project n report due 
Break / holiday / end of semester 


Homework topics 

1 Units 

2 Scrutiny, Excel, statistics 

3 Forces & torques 

4 Materials and stresses 

5 Fluid flows 

6 Energy and thermal systems 
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Design projects 

1 King of the Hill (force, torque, mechanical design) — week 7 

2 Spaghetti bridge (materials, stresses) — week 14 

Note: Design teams will be assigned at random and different for each project in order for students 
to become better acquainted with each other and to avoid the “A-list, B-list, C-list” group dynamics 


(Possibly) useful information and disclaimers 

1. Exams will mirror lectures; be sure you understand the lectures. Please ask questions inside 
and outside class! (If you choose to buy the textbook, please understand that it’s just an 
additional reference, not something that I will follow closely.) 

2. I will call on students in class. This is not a popular practice with students, but I do it 
anyway because (a) it encourages students to attend class (though I don’t take attendance); 
(b) it encourages students to pay attention in class and (c) it helps me to get to know the 
students, and the students to get to know each other by name - many of you will be together 
for 4 years, so why not get acquainted now? 

3. This is my sixth time teaching this course, and since (in my humble opinion) all available 
textbooks are weak at best, my plan is to turn my lecture notes into a textbook. Thus, 
constructive suggestions are most welcome! Note: the lecture notes are still a work in 
progress, so there will be changes. I will do my utmost to inform students of changes and 
make the updating as painless as possible. 

4. This course is sort of like engineering boot camp; not always popular but students do come 
back in a year or two and tell me that what they learned in this class was useful and made 
their subsequent classes easier. 
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Class objectives 


• Furnish you with some basic tools of engineering 

• Units — English and metric system 

• “Engineering scrutiny” 

• Statistics 

• Approaches to problem-solving and teamwork 

• Provide introductory knowledge of engineering topics 

• Forces and torques 

• Fluid flows 

• Materials and stresses 

• Thermal and energy systems 

• Provide introductory knowledge of engineering graphics (laboratory section) 

• Solid modeling 

• Views and shading 

• Dimensions 

• Fillets, rounds, patterns 

• Assemblies 

• Computer Numerical Control (CNC) milling 

• Retention-related objectives 

• Provide a “roadmap” of why subjects you will be learning, and what will you do in the future 
with the knowledge gained 

• Making an intelligent choice of major - make your first engineering class a positive enough 
experience that you make a choice based on knowledge, not fear or intimidation 

• Develop confidence in your ability — “pride of ownership” of knowledge gained 

• Topics NOT covered in this class (but should be) 

• Electrical circuits 

• Ethics (covered to some extent in WRIT 130 and 340) 

• Computer animation (covered in AME 308) 

• History of engineering 

• Philosophy of engineering 

• Written and oral reporting 

Hidden agenda: To start teaching you to think like engineers. Over and over, engineering faculty 
hear from practicing engineers and corporate recruiters words like, “teach the students how to think 
and we’ll teach them the rest.” 

“You come in here with a skull full of mush and if you survive you leave thinking like a 
lawyer” - Actor John Houseman, portraying Harvard Law School Professor Charles 
Kingsfield in The Paper Chase (1973). [Substitute ‘engineer’ for ‘laywer’.] 
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USC and the Viterbi School 


Why USC engineering? 

• Aggressive, proactive leadership — buildings, rankings 

• Engineering has a high priority from the USC central administration 

• Student services and programs (Merit research, work study, counseling and tutoring, 
professional organizations, under-represented group organizations, ...) 

• Breadth of courses and escape routes for those who decide engineering is not in their future 

• Class sizes and faculty to student ratios 

• But it’s up to you to take advantage of all the opportunities and not develop “senioritis” in 
your freshman year 

USC Viterbi School of Engineering mission statement 

“The School of Engineering seeks to provide undergraduate and graduate programs of instruction 
for qualified students leading to academic degrees in engineering; to extend the frontiers of 
engineering knowledge by encouraging and assisting faculty in the pursuit and publication of 
research; to stimulate and encourage in its students those qualities of scholarship, leadership, and 
character that mark the true academic and professional engineer; to serve California and the nation 
in providing for the continuing education of engineering and scientific personnel; and to provide 
professional engineering leadership in the solution of community, regional, national and global 
problems.” 

Who’s in charge here??? 

1. The USC Board of Trustees has the ultimate say in what happens on campus. “As a private 
corporation, USC is governed by a board of trustees which has approximately 50 voting 
members. The board is a self-perpetuating body, electing one-fifth of its members each year 
for a five-year term of office.” 

2. President Max Nikias, Professor of Electrical Engineering - Systems — sets policy and directs 
others to execute that policy — not unlike the role of the U.S. President 

3. Provost (to be named) — the single person most responsible for making the vision of the 
President actually happen — role similar to that of “chief executive officer” of a corporation 

4. Dean of Engineering Yannis Yortsos — overall responsibility for the operation of the School 
of Engineering 

o Senior Associate Dean for Academic Affairs John O’Brien (Professor of Electrical 
Engineering — Electrophysics) — responsible for the integrity and operation of the 
academic program including teaching, accreditation, promotions and tenure of 
faculty, etc. 

o Associate Dean, Admissions and Student Affairs Louise Yates — you know her 
o 7 other Associate Deans — see http://viterbi.usc.edu/about/administration/ 

• Chairman Geoff Spedding, Department of Aerospace and Mechanical Engineering (AME) — 
overall responsibility for the operation of AME 

• AME faculty — 27 and growing 

• AME students — — 125 freshman - In what ways are you in charge? 

o Participate in aforementioned activities 
o Teaching evaluations 
o Directed research 
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o (Someday) alumni activities 
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ABET 


Engineering programs are accredited by the Accreditation Board for Engineering and Technology 
(ABET) (http:/7www.abet.org) . Each course is expected to have a “course objective” and a list of 
“course outcomes.” At the end of the semester, there will be a survey passed out to all students 
asking to what extent (on a 1 — 5 scale) the course outcomes were or were not met. 

Course objective for AME 101: 

To introduce the student to the science and art of Mechanical Engineering by providing (1) basic 
tools of engineering practice, (2) introductory knowledge of engineering topics, (3) facility with 
Computer-Aided Design software and (4) a perspective on how the large number of subjects 
covered in the mechanical engineering curriculum are inter-related. 

Course outcomes for AME 101: 

By the end of the course, the student will 

1. Understand the courses required for his/her Mechanical Engineering education at USC and 
why these courses are useful 

2. Understand and manipulate the units of engineered systems 

3. Scrutinize a calculated or measured result for “obvious” mistakes 

4. Be able to work productively as part of an engineering team working toward a common 
objective 

5. Create simple 2-D and 3-D models of parts and assemblies using Computer-Aided Design 
(CAD) software such as Solid Edge 

6. Have a basic understanding of the forces and torques on rigid, solid objects 

7. Have a basic understanding of engineered materials and the stresses they can withstand 

8. Have a basic understanding of the flow of fluids and the forces they exert on structures 

9. Have a basic understanding of thermodynamics, in particular application of the principle of 
conservation of energy to very simple systems. 

10. Have a basic understanding of the three modes of heat transfer and be able to apply the 
basic equations of heat transfer to very simple systems. 


ABET Program Objectives 

In addition to course-specific objectives and outcomes, ABET also specifies a set of “Program 
objectives” which are broad statements that describe the career and professional accomplishments 
that the program (in your case, Mechanical Engineering at USC) is preparing the graduates to 
achieve. For all engineering disciplines, the Program Objectives are: 

1. Graduates will be professionals working in engineering or in related areas such as computer 
science, business, law, medicine or public service, at both large- and small-scale businesses. 

2. Graduates will engage in lifelong learning, such as continuing their education through 
graduate school or professional development courses. 

3. Graduates will make use of modern and cutting-edge tools, such as advanced computer 
software and state-of-the-art laboratory equipment. 
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4. Graduates will be both competent technical innovators and industrial leaders. 

5. Graduates will incorporate societal, ethical and environmental considerations into technical 
decisions. 

6. Graduates will effectively communicate and work with persons and teams of diverse 
technical and non-technical backgrounds. 

ABET Program Outcomes 

Again at the “Program” level, ABET also specifies a set of “Program Outcomes” which are 
narrower statements that describe what students are expected to know and be able to do by the time 
of graduation. For all engineering disciplines these Program Outcomes are that the student should 
have 

1. an ability to apply knowledge of mathematics, science, and engineering 

2. an ability to design and conduct experiments, as well as to analyze and interpret data 

3. an ability to design a system, component, or process to meet desired needs within realistic 
constraints such as economic, environmental, social, political, ethical, health and safety, 
manufacturability, and sustainability 

4. an ability to function on multidisciplinary teams 

5. an ability to identify, formulate, and solve engineering problems 

6. an understanding of professional and ethical responsibility 

7. an ability to communicate effectively 

8. the broad education necessary to understand the impact of engineering solutions in a global, 
economic, environmental, and societal context 

9. a recognition of the need for, and an ability to engage in life-long learning 

10. a knowledge of contemporary issues 

11. an ability to use the techniques, skills, and modern engineering tools necessary for 
engineering practice. 

For Mechanical Engineering, the USC AME department has developed a more specific set of 
Program Outcomes, name that the student should have: 

1. a knowledge of chemistry and calculus-based physics with depth in at least one 

2. an ability to apply advanced mathematics through multivariate calculus and differential 
equations 

3. a familiarity with statistics and linear algebra 

4. an ability to work professionally in both thermal and mechanical systems areas including the 
design and realization of such systems 

5. (Petroleum concentration only) a knowledge of petroleum engineering topics 
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Appendix B. Design projects 

Generic information about the design projects 

How to run a meeting (PDR’s philosophy...) 

Every meeting must have three things, all in writing: 

• An agenda. What is it that needs to be discussed at the meeting? If it isn't written down, 
some items will be forgotten or will get dropped as the meeting runs over its time limit, so a 
written agenda is usually needed. Sometimes everyone knows the agenda items (as in a 
weekly meeting, for example), or the list of things to be discussed is very short, so a written 
agenda isn't needed 

• Minutes. What was said and what was decided at the meeting? There definitely needs to be a 
permanent record of this, because you WON'T remember a week later what was said or 
what was decided. (More likely, you will remember but your recollection will be different 
from everyone else's.) 

• Action items. Who will do what as a result of the meeting? When is it needed? What will 
people do that is different than what they would have done without the meeting? Think 
about that last question — if no one is going to do anything differently as a result of the meeting, what was 
the purpose of the meeting? 

If you don't have all three of these items, then you have to ask yourself, why did you meet? What 
were you trying to accomplish by meeting? Was it a meeting or just a party, seminar, etc.? 

Many groups choose to start a Facebook (or Myspace, or Yahoo) page just for their project; this is 
acceptable and in fact encouraged; clear, swift and accurate communication is of importance even in 
small projects like this one and it large, real-world engineering projects it is absolutely essential. 

Suggestions for the written report 

Please no binders for these reports, just pages of plain paper stapled together; all those 3- 
ring binders take up too much space! 

The report should include 

1. Cover page - a title page is nice to have on a report (yeah, it's a waste of a sheet of paper, agreed, 
but it looks nice.) I do succumb to enjoying project titles with a twist, e.g. “playing with your 
food” for the spaghetti bridge project. It breaks up the monotony of grading all those reports! 

2. Table of contents - Some sort of organization of the pages (e.g. chronological, by type of 
document (meetings, test results, postmortem) makes it easier to read and understand. 

3. Body of the report including 

a) Statement of objective(s) 

b) Drawings of preliminary design concepts and critiques of these designs 

c) Test data for preliminary designs 

d) Explanation of why you chose your final design 

e) Construction of your device 
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f) Results of the “official” test 

g) What you would do differently if you built another device 

4. Optional appendices(s) including (these may be incorporated into the body of the report rather 
than an appendix if you feel it improves the organization and readability of the report.) 

a) Meetings: agenda, minutes, action items 

b) Email exchanges 

c) Pictures of your test apparatus, construction and the final device 

d) Anything else you think is appropriate 

General comments about the reports: 

• Printouts look a lot nicer than handwritten notes. I really can't do any writing by hand any 
more; my handwriting is illegible even to myself and even worse, my thought processes have 
become so jumpy I can't write anything, even one paragraph, from start to finish. 

• Pictures are also nice. I think a picture is worth way more than 1000 words, because who wants 
to read 1000 words? (For PowerPoint presentations, I have a corollary - a video is worth a 
thousand pictures.) But, a collage of pictures at the end of the report isn’t very useful. 
Include the pictures in the body of the report, and every picture needs a caption so the reader 
knows what the picture shows and what the reader is supposed to learn from looking at it. 

• Make sure throughout that whatever meetings, background research, testing, etc. you do, you are 
focused on the specific objectives of the project. For example, in the spaghetti bridge project, 
just stating that such-and-such an idea is a good one because it should make the bridge strong is 
not very useful. Remember that the objective in that case was a high strength to weight ratio, 
not just high strength, so everything you do should revolve around that fact. 

• Just showing sketches of ideas isn't terribly useful unless you then explain which ideas you 
embraced, which you rejected and why. (In a more formal setting you would document all the 
team members opinions and ideas, since this may lead to patents later on.) Also, just showing a 
bunch of stuff downloaded from the internet to pad the report isn't useful. Think of a report 
as if you were making a movie - you have a story to tell, you want to tell it in the most 
compelling way possible, and anything that doesn’t help tell the story should be left on 
the cutting room floor. 

• MOST IMPORTANT POINT. What convinces me more than anything else that (independent 
of the outcome of the contest) you are serious about the project is the TESTING part. Think 
about it - if you don't test anything, it says you're just going to accept whatever the first attempt 
at the design/build/test cycle gives you. And if you don't document the testing, did you really 
do any testing? (A variation of the proverbial question, "if a tree fell in the forest and no one 
heard it, did the tree make a sound?") So presenting test results is HIGHLY ENCOURAGED. 
Test results can be from either physical testing (e.g. spaghetti strength) or modeling (e.g. 
SolidWorks stress analysis of your spaghetti bridge). 

• Reference your statements. For example, if you say, “our research has shown that spaghetti is 
stronger in tension than compression,” state what test you did or what the source is from which 
you obtained that information. (Note: a comment by Joel2345 on blogger.com doesn’t qualify 
as a legitimate source.) 

• Make backup copies of your data. I'm not very sympathetic to "my computer crashed," which is 
the modern equivalent of "the dog ate my homework." All hard disk drives have anxiety and 
deadline sensors and are pre-programmed to crash when your anxiety level is highest or the 
deadline is nearest. So put your data on USB drives, Dropbox folders or just email everything 
back and forth to your partners so there are no "single-point failures" regarding data security. 
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Candle-powered boat 


Objective 

The goal of this contest is to design and construct a boat that uses a single candle as an 
energy source for propulsion and travels as fast as possible using only this energy source. You will 
race your boat against those built by other teams. 

Design rules 

1. The boat will use a single “tealight” candle, supplied by the “race officials/ 5 similar to these: 

http://www.walmart.com/ip/Light-In-the-Dark-Unscented-Tealight-Candles/22697463 . No 

other stored energy sources (e.g. other fuels, batteries, rubber bands, gravitational 
potential energy, “general purpose heat sources”, etc.) are allowed. You may use 
naturally occurring wind or solar power instead of or in addition to the candle, however, 
there is no guarantee that there will be any direct sunlight or wind on the day of the contest. 

2. You may use a “Pop Pop Boat 55 type of design, and in fact you may use a store-bought boat, 
however, if that’s all you do you won’t be getting a very good score on the project report due 
to a lack of creativity, originality and level of effort. I do, however, suggest buying a Pop 
Pop Boat as a starting point for your own design. There is no requirement, however, that 
you use a Pop Pop type of boat. For example you could use the candle and sunlight to heat 
a thermoelectric module on one side, use the water to cool the other side, and use the 
electricity generated by the thermoelectric module to drive a motor / propeller system. 

3. You can build and use as many different boats as you wish on race day; you do not need to 
use the same boat for each race, but of course in any given race you can only use one boat. 
It is suggested that you bring at least one back-up boat on contest day so that if your “A- 
team” boat fails you’ll at least have something else to enter in the remaining races. 

4. The boat must be autonomous. No remote power, control wires or radio-control links 
allowed. 

5. The boat’s exterior dimensions must be no larger than 1 foot long x 6 inches wide and have 
a draft (maximum depth under the water’s surface) of 6 inches. There is no limit on the 
height or weight of the boat. The “boat” does not even have to float, but keep in mind that 
the depth of the fountain is not constant (it varies from about 6 inches to 18 inches) so a 
wheeled vehicle that touches the bottom would need to have adjustable wheel struts. 

6. To travel in a straight line from the starting line to the finish line the boat may connect to a 
guide string (provided by the judges) that will be strung about 3 inches above the water’s 
surface. The height of the string above the water line is not known exactly because it 
depends on the amount of water in the fountain that day; it may be anywhere between 2 and 
6 inches, so you’ll need to have a way of adjusting your guide mechanism depending on the 
string height on the day of the contest. It is up to the team to decide how to connect to the 
string and avoid getting snagged by the string. 

7. Time will be scheduled in BHE 310 for groups to use the laboratory facilities and tools for 
making and testing their vehicles. The lab gets VERY crowded on the last work day! I 
strongly recommend that you do at least some of your building and testing during the earlier 
work days. The rules for using the BHE 310 lab are as follows: 

• You can only work in the BHE 310 lab during scheduled work hours 

• You must have completed the safety training given during the first few weeks of class 

• You must wear closed-toed shoes to enter the lab 

• You must put tools back where you found them when you’re done! 
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• There are some bins of materials (wood, aluminum, plexiglass, etc.) as well as wire, 
solder, tape, etc. are available in BHE 310 but DO NOT use any materials in the 
wooden drawers assigned to the AME 441 groups! There is no guarantee that any 
particular type of material will be available for you. 

• “Official” candles will be available for testing, however, fresh candles will be given on 
race day. This is to ensure that groups do not “juice up” their candles with gasoline, 
gunpowder, plutonium-238, etc. As a result your design must allow for installation of a 
fresh candle on race day. You may continue to use the same candle for all contests or 
request a fresh one. 

The contest 

1. The contest will be held on Thursday October 9 at 8:30 AM (i.e., during class time) at the 
fountain in the E-quad. Since many of you have classes during the time block (8:00 am — 
9:30 am or 9:30 am — 11:00 am) when you don’t have AME 101, your group may need to 
staff the contest in shifts. 

2. Boats will compete in a 5-round “Swiss system tournament” similar to that typically used in 
chess contests (see for example http://en.wikipedia.org/wiki/Swiss system tournament) so 
that every team will participate in the same number of rounds. Modifications to the boat(s) 
are permitted between (but obviously not during) runs. 

3. Boats will be allowed up to 2 minutes to travel across the fountain. The first one to reach 
the other side is the winner of that contest. If neither boat reaches the finish line in the 
allotted time, the boat closest to the finish line will be declared the winner of the contest. 

4. Three 1-on-l contests will be run at the same time, on the west, south and east sides of the 
fountain. The lengths of the three “race courses” are about 25 feet but not exactly the same. 
The north side of the fountain can be used for testing. 

5. The order of the 1-on-l contests will be selected at random and published in advance. Rigid 
adherence to schedule is necessary to conduct all 55 or so required contests during the 
available time. After the previous contest is completed, you will have 3 minutes to ready 
your vehicle for the event. If you are not ready at the time of your event, your team will be 
disqualified for that specific 1-on-l contest. 

6. A team member will place the boat in the water and perform any required start-up 
operations (e.g. lighting the candle, adjusting the heat exchanger). To ensure fairness, an 
impartial “race official,” not a team member, will release the boat at the start of the race. 
The only thing the “race official” will do is release the boat, e.g. not flip a switch or turn a 
valve. Contestants may not touch the boat once the race begins. If the candle goes out, a 
race official will relight it, but only one relight per race will be allowed. 

7. Project reports will be due on Monday, October 20, at 4:30 pm in OHE 430N. 

Organization and grading 

You will work in teams of 4 people, assigned at random, nominally with 2 people from each 
section of the class. Each team must write a report about their work, with one report per team, not 
per student. The report will be the primary means of grading the projects (2/3). Your level of 
success in the competition also counts (1/3) according to the formula: 

Competition score = 50 * (N-P)/(N-l) + 50 

where P is your place at the end of the Swiss tournament and N is the total number of teams 
competing. So even the last place team (N = P) gets 50% on the competition score. Also, the 
report score counts more than the competition score; in other words more weight is given to the 
process than the result, which is just the opposite of how life really works where you’re judged 100% 
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on the result and 0% on the process. 
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King of the Hill 

Objective 

The goal of this contest is to design and construct a vehicle that can climb a “hill” under its 
own power, stop at the top of the hill and defend its position against an opposing vehicle coming up 
from the other side of the hill. The “hill” will be represented by two 5 foot long x 1 foot wide 
ramps separated by a 1 foot x 1 foot mesa at the top. The hill altitude is 33 inches from base to 
summit. This is a very steep hill, but the ramps are covered with a non-slip material like you might 
find in a bathtub. Experience has shown that the most important thing is to get your vehicle 
to the top. Many vehicles fail because they don’t have enough torque (i.e. the wheels or tracks don’t 
move) or traction (i.e. the wheels or tracks spin but the vehicle doesn’t move) to climb the steep 
ramp. Offensive or defensive tactics for use on top of the ramp are of secondary importance. I 
strongly recommend you don’t buy a radio-controlled racecar and cannibalize it, almost all of them 
are designed for speed, not hill-climbing ability. However, experience also shows that possession is 
9/10 of the law and usually the vehicle that reaches the top of the ramp first wins because it’s 
difficult for a moving vehicle to displace a stationary one. So you’ll have to decide on some optimal 
combination of speed, torque, traction, etc. 

A suggested source of motors, gearboxes, battery boxes, tracks, etc. is 
http://www.pololu.com , particularly the products made by Tamiya, for example these gearboxes: 
http://www.pololu.com/catalog/product/61 or http://www.pololu.com/catalog/product/74 and 
these tracks and sprocket sets: http://www.pololu.com/catalog/product/415 , 

http:/ /www.pololu.com/catalog/product/224 . Tilt switches to shut your vehicle off when it 
reaches the top of the hill can be obtained from many sources. Keep in mind that there are tilt 
sensors (like your iPhone has) which can carry almost no current, you need a tilt switch. Check the 
specifications of your switch to make sure it can carry the current your motor draws. (How much 
current does your motor draw? Duh, read the motor specifications!) Another source for pulleys, 
belts, gears, nuts and bolts, etc. is http://www.smallparts.com . A fairly complete kit (motor, 
gearbox, switch, treads, chassis) is the “Tamiya 70108 Tracked Vehicle Chassis Kit” which you can 
find online. Fry’s Electronics stocks the Tamiya kit and tilt switches. Peter Ronney (age 10) will be 
an unofficial contestant using this kit. 

Time will be scheduled in BHE 310 for groups to use the laboratory facilities and tools for 
making and testing their vehicles. Also, the “hill” may be inspected and test runs performed at that 
time. The lab gets VERY crowded on the last work day! I strongly recommend that you do at 
least some of your building and testing during the earlier work days. The rules for using the BHE 
310 lab are as follows: 

• You can only work in the BHE 310 lab during scheduled work hours, but the ramps are 
available for test runs any time the BHE building is open 

• Basic hand tools are available (screwdrivers, pliers, wire cutters, etc.) but put the tools back 
where you found them when you’re done! 

• Bins of materials (wood, aluminum, plexiglass, etc.) as well as wire, solder, tape, etc. are 
available in BHE 310 but DO NOT use any materials in the wooden drawers assigned to the 
AME 441 groups! 

• You can use the machine tools (drill press, lathe, bandsaw, grinder, etc.) but you have to be 
checked out by the supervisor on duty (PDR, AME 101 TA or AME lab technician.) Of 
course you MUST wear safety goggles when using any machine tool. 
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Design rules 

1. The vehicle must be completely autonomous. No remote power, control wires or wireless 
remote-control links are allowed. 

2. The vehicle's exterior dimensions at the start of each run must not extend beyond the sides 
of an imaginary cube 1 foot on a side. A device, such as a ram, may extend beyond this limit 
once activated, but cannot be activated before the start of the run. 

3. The vehicle's total mass must not exceed 1 kg. 

4. The vehicle can start either on the ramp or at the base, on flat ground, your choice. The 
back of your vehicle must be even with the bottom of the ramp (which means that the front 
of your car cannot be more than 1 foot up the ramp). 

5. The vehicle must be started by one single activation device (e.g., a switch or mechanical 
release) on the vehicle. Team members may not activate any device before the start of the 
contest and the vehicles may not have their “motors” running before the start. A team 
member will place the vehicle on the ramp and an “official,” not a team member, will 
start the vehicle, so you must have a simple way to activate the vehicle. 

6. The vehicle can be powered by the following energy sources only, either individually or 
combined (but still only one single activation device is allowed for all energy storage devices 
combined, and this activation device can only be used at the start of the contest). Note that 
these energy sources can be used for propulsion, offensive or defensive purposes. 

• Batteries — limited to 2 AA batteries. Experience has shown that some groups have lost 
contests because their batteries ran out at the most inopportune moment. Use fresh batteries at the 
start of the contest and bring extras. 

• Mousetraps having a maximum spring size of 5 mm outside diameter x 5 cm long 
(use as many as you want). This basically corresponds to a standard mouse trap; the 
specifications preclude the use of rat traps. 

• Rubber bands having a maximum size of 4 mm width x 10 cm length in their 
unstretched state (use as many as you want). 

7. While it shouldn’t be necessary to say so in view of rule 6, just for completeness here it is. 
The vehicle may not use fuels , explosives, compressed gases or any dangerous / hazardous materials. Also, if 
batteries are used, they may not be used for a “thermal protection device” (i.e. a device that 
attempts to set the opposing vehicle on fire.) If in doubt, ask PDR what is/is not acceptable. 

8. Prior to the start of the competition, vehicles will be measured and weighed to ensure 
compliance with size and mass limits. If the vehicle is modified in any way during the course 
of the competition, it will be reweighed and remeasured. 

9. The vehicle must run within the 1-foot-wide track. The vehicle may not run on top of the 
guide rails, but parts of the vehicle may hang over the guardrail. 

The contest 

1. Vehicles will compete in a 5-round “Swiss system tournament” similar to that typically used 
in chess contests (see for example http://en.wikipedia.org/wiki/Swiss system tournament) 
so that every team will participate in the same number of rounds. Modifications to the 
vehicle are permitted between (but obviously not during) runs. 

2. The contest will be held on October 6 at 8:30 AM (i.e., during class time) in the E-quad. 
Since many of you have classes during the time block (8:00 am — 9:30 am or 9:30 am — 11:00 
am) when you don’t have AME 101, your group may have to staff the contest in shifts. 
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3. Teams may inspect the ramps before the contest and ask for modifications if they believe the 
ramps are not level, have undesirable obstacles, etc. 

4. The order of the 1-on-l contests will be selected at random and published in advance. If 
you are not ready at the time of your event, your team will be disqualified for that specific 1- 
on-1 contest. Rigid adherence to schedule is necessary to conduct all 55 or so 
contests during the available time. After the previous contest is completed, you will 
have 1.5 minutes to ready your vehicle for the event. 

5. To ensure fairness, the vehicles will be started by an “impartial” judge, not a team member. 

6. Vehicles will be given 30 seconds to climb to the top of the hill and compete for possession 
of the hill. They will start at the bottom of the ramp, on the sloped part, rather than on the 
flat ground at the base of the ramp. It is possible that a vehicle may not make it all the way 
to the top or may travel across the top to the other side. At the end of the 30 second 
period, the vehicle whose farthest point is closest (in 3 dimensions) to the center of 
the top mesa will be declared the winner of that particular trial. Any parts such as 
telescoping probes, anchors, oil slicks or projectiles count as part of the vehicle. I will make 
exceptions if in my judgment a part unintentionally breaks off, that won’t count as part of 
the vehicle. 

7. If at the end of the contest one vehicle is off the ramp and another is on the ramp, then the 
vehicle on the ramp wins, even if the other vehicle is closer to the top mesa. 

8. If at the end of the contest the judges cannot declare a clear winner (because the difference 
in distances if very small), the contest will be repeated. 

9. It is entirely acceptable (in fact it’s really the whole point of the contest) to incorporate 
means to remove the opposing vehicle from the top of the hill, as well as defensive measures 
to prevent your vehicle from being removed from the top of the hill and/or being damaged 
by the opposing vehicle. 

10. Project reports will be due on Monday, October 21, at 4:30 pm in OHE 430N. 
Organization and grading 

You will work in teams of 4 people, assigned at random, nominally with 2 people from each 
section of the class. Each team must keep a report of their work. The report will be the primary 
means of grading the projects (2/3). Your level of success in the competition also counts (1/3) 
according to the formula: 

Competition score = 50 * (N-P)/(N-l) +50 

where P is your place at the end of the Swiss tournament and N is the total number of teams 
competing. So even the last place team (N = P) gets 50% on the competition score. Also, the 
report/documentation score counts more than the competition score; in other words more weight is 
given to the process than the result, which is just the opposite of how life really works where you’re 
judged 100% on the result and 0% on the process. 
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Spaghetti Bridge 


The goal of this project is to design, construct and test (to failure) a bridge made out of spaghetti. 

Your score on the performance part of the project is based entirely on how much weight your 

bridge holds before failure. The testing will be conducted in the E-quad during class time on 

November 25, 2014. 

Rules: 

1. The only construction materials allows are (1) spaghetti (any brand, uncooked or cooked) 
and (2) Elmer's white school glue. No other materials are allowed. In particular, if you 
use any glue other than white glue (Elmer’s or generic equivalent) your bridge will 
be disqualified. The judges will keep a sample of your (broken) bridge after the contest to 
confirm adherence to the rules. 

2. Your bridge should weigh no more than V 2 pound. There is no bonus for having a bridge 
lighter than V 2 pound but if your bridge is overweight, the “official” weight your bridge held 
will be reduced according to the formula: 

Weight (official) = Weight (actual) x (0.5 Ibf / Bridge weight) 2 

In other words, you will be penalized not just linearly, but in proportion to the square of 
how much you’re overweight. 

3. Spaghetti is defined as a generally circular cross-section pasta not more than 2 mm in 
diameter. If in doubt, show me your pasta before building your bridge. (Obviously ziti or 
some hollow structure would have a better strength to weight ratio, but you can laminate the 
spaghetti into whatever shapes you want.) If you want, you can use cooked spaghetti to 
build a suspension bridge with spaghetti “cables.” (Of course the weight of the water in the 
cooked spaghetti is counted in your total bridge weight.) 



Figure 26. Diagram of spaghetti bridge contest 
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4. The bridge must be able to span a gap of 25 inches between two wooden dowels (cylinders), 
each 1 inch in diameter. The bridge cannot have any structure below the plane of the top of 
supporting dowels, nor can anything be wrapped around the dowels. Loads will be applied 
by suspending a bucket below the center of the bridge from another 1 inch diameter dowel 
and slowly filling the bucket with sand until the bridge fails. 

You will work in teams of 4 ± 1 people, assigned at random. Each team will keep a report of their 
work. Keep in mind that “real” engineered systems always start with a statement of specifications, 
followed by brainstorming, then analysis, testing and construction, followed by an evaluation of the 
performance of the final product. So your report should reflect “real” engineering practice. 

The report will be the primary means of grading the projects (2/3). Your level of success in the 
competition also counts (1/3). The performance score (out of 100 possible) will be computed as 
follows: 

Performance score = 50(1 + W/M) 

where M = maximum weight held by any bridge and W = weight held by your bridge, so your 
group’s score will be between 50 and 100. So more weight (no pun intended) is given to the process 
(i.e. your report) than the result, which is just the opposite of how life really works. The report 
should include: 

• Drawings of preliminary design concepts and critiques of these designs - In the report you 
need to explain your design in terms of maximizing the strength for the given maximum 
weight. So talk about what you did that increased strength more than anything else you could 
have done. 

• Results of analysis using SolidWorks 

• Test data for preliminary designs. At a minimum you should test different brands of spaghetti 
for their strength to weight ratio (NOT just strength), and test different ways of gluing the 
pieces together and laminating multiple strands together. Simply saying that you did some 
testing without showing the results isn’t very enlightening. 

• Reports with no testing and no analysis will get lower scores for the report, regardless 
of how well the bridge did in competition. 

• Construction techniques 

• Results of the “official” test 

• “Post-mortem” of your bridge (why it failed). If you used Solidworks or a bridge design 
program, did the bridge fail at the location where the predicted stress was highest? 

• What you would do differently if you built another bridge 

• Meetings — agenda, minutes, action items 

• Copies of e-mail exchanges 

• Whatever else you think is appropriate — pictures of your test apparatus, construction 
“shop” and the final bridge are very nice to have. Some people show pictures of “real” 
bridges or other people’s bridges, which isn’t very insightful. Also, I prefer to have the 
pictures in the body of the report where they are referred to in the text, not at the end where I 
have to keep referring back and forth. 

• Project reports will be due on Friday, December 5, 2014 at 4:00 pm in OHE 430N. 
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Please no binders for these reports, just pages of plain paper stapled together; 20+ 3-ring 
binders take up too much space! 

Suggestions for the design/construction 

• Do some strength testing of various brands of spaghetti and pick the best one before 
building 

• Keep in mind that spaghetti may be much stronger in compression than tension or vice 
versa — only testing will determine if this is true or not. This means that the optimal design of a 
spaghetti bridge could be very different from a bridge made of steel, which has similar strength in tension and 
compression. Overall, there needs to be as much compressive as tensile force on the bridge 
(otherwise it would be moving!) so that if you find that spaghetti is stronger in compression 
then members in tension need to be thicker than those in compression or vice versa. But 
also keep in mind that members in compression can buckle, so they should have side 
supports. 

• Make sure your bridge is a little longer than 25 inches so that when the load is applied and it 
starts to bend, it won’t slip between the supports! 

• Construction tolerances seem to be the most frequently overlooked aspect. The bridges 
never come out as straight and well balanced as one would hope due to the lack of 
uniformity and straightness of the spaghetti as well as warping by the glue, which greatly 
softens the spaghetti before it dries. Consider using a “scaffolding” to hold the bridge in 
place while it dries. If the bridge is warped, when the load is applied one side will fail much 
sooner than the other side would have. 

• Complicated designs with many individual pieces of spaghetti glued together (e.g. in 
tetrahedral patterns) should be able to hold more weight, but in practice it seems that the 
difficulty in construction outweighs the potential advantages. 

• Many people find that spaghetti is more brittle after being glued. (And everyone finds that it 
takes a LONG time for the glue to dry.) So you may want to minimize the amount of glue 
you use, and when you laminate several strands of spaghetti together to form a beam, you 
might want to “spot glue” them at intervals on the beam rather than gluing the entire beam. 
Again, only testing will determine what works best for you. 

• You may find it useful to build a simple scale-model bridge to test out your ideas before 
building an elaborate, full-scale bridge. You may also want to build a back-up full scale 
bridge in case a problem (e.g. last minute breakage of your primary bridge.) 

• If you do use SolidWorks or another program to help design your bridge (highly 
recommended), you need to explain how you used the results to improve your design. 

• You can also construct a virtual bridge using SolidWorks and test it using the stress analysis 
feature to estimate its breaking point. In particular, if the program shows that some 
elements are under more stress than others, you should strengthen that element, for example 
by increasing the number of strands of spaghetti in that element. Conversely, elements with 
little stress should be weakened or eliminated completely. Show how you use the results 
to maximize your strength to weight ratio. Just running the program without using 
the results isn’t very insightful. What is really cool and will get you a really good 
grade on the report is if you actually use SolidWorks or another program to predict 
the load at which the bridge will fail and compare the actual failure load to the 
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predicted one. I don’t expect that the agreement will be very good (but sometimes it is, 
surprisingly) and having a poor agreement won’t hurt your report score, but if the agreement 
isn’t good you should list some possible reasons for the discrepancy (i.e., what happens in 
reality that isn’t modeled by SolidWorks?) 

If you want a really awesome grade on the report, you’ll construct a 3D SolidWorks model 
and intentionally make it not quite symmetrical (just like your real bridge will be). That way, 
any twisting of the bridge as you load it can be predicted and you can add cross-braces to 
avoid twisting which will otherwise greatly reduce your bridge’s strength. 

If you do the modeling you’ll find this out anyway, but inevitably a bridge with a 
superstructure (i.e. a truss bridge) will be stronger than a flat bridge for the same weight. 
This is a natural consequence of the fact that you’ll get more moment of inertia (I) for the 
same weight with a tall structure than a flat one (think about the discussion of I-beams). A 
flat bridge is highly un-recommended. Also, you need similar amounts of material at the 
roadbed level and at the top of the structure in order to have comparable stresses in tension 
and compression. Adding a few strands of superstructure isn’t going to help much. 

As simple as this sounds, transportation is a problem. Every year several bridges fail 
before the contest because it broke during transportation from the dorm room to the 
contest site. Be careful! 
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Plaster of Paris Bridge 

The goal of this project is to design, construct and test (to failure) a bridge made out of Plaster of 
Paris. Your score on the performance part of the project is based entirely on how much weight your 
bridge holds before failure. The testing will be conducted in the E-quad during class time on 

November 26, 2013. 

Rules: 

i. 

1. The only construction material allowed is Plaster of Paris. Since there are a lot of products 
with “plaster” in the name, you material must have the words “Plaster” and “Paris” on the 
package. You can get a 25 lbm bag at Home Depot or Walmart for about $16, e.g.: 

http://www.walmart.com/ip/Dap-10312-25-Lb-Plaster-of-Paris-Exterior/21798112. 

No other materials are allowed. If you use any other material your bridge will be 
disqualified. The judges will keep a sample of your (broken) bridge after the contest to 
confirm adherence to the rules. 

2. Your bridge must weigh no more than 2 pounds. There is no bonus for having a bridge 
lighter than 2 pounds but if your bridge is overweight, the “official” weight your bridge held 
will be reduced according to the formula: 

Weight (official) = Weight (actual) x (2 Ibf / Bridge weight) 2 

In other words, you will be penalized not just linearly, but in proportion to the square of 
how much you’re overweight. 

3. The bridge must span a gap of 25 inches between two tables. The bridge cannot have any 
structure below the plane of the tables, nor can anything be attached to the tables. 

4. Loads will be applied by suspending a bucket below the center of the bridge from a 1 inch 
diameter dowel and slowly filling the bucket with sand until the bridge fails. You can ask the 
contest officials to place the dowel anywhere near the center of the bridge. You should 
make a “cradle” for the dowel so it doesn’t slip or roll away from the center. 
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You will work in teams of 3 or 4 people, assigned at random. Each team will keep a report of their 
work. Keep in mind that “real” engineered systems always start with a statement of specifications, 
followed by brainstorming, then analysis, testing and construction, followed by an evaluation of the 
performance of the final product. So your report should reflect “real” engineering practice. 

The report will be the primary means of grading the projects (2/3). Your level of success in the 
competition also counts (1/3). The performance score (out of 100 possible) will be computed as 
follows: 

Performance score = 50(1 + W/M) 

where M = maximum weight held by any bridge and W = weight held by your bridge, so your 
group’s score will be between 50 and 100. So more weight (no pun intended) is given to the process 
(i.e. your report) than the result, which is just the opposite of how life really works. The report 
should include: 

• Drawings of preliminary design concepts and critiques of these designs - In the report you 
need to explain your design in terms of maximizing the strength for the given maximum 
weight. So talk about what you did that increased strength more than anything else you could 
have done. 

• Results of analysis using SolidWorks or other software 

• Test data for preliminary designs. At a minimum you determine the strength to weight ratio 
(NOT just strength) for different ways of preparing the Plaster of Paris (i.e. different amounts 
of water) and different ways of assembling individual pieces together (unless you cast your 
whole bridge in one piece). Simply saying that you did some testing without showing the 
results isn’t very enlightening. 

• Reports with no testing and no analysis will get lower scores for the report, regardless 
of how well the bridge did in competition. 

• Construction techniques 

• Results of the “official” test 

• “Post-mortem” of your bridge (why it failed). If you used SolidWorks, did the bridge fail at 
the location where the predicted stress was highest? 

• What you would do differently if you built another bridge 

• Meetings — agenda, minutes, action items 

• Copies of e-mail exchanges 

• Whatever else you think is appropriate — pictures of your test apparatus, construction 
“shop” and the final bridge are very nice to have. Some people show pictures of “real” 
bridges or other people’s bridges, which isn’t very insightful. Also, I prefer to have the 
pictures in the body of the report where they are referred to in the text, not at the end where I 
have to keep referring back and forth. 

• Project reports will be due on Tuesday, December 10, at 12:00 noon in OHE 430N. 

Please no binders for these reports, just pages of plain paper stapled together; 20+ 3-ring 
binders take up too much space! 

Suggestions for the design/construction 

• Do some strength testing and pick the best construction techniques before building your 
“real” bridge. 
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Keep in mind that Plaster of Paris may be much stronger in compression than tension or 
vice versa — only testing will determine if this is true or not. This means that the optimal design of 
a Plaster of Paris bridge could be very different from a bridge made of steely which has similar strength in 
tension and compression. Overall, there needs to be as much compressive as tensile force on the 
bridge (otherwise it would be moving!) so that if you find that Plaster of Paris is stronger in 
compression then members in tension need to be thicker than those in compression or vice 
versa. But also keep in mind that members in compression can buckle, so they should have 
side supports. 

Make sure your bridge is a little longer than 25 inches so that when the load is applied and it 
starts to bend, it won’t slip between the supports! 

Construction tolerances seem to be the most frequently overlooked aspect. The bridges 
never come out as straight and well balanced as one would hope. Consider using a 
“scaffolding” to hold the bridge in place while it dries. If the bridge is warped, when the 
load is applied one side will fail much sooner than the other side would have. 

Complicated designs with many individual pieces (e.g. in tetrahedral patterns) should be able 
to hold more weight, but in practice it seems that the difficulty in construction outweighs the 
potential advantages. 

You may find it useful to build a simple scale-model bridge to test out your ideas before 
building an elaborate, full-scale bridge. You may also want to build a back-up full scale 
bridge in case a problem (e.g. last minute breakage of your primary bridge.) 

You can construct a virtual bridge using SolidWorks and test it using the stress analysis 
feature to estimate its breaking point. If you use this approach, show how you used to to 
improve/optimize the design. What is really cool and will get you a really good grade 
on the report is if you actually use one of these programs to predict the load at which 
the bridge will fail and compare the actual failure load to the predicted one. I don’t 
expect that the agreement will be very good (but sometimes it is, surprisingly) and having a 
poor agreement won’t hurt your report score, but if the agreement isn’t good you should list 
some possible reasons for the discrepancy (i.e., what happens in reality that isn’t modeled by 
SolidWorks?) 

If you want a really awesome grade on the report, you’ll construct a 3D SolidWorks model 
and intentionally make it not quite symmetrical (just like your real bridge will be). That way, 
any unsymmetrical loading of the bridge as you load it can be predicted and you can add 
cross-braces to avoid twisting which will otherwise greatly reduce your bridge’s strength. 

If you do the modeling you’ll find this out anyway, but inevitably a bridge with a 
superstructure (i.e. a truss bridge) will be stronger than a flat bridge for the same weight. 
This is a natural consequence of the fact that you’ll get more moment of inertia (I) for the 
same weight with a tall structure than a flat one (think about the discussion of I-beams). A 
flat bridge is highly un-recommended. Also, you need similar amounts of material at the 
roadbed level and at the top of the structure in order to have comparable stresses in tension 
and compression. Adding a thin superstructure isn’t going to help much. 

As simple as this sounds, transportation is a problem. Every year several bridges fail 
before the contest because they broke during transportation from the dorm room to the 
contest site. Be careful! 
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Hydro power 


The purpose of the second design project is to construct a “gadget” that quickly and efficiendy 
converts the gravitational potential energy of water into mechanical power. The power generated 
will be determined by the time (t) required to raise a 1.5 kilogram mass a distance of 4 meters, 
starting at rest. The volume of water (V) needed to raise the mass, which is a measure of the 
efficiency of energy conversion, will also be measured. The total performance score will be 
computed as follows: 

t . V • 

Score = 50^+50^^ 

t V 

where t min and V min are the minimum time and volume among all the teams. 

1. The “Kit.” Each team will receive a Turbine Kit. Any or all of the parts in this Kit may be 
used. (See Figure 28 for supplied components.) At a minimum, one of the supplied turbine 
wheels must be used. You may use 2 or even all 3 turbine wheels if you want, but only these 3 
and no others. No modifications can be made to any of the parts of the Kit (so we can re-use 
them in future years). Pick up your parts from Sylvana in RRB 101. 





Figure 28. Turbine kit containing 3 turbines, 2 nobles, 

1 hose coupling and 1 connector section 

2. Other stuff you can/must use. Other than turbine wheels, you can add any parts you want to 
the gadget; in particular you may want to build an enclosure for your turbine wheel(s). The 
gadget may use shafts, gears, pulleys or other mechanisms to convert the (rotational) turbine 
power into (linear) mechanical movement. Each team will have a $50 budget from McMaster- 
Carr (http://www.mcmaster.com) to buy such parts. Use the “ work order ” below to purchase 
parts; fill out the form and give it to Sylvana in RRB 101. 

3. No energy storage allowed! All motive force produced by the gadget for lifting the weight 
must come from the potential and kinetic energy of the water that flows through the gadget. 
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No means of storing mechanical, electrical, etc. energy in the gadget or take-up line prior to its 
operation by water is allowed. 

4. Mounting your gadget to a board compatible with the test stand. The gadget must be 
attached to a mounting board ( Figure 29) (provided by AME) that is 20 inches wide and 15 
inches in height. The mounting board will have two 1/2 inch holes drilled through it. The 
holes will be centered at a point 1-1/2 inches from the edges of the board nearest the two upper 
corners. These holes will be used to attach your board to support brackets on the test fixture 
(see Figure 30). No part of your gadget can protrude beyond the 20 inch x 15 inch 
envelope of the board. You can only use one side of the board. 


Mounting 

hole 



Mounting 

hole 


Hose 
connector 
mounts here 


Figure 29. Picture of mounting board, with some turbine wheels shown for scale. 


5. Plumbing. The gadget will be connected to the test apparatus using a standard garden hose 
thread female connector (on the gadget) that connects to a male garden hose connector (on the 
test fixture) (see Figure 30). The male garden hose can be attached to your gadget anywhere on 
the board. 

6. The testing procedure. Up to 30 liters of water will be allowed to flow, starting at a 4-meter 
elevation above the control valve centerline, through a conduit (a 1-inch ID (inside diameter) x 
1.5-meter-long polyethylene tube) to a flow meter. From the flow meter, the water will run 
through another conduit (a 1-inch ID x 2.5-meter-long polyethylene tube) to the control valve, 
and then to your gadget. Your gadget will lift the 1.5 kg mass vertically 4 meters. A sump 
beneath the testing apparatus will collect water flowing out of the gadget; you don’t need to do 
“water management.” Figure 30 is a sketch of the test stand. 
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Figure 30. Schematic of test stand. Note: some of the dimensions are wrong in this figure; see item 6 above for the 

correct dimensions. 

To “seed” your brainstorming process, Figure 31 shows just about the simplest possible arrangement 
for a turbine-powered lifting gadget. Water strikes a turbine wheel that is fixed to the same shaft as 
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the take-up spool. As the turbine rotates due to the force of the water, the spool also turns thereby 
winding the take-up line. 



Figure 31 . Simple arrangement for turbine-powered lifting gadget. 


Teaming. You can work in teams of between 2 and 4 people, either the same group of people as 
the first design project or a different group. Each team must keep a report of their work. The 
report will be the primary means of grading the projects (2/3). Your score in the competition also 
counts (1/3). So more weight (no pun intended) is given to the process than the result, which is just 
the opposite of how life really works. The report should include: 

• Meetings — agenda, minutes, action items (see below...) 

• Copies of e-mail exchanges 

• Drawings of preliminary design concepts and critiques of these designs 

• Test data for preliminary designs 

• Results of the “official” test 

• “Post-mortem” - what you would do differently if you built another hydropower system 

• Whatever else you think is appropriate 
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Appendix C. Problem-solving methodology 


1. Make a clean start (clean sheet of paper) 

2. Draw a picture 

3. State givens, state unknowns (in real life, in most cases you won’t have enough givens to 
determine the unknowns so you’ll have to turn some of the unknowns into givens by making 
“plausible” assumptions.) You need to have as many equations as unknowns. If you 
can state what the equations and the unknowns are, you’re 90% of the way to the 
solution. 

4. Think, then write (I don’t necessarily agree with this...) But it is essential that you STATE 
YOUR ASSUMPTIONS AND WRITE DOWN THE EQUATIONS THAT YOU USE 
BEFORE YOU PLUG NUMBERS INTO SAID EQUATIONS. Why is this so 
important? To make sure that your equations are valid for the problem assumptions. For 
example, Bernoulli’s equation is valid only for steady, incompressible (constant-density) flow 
but on the panic of an exam you’ll try to apply it to a gas at high Mach number. Another 
example is Hooke’s law, which applies only to an elastic material, not Play-Doh. 

5. Be coordinated — show your coordinate system on you picture and follow through with this 
coordinate system in your equations 

6. Neatness counts 

7. Units 

8. Significant figures 

9. Box your answer 

10. Interpret the result — is it reasonable? 
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Appendix D. Excel tutorial 


Excel is a fairly powerful tool for data analysis and computation. It is primarily oriented toward 
financial calculations but it works reasonably well for engineering and scientific calculations as well. 
It is certainly not as capable as programs like Matlab, Mathematica, TKSolver, etc. but most people 
have free access to Excel. And you can imbed an Excel spreadsheet in a Word or Powerpoint 
document such that all you need to do to open the spreadsheet is to click on the figure. Plus, a lot 
of people (PDR included) have created spreadsheets in Excel for solving various types of problems 
that can be downloaded from the internet. Some of my spreadsheets are available on-line at 

http://ronney.usc.edu/excel-spreadsheets/ 

This short tutorial is aimed to give you a few pointers at how to use Excel for engineering problems. 
Of course, there’s no substitute for actually playing with the program — reading about how to use 
software is about as useful as reading about how to ride a bicycle. 

Cells, rows and columns 

Spreadsheets are organized into cells arranged in rows and columns of information. In Excel the 
columns are A, B, C, ... and the rows are 1, 2, 3, — So the address of the cell in 4 th column, 5 th row 
would be D5. 

Formulas 

Each cell may contain raw data (i.e. just a plain number) or a formula, or text. The formula cells 
generally refer to other cells, for example if cell A1 had a number in it, and B1 had another number, 
and you wanted to know the sum of those 2 values, you could enter =A1+B1 into another cell (not 
A1 or Bl) for example Cl: (If you’re viewing the Word version of this document (not the online 
web version, not the pdf version), you can double click the table to open the Excel spreadsheet) 


1st number 

2nd number 

Sum 

7 

5 

12 


You can look at the formula entered into a cell by clicking on that cell and looking at the “formula 
bar” at the top of the screen. There’s a zillion different functions you can use in Excel, e.g. addition 
(+), multiplication (*), subtraction (-), division(/), exponentiation ( /x ), ln(), exp(), sin(), etc. Pick 
“function” from the “insert” menu to see the available functions. Some functions like SUM, 
AVERAGE, STDEV, etc. refer to an array of cells rather than an individual cell, in which case the 
formula is of the form =SUM(A1 :B10). (Note that the array of cells can be a vertical column, a 
horizontal row, or a block more than one cell wide in both the horizontal and vertical directions). 

Also, sometimes you want to create a formula in one cell then copy/paste the same formula into 
other cells, e.g. E = me 2 . If your cell contains a constant, when you copy/paste, you’ll get the 
constant in all the cells into which you paste. If your cell contains a formula, when you copy/paste, 
you’ll get that formula in the other cells, but the cells to which the formula refers will be adjusted 
accordingly. For example, if cell Cl contains the formula =A1+B1, if you copy/paste this formula 
into cell E7 (2 columns to the right and 6 rows down), the formula in cell E7 will read =C7+D7 (each 
cell reference is changed by 2 columns to the right and 6 rows down). This is extremely convenient 
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for calculating E = me" for a large set of masses (m), but really you only want to enter c (speed of 
light) in one cell, and have all formulas refer to the value of c in that cell which would be (in the 
example below) cell A2. In that case you can use an “absolute reference” to that cell which is of the 
form $A$2 rather than just A2. Without the dollar signs formulas use “relative references” and thus 
change when you copy/paste them. (If on the other hand you cut and paste rather than copy and 
paste, meaning you’re moving a cell or cells from one location to another rather than creating new 
formulas in new cells, then the referencing doesn’t change, that is, A2 stays A2.) 


c (m/s) 


3.00E+08 

m (kg) 

1 

10 

100 

1000 

E (Joules) 

3.00E+08 

3.00E+09 

3.00E+10 

3.00E+11 


There are also array formulas that are especially useful for solving a set of simultaneous linear 
equations. They are rather cryptic to create, almost like a secret handshake, so I’ll just give you a 
“template” you can use: 


1.01E+00 

2 

3 

4 

10 

1 

2 

3 

4 

10 

1 

-1 

1 

-1 

0 

4 

5 

2 

0.00E+00 

12 

-1.003E-14 

1.64102564 

1.8974359 

0.25641026 



This solves the set of equations 
1 X t + 2 X 2 + 3 X 3 + 4 X 4 = 10 
8 X, + 7 X 2 + 6 X 3 + 5 X 4 = 26 

ix r ix 2 + ix 3 + -ix 4 = o 

4 X , + 5 X 2 + 2 X 3 + 1 X 4 = 12 

which has the solution X 1 = X 2 = X 3 = X 4 = 1. 

Another useful function is “Goal Seek” from the “Tools” menu, for which you can ask Excel to 
modify the value in one cell until another cell has a specific value. For example, you could input the 
formula for the left-hand side of an equation in one cell, input the formula for the right-hand side of 
the equation into another cell, then set another cell to compute the difference between the right and 
left-hand sides, and use Goal Seek to find the solution. Let’s suppose you want to find x such that 
150 sin(x) e x = 12 ln(x) + 7x" (x in radians for the sin(x) term). There’s no way to solve this 
analytically, so you have to do tedious trial and error to find the solution, so set up the spreadsheet: 


X 

150*sin(x)*exp(x) 

12*ln(x) + 7x A 2 

LHS-RHS 

1.88396438 

32.44579276 

32.44579189 

-0.000875599 


In this case use Tools / Goal Seek / Set cell: D2 / To value: 0 / By changing cell: A2 

Of course, your initial guess of x has to be good enough that Excel can converge on the solution. 
There are also other selections under the Tools menu such as “Solver” that has more options (like 
changing multiple cells to find the solution, or find the maximum or minimum rather than a specific 
value, optionally subject to constraints such as certain cells have to be greater than zero) but in my 
experience Solver is less reliable for simple problems — use Goal Seek if it will do what you need. 
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A very powerful “dirty trick” within Excel is the “iterate” feature. Select Preferences / Calculation 
and check “Manual.” Also select “Iteration” and set “Maximum iterations” = 1. With this, Excel 
(a) does not update the calculations automatically, but only when you type Cmd = (on the Mac, or 
something similar on the PC) and (b) Excel doesn’t complain when cells refer to each other (circular 
references, like if you had the formula “=A1” in cell Bl, and “=A1” in cell Bl. This might not seem 
like anything useful, but in most scientific calculations, one has a large set of simultaneous, non¬ 
linear equations and the only way to solve them is iteratively. Each time you type Cmd =, the 
calculation advances by one iteration towards the solution. A trivial example of this is to put the 
formula “=A1 + 1” in cell Al: 



Every time you hit Cmd -, the value of this cell will increase by 1. 

Also, for time-dependent problems, you can use the iteration feature and each iteration will 
increment the solution by one time step. I have written a fairly elaborate sheet for use in heat 
conduction problems: 

http://ronney.usc.edu/spreadsheets/Unsteady 2D conduction.xls 

You can also plot data sets by highlighting the cells and selecting “Chart” from the “Insert” menu 
and you’ll get a bunch of options of what to plot and how to plot it. Excel doesn’t make very good 
quality plots suitable for publication in journals, but they’re adequate for homework, internal reports, 
etc. If you click on the chart you created and select the “Add Trendline” option from the Chart 
menu, you can add a least-squares fit to the data in the form of a line, polynomial, power law, etc. 

One of my favorite examples of a fairly complete spreadsheet package including plotting is the one I 
wrote for analyzing internal combustion engine cycles including the effects of compression and 
expansion, heat losses, the rate of combustion, the exhaust gases trapped in the cylinder after the 
end of the exhaust stroke, etc.: 

http://ronney.usc.edu/spreadsheets/AirCycles4Recips.xls 
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